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I.  Overview  -  L.F.  Eastman 

Both  experimental  and  theoretical  studies  have  been  made  on  the  use  of 
AlAs/GaAs  short  period  superlattices  to  suppress  phonon-electron 
scattering,  for  improved  electron  transport  parallel  to  the  heterojunctions. 
The  structures  would  be  applicable  to  MODFET  devices.  These  two 
efforts,  experimental  and  theoretical,  are  reported  separately  in  the 
following  two  sections.  In  spite  of  the  confmement  of  the  GaAs  polar 
optical  phonons  in  the  short  period  superlattice  channels,  the  electrons  had 
lower  mobility  at  room  temperature  than  they  would  have  in  a  single,  thick 
GaAs  channel.  The  added  scattering  by  phonons  due  to  the  added  abrupt 
heterojunction  interfaces  was  determined  to  be  the  cause  of  the  reduced 
mobility.  The  abrupt  interfaces  that  occur,  even  with  a  single,  uniform 
composition  channel,  also  limit  the  room-temperature  mobility.  In  order 
to  concentrate  the  electrons  in  the  region  of  the  channel  away  from  the 
interfaces,  graded  InGaAs,  with  a  higher  Indium  fraction  in  the  center  of 
the  channel  was  studied.  This  yielded  10%  higher  mobility  at  room 
temperature,  and  MODFET's  fabricated  from  the  graded-channel  material 
showed  other  important  improvements  in  breakdown  voltage  and  gate 
leakage.  A  reprint  by  a  visitor  is  attached  to  cover  this  important  new 
concept,  for  which  an  invention  disclosure  has  been  filed  at  Cornell.  A 
new  project  on  the  graded-composition  channel  has  been  initiated  with 
ARO  and  ARL  support. 


II.  On  the  Suppression  of  Phonon-Electron  Scattering  in  Short  Periodic 

AlAs/GaAs  Multiple  Quantum  Well  Structures  -  T.  T.  LeTran,  W.J. 

Schaff,  B.K.  Ridley,  Y.P.  Chen,  A.  Clark,  S.  O'Keefe,  and  L.  F. 

Eastman 

2.1  Abstract 

The  suppression  of  longitudinal  optical  phonon  (LOP)-electron 
scattering  was  sought  in  multiple  quantum  well  (MQW)  structures.  The 
structures  had  GaAs  well  widths  -  12,  15  and  20  mono-layers  (ML)  and 
AlAs  barrier  widths  =  2  and  4  ML.  The  MQWs  were  grown  in  the  channel 
of  GaAs/Alo.sGao.vAs  modulation  doped  field  effect  transistors 
(MODFETs)  without  gates.  The  Hall  mobility  and  carrier  sheet  density 
were  measured  by  the  van  der  Pauw  method.  The  Hall  mobility  of  the 
MQW  samples  were  found  to  be  less  than  that  of  the  control  samples 
(without  MQW)  at  room  temperature,  but  was  better  at  temperatures  lower 
than  50  K.  The  reduction  of  the  room  temperature  mobility  was  due  to 
interaction  of  the  well  electrons  with  the  interface  polaritions  from  the 
AlAs  barriers.  The  increase  of  the  low  temperature  mobility  was  due  to 
reduced  remote  ionized  impurity  scattering  of  the  well  electrons.  The 
evidence  of  performance  improvement  of  MQW  devices  at  room 
temperature  due  to  suppression  of  electron-LOP  scattering  is  thus  disputed 
by  this  study.  The  results  of  an  experiment  made  elsewhere,  which 
appeared  to  show  the  contrary,  can  be  interpreted  using  arguments  other 
than  the  suppression  of  electron-LOP  scattering  in  MWQs. 

2.2  Introduction 

In  recent  years,  the  possibility  of  electron  mobility  enhancement  by 
suppression  of  phonon-electron  interaction  in  multi-quantum  well  (MQW) 
structures  has  raised  considerable  interest  1-5.  The  suppression  is  sought 
by  the  confinement  of  longitudinal  optical  phonons  (LOP)  inside  the 


quantum  wells  that  have  different  coupling  constant  and  ion  masses  than  the 
barriers.  The  interaction  of  2D-LOP/2D-eiectrons  gets  weaker  as  the  well 
width  decreases,  in  contrast  to  the  3D-LOP/2D-electron  system.2  On  the 
other  hand,  interface  polariton  (IP)  interaction  with  free  carriers  becomes 
more  important  as  the  well  width  decreases  in  systems  of  finite  barriers.3-5 

A  study  of  several  different  GaAs/AlAs  multi-quantum  well  samples  is 
reported  here.  It  indicates  that  there  may  be  an  upper  limit  of  electron 
mobility  parallel  to  the  AlAs/GaAs  MQW  plane  due  to  the  interplay 
between  LOP  and  IP  scattering  of  electrons.  This  study  confirms  the 
importance  of  IP  scattering  in  very  narrow  wells  ( 1 2  to  20  atomic  layers). 
The  latter  is  shown  to  outweight  the  benefit  of  suppression  of  LOP 
scattering. 

2.3  Samples 

Two  control  samples  and  MQW  stmctures  of  different  well  and  barrier 
widths  were  examined.  Fig.  2.1  shows  the  first  control  sample’s  (REF) 
structure.  The  2000A  n-GaAs  cap  layer,  the  SOOA  Alo.3Gao.7As  uniform 
donor  layer,  the  150 A  spacer  layers  and  the  GaAs  channel  were  grown  by 
MBE  on  (100)  GaAs  at  620®C.  The  donor  concentration  was  2x818  cm-3. 
The  second  control  sample,  REF+4,  is  similar  to  the  REF,  except  that  there 
is  a  4  ML  thick  ALAs-barrier  between  the  spacer  and  the  channel.  All 
parameters  in  the  other  samples  were  kept  the  same  as  the  REF's,  except 
that  they  have  AlAs/GaAs  multiple  quantum  well  in  their  channel.  Six  such 
samples  were  studied:  MQW(12/2),  MQW(15/2),  MQW(20/2), 
MQW(12/4),  MQW(15/4),  MQW(20/4),  where  the  numbers  in  brackets 
represent  the  well  and  barrier  widths  in  number  of  ML,  respectively.  Fig. 
2.2  shows  the  band  structure  of  the  MQW(15/4)  sample.  Note  the  large 
density  of  the  low  mobility  electrons  in  the  cap  and  donor  layers. 
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the  bulk  density  of  free  electrons  is  shown  by  the  dotted  line. 


The  barriers  in  the  MQW  samples  were  kept  thin  to  minimize  the 
possibility  of  carriers  being  scattered  into  the  X-valleys  of  AlAs,  while  still 
being  large  enough  to  confine  LOP  to  within  the  wells.  For  example,  the 
energy  gap  between  the  X-level  of  the  AlAs  barrier  and  the  first  F-level  of 
the  GaAs  well  is  0.38  eV  in  the  MQW(15/2)  sample.  The  fraction  of  the 
well's  electrons  being  scattered  into  the  barriers  is  approximately  lO'^.  In 
the  MQW(15/4),  the  fraction  is  just  1%.  The  upper  limit  of  the  well  width 
is  roughly  30  ML,  above  which  the  confinement  of  LOP  may  not  be 
appreciable. 2  If  the  wells  were  narrower  than  10  ML,  the  electronic 
subband  levels  would  rise  above  the  barriers  causing  the  carriers  to  be 
deconfined  by  the  MQW,  as  shown  by  computer  simulations  that  solve  the 
Schrddinger  and  Poisson  equations  self-consistently. 

2.4  Data  and  Discussions 

The  sample's  Hall  mobility  and  carrier  sheet  density  were  measured  by 
the  van  der  Pauw  method.  The  magnetic  field  was  2  kiloGauss.  The 
samples  were  exposed  only  to  diffused  ambient  light  shortly  before  the 
measurements.  To  minimize  the  contribution  of  the  parallel  conduction  in 
the  donor  layer  to  the  total  mobility,  the  samples  were  etched  with  the 
solution  H3P04:H202:H20  =  3:1:50  until  their  individual  Hall  mobility 
reached  a  maximum  value. 

The  mobilities  after  etching  and  the  carrier  sheet  densities  are  shown  in 
Table  2.1.  The  REF's  mobility  was  7590  cm^V-ls-^  at  room  temperature, 
comparable  to  the  value  of  REF+4,  taking  into  account  its  lower  carrier 
sheet  density.  The  MQW  samples,  however,  all  showed  lower  room 
temperature  mobilities  than  the  control  samples.  Among  them,  the 
MQW(15/x)  samples  had  the  highest  mobilities.  At  T  =  77K,  the  mobility 


6 


of  the  MQW  samples  approached  that  of  tl«  control  samples,  except  for  the 
MQW(20x)  samples 


Table  2. 1 .  The  HaU  mobility  and  carrier  sheet  density  of  the  samples  at 

room  and  liquid  nitrogen  temperature.  The  control  sample  REF  has 
the  basic  MODFET  structure;  REF+4  has  an  additional  4  ML  thick 
AlAs  between  the  AlGaAs  spacer  and  the  GaAs  channel.  The 
samples  MQW(x/y)  have  x  ML  of  GaAs  as  well  and  y  ML  of  AlAs  as 
barriers. 


T=300K 

fcm^V’ls-l 

Ns 

T=77K 

^[103 

cm2v-lsM 

Ns 

flO^  ^  cm"2] 

REF 

7590 

5.44 

189 

5.52 

REF+4 

8200 

4.12 

193 

4.15 

MQW(12/2) 

5287 

2.94 

191 

1.80 

MQW(12/4) 

5189 

2.00 

180 

1.48 

MQW(15/2) 

6318 

3.16 

198 

2.20 

MQW(15/4) 

6682 

2.49 

181 

2.66 

MQW(20/2) 

6024 

3.72 

124 

2.48 

MQW(20/4) 

5055 

3.31 

153 

1.90 

The  true  mobility  fo  the  GaAs  channel  (and  its  MQW)  is  given  by  the 
Petriz  model:^ 


i-4^n 

^AlGaAs  ^channel 

where  |i  amd  M  are  the  effective  mobility  and  carrier  sheet  density  at  T  = 
300K,  respectively;  Nchannel  is  the  sheet  density  measured  at  T  <  50K; 
Itchannel  and  Nchannel  are  the  true  mobility  and  sheet  density  at  T  =  300K, 
respectively.  jiAlGaAs  is  the  mobility  of  the  free  electrons  inside  the  donor 
layer  at  300K;  its  value  is  conventionally  taken  to  be  500  to  1000  cm^V'l 
S'T  Thus  lichannel  =  9600  cm^V-ls'l  and  87(X)  cm^V  ls*^  for  the  REF  and 
MQW  (15/x)  samples,  respectively. 

The  lower  mobility  of  the  MQW  samples  is  readily  understood  by  the 
hybrid  model  of  phonons  in  superlattice  structures.  This  model  is 
presented  in  a  paper  submitted  elsewhere.7  In  the  case  of  the  AlAs/GaAs 
MQW  structure,  the  LOP-like  hybrid  modes  of  the  AlAs  barriers  are  not 
coupled  to  the  interface  polaritons  (IP).  The  latter  thus  interact  strongly 
with  the  electrons  in  the  well  region.  Although  the  interaction  of  electrons 
and  LOP  is  reduced  in  the  MQW  samples,  the  increasing  strength  of  the 
electrons-barrier  IP  scattering  enhances  the  overall  scattering  rate  as  the 
well  width  decreases,  and  hence  reduces  the  mobility.  The  channel 
mobility  calculated  by  this  model  was  6120  cm^V-Js’^  for  an  infmite  MQW 
with  15  ML  thick  GaAs  wells  and  2  ML  thick  ALAs  barriers,  compared  to 
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the  actual  value  of  8700  cm^V-ls-l  of  the  (short  periodic)  MQW(15/2) 
sample. 

The  dependence  of  the  Hall  mobility  and  carrier  sheet  density  on 
temperature  was  also  investigated.  To  minimize  the  persistent  photo¬ 
conductivity  effect,  the  samples  were  kept  in  the  dark  for  at  least  10  hours 
before  being  measured  in  a  continuous  flow  cryostat  fed  by  liquid  helium. 
The  results  for  the  REF,  MQW(15/2)  and  MQW(15/4)  are  shown  in  Fig. 
2.3.  One  observes  that  the  mobility  of  the  MQW  samples  raised  above  that 
of  the  REF  at  T  <  50K.  The  MQW(l5/2ys  mobility  was  1.4x10^  cm2  y-l 
s-l,  twice  as  high  as  the  REFs  (6.9  x  105  cm^V-ls-^)  at  T  =  IIK.  The 
sheet  desntiy  of  the  MQW  samples  were  considerably  lower  than  that  of  the 
REF  because  they  were  etched,  as  mentioned  above. 

The  improved  mobility  of  the  MQW  samples  indicates  that  the  electron- 
IP  interaction  is  no  longer  the  dominant  scattering  mechanism  at  low 
temperatures.  Our  low  temperature  simulations  showed  a  smaller  field 
across  the  spacer  layer  in  the  MQW  samples  than  in  the  REF,  which 
reduced  the  effect  of  remote  ionized  impurity  scattering.  More  important, 
the  improved  mobility  shows  that  the  interface  roughness  scattering  is  not  a 
factor  in  lowering  the  MQW  sample’s  mobility  at  room  temperature. 

2.5  Conclusions 

We  have  tested  several  short  periodic  MQW  structures  for  the 
enhancement  of  their  room  temperature  mobility  due  to  the  suppression  of 
electron-LOP  scattering.  The  results  were  interpreted  as  negative.  The 
best  room  temperature  mobility  was  achieved  by  the  MQW  with  15  ML 
thick  GaAs  wells  and  2  or  4  ML  thick  barriers.  It  was  roughly  10%  lower 
than  that  of  the  samples  without  the  MQW.  This  result  was  readily 
explained  by  the  strong  interaction  of  the  well  electrons  with  the  IP-like 
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|t  o(  REF 
H  of  MQW(15/2) 
ji  of  MQW(15/4) 


I 


hybrid  modes  from  the  AlAs  barriers.  At  low  temperatures,  the  MQW 
samples'  mobility  was  higher  than  that  of  the  control  samples,  due  to 
reduced  remote  ionized  impurity  scattering. 

To  the  best  of  our  knowledge,  the  enhancement  of  room  temperature 
mobility  by  suppression  of  electron-LOP  scattering  has  yet  to  be 
experimentally  proven.  The  results  of  Zhu  et.al.^  that  appeared  to  confirm 
the  opposite  may  actually  be  explained  using  some  alternative  arguments. 
The  energy  band  structure  and  the  bulk  carrier  density  of  the  samples  are 
shown  in  Fig.  2.4.  The  MQW’s  wall  and  barrier  thickness  are  8  ML  and 
4ML,  respectively.  These  narrow  wells  raise  the  electron  subbands  over 
the  barriers,  spiUing  most  of  the  electrons  onto  the  AlLAs/GaAs-substrate 
interface,  which  is  both  smooth  he  ionized  donors.  The  major  contribution 
to  the  total  mobility  is  from  the  2DEG  at  this  interface,  not  from  the 
electrons  in  the  MQW  structure.  In  case  of  the  control  sample,  the 
electrons  stay  inside  the  GaAs  well  and  are  subjected  to  more  interface 
roughness  and  ionized  impurity  scattering,  as  compared  to  the  majority  of 
the  electrons  in  the  MQW  sample;  hence  the  mobility  is  considerably 
lower.  Please  not  that  both  of  these  samples  have  lower  room  temperature 
mobility  than  bulk  GaAs. 
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with  its  GaAs  well  at  the  center,  where  the 
right  is  the  AlAs/GaAs-substrate  interface. 
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electrons  are  spilled  into  the  smooth  AlAs/GaAs-substrate  interface. 
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III.  Optical  Mode  Hybrids  in  Quantum  Wells,  Suj)erlattices  and  Slabs  - 
B.K.  Ridley 

3.1  Abstract 

The  presence  of  interfaces  causes  the  hybridization  of  LO  and  TO 
modes  in  non-polar  material  and  of  LO,  TO  and  interface  polaritons  (IP) 
in  polar  material.  Hybridization  arises  as  a  consequence  of  the  necessity  to 
satisfy  simultaneously  electromagnetic  and  mechanical  boundary 
conditions.  Optical  mode  analogues  of  surface  Rayleigh  waves  appear  in 
non-polar  slabs,  and  in  polar  superlattices  the  involvement  of  interface 
polaritons  causes  significant  dispersion  and  anisotropy.  Results  from  a 
continuum  theory  are  shown  to  describe  experimental  results  reasonably 
well.  Hybridization  affects  the  quantum  properties  due  to  the  coherent 
mixing  of  LO,  TO  and  IP,  and  this,  in  turn,  affects  the  interaction  with 
electrons.  Estimates  of  the  electron-hybridon  scattering  rate  in  GaAs 
quantum  wells  are  presented. 

3.2  Introduction 

The  advent  of  layered  semiconductor  structures  has  generated  the  need 
for  a  simple  analytical  description  of  the  way  in  which  confinement  of 
electrons  and  phonons  affects  their  interaction.  While  effective-mass 
theory  has  been  found  to  be  adequate  for  electron  confinement,  obtaining 
an  equivalent  theory  for  optical  phonons  has  been  more  problematic, 
particularly  in  polar  material.  The  main  difficulties  have  been  to  do  with 
boundary  conditions  and  the  nature  of  the  basic  interaction  of  confined 
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modes  with  electrons.  Recently,  continuum  models  of  confined  optical 
modes  which  overcome  these  difficulties  to  varying  degrees  have  begun  to 
emerge  [1-3]  ,  and  to  show  agreement  with  Raman  scattering  data  [4].  A 
subject  of  recent  discussion  [5,6]  has  been  the  nature  of  the  electron- 
polariton  interaction  -  whether  it  is  via  a  vector  or  a  scalar  potential  -  and 
this  has  now  been  resolved  [7].  An  outline  of  a  continuum  model  for 
confined  optical  vibrations  is  presented  here,  along  with  some  of  the  results 
generated  by  it. 

Section  3.3  reviews  the  basic  equations  which  describe  bulk  waves  in  an 
elastically  isotropic  continuum.  Section  3.4  discusses  the  boundary 
conditions  which  forces  the  mixing  of  optical  modes  and  Section  3.5 
describes  the  form  of  the  hybrid  modes  produced.  The  interaction  of  these 
hybrid  modes  with  electrons  is  discussed  in  Section  3.6.  Some  results  of 
the  theory  are  presented  in  Section  3.7  and  the  dispersion  predicted  in  an 
AlAs/GaAs  superlattice  is  compared  with  experiment.  It  is  also  shown  that 
hybrid  modes  can  have  either  LO-like  or  IP  (interface  polariton)-Iike 
properties,  and  the  strength  of  interaction  with  electrons  can  be  highly 
anisotropic.  The  paper  concludes  with  a  brief  discussion  in  Section  3.8  of 
some  of  the  problems. 

3.3  Continuum  Optical  Modes 

For  simplicity,  the  medium  supporting  optical  vibrations  is  considered 
to  be  an  isotropic  elastic  continuum  in  which  long-wavelength  modes  obey 
an  equation  of  motion  of  the  form; 


(1) 


where  u  is  the  relative  displacement  of  the  ions,  p  is  the  reduced  density, 
(Oto  is  the  natural  frequency  of  vibration,  ei  is  the  effective  ionic  charge 
density  (zero  for  non-polar  material),  E  is  the  macroscopic  electric  field, 
and  cii,  C44  are  optical-mode  elastic  constants  which  quantify  the  stress  in 
the  medium  and,  associated  with  it,  the  dispersion  relation.  In  the  absence 
of  free  charge  the  electric  field  must  satisfy 

V£(£0)E  =  0  (2) 


and  any  transverse  part  of  the  field  will  obey  the  electromagnetic  wave 
equation: 


_ T_  _  ^ 

(9/2  pi  £(co) 
0 


(3) 


where  e((o)  is  the  permittivity  (m  general,  frequency  dependent)  and  po  is 
the  permeability  of  free  space  in  the  case,  as  assumed  here,  of  non¬ 
magnetic  material.  The  ionic  displacement  can  be  factorized  into 
longitudinal  and  transverse  parts: 

u  =  u^  +  u^,  Vxu^=0,  V  Uy.=0,  (4) 


and  travelling-wave  solutions  for  longitudinal  optical  (LO),  transverse 
optical  (TO)  and  polaritons  (P)  emerge,  with  the  properties: 


LO  Qr  =  (olQ-wlk",  £(o))  =  0,  ^l  =  -;^^l  =  -[p(^Io-^to]  “l 
(a) 


TO  0)^  =  (0^-v^k^,  £(o))  =  oo,  E^  =  0, 


(b)(5) 


lf2  (qj^  —CO^) 

P  (0^= - .  £(«y)  =  f_ — ^ E^  =  sE 


“(6)2  -o)^)  ^ 

'to  ^ 


(C) 


where  ci>lo““to’  ^,and£_  are  the  zero-  and  high-frequency 

permittivities, 

S  =  (6)2 -6)2  )/(6)2  —6)2  ). 

These  three  types  of  mode  provide  the  basis  of  a  continuum  model  of 
confinement. 

3.4  Boundary  Conditions 

At  each  interface  both  electromagnetic  and  mechanical  connection  rules 
must  be  satisfied.  (Fundamentally,  there  are  only  electromagnetic 
conditions,  but  it  is  convenient  to  distinguish  short  range  elements  by  the 
label  "mechanical".)  The  EM  conditions  are  the  standard  ones,  i.e,, 
continuity  of  tangential  field  and  of  perpendicular  electric  displacement. 
The  mechanical  conditions  are  less  obvious.  For  acoustic  waves  these 
conditions  would  amount  to  entailing  the  continuity  of  displacement  and  of 
stress  for  equilibrium  to  be  maintained.  Macroscopic  stress  must  enter  the 
conditions  for  optical  modes,  but  since  the  latter  involve  internal  motion 
within  a  primitive  unit  cell,  the  situation  is  more  complicated,  as  has  been 
demonstrated  using  a  linear  chain  model  [8].  However,  in  linear  chain 
models  the  interface  is  assumed  to  be  located  with  precision,  but  its  precise 
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location  affects  the  transfer  matrix,  and  so  in  reality  some  averaging  is 
necessary,  and  this  reduces  the  importance  of  microscopic  stresses.  Thus 
taking  a  boundary  condition  based  upon  macroscopic  stress  only  is  likely  to 
be  reasonably  realistic.  The  mechanical  connection  rule  is  therefore  taken 
to  be  the  continuity  of  macroscopic  stress  and  displacement,  in  an 
"acoustic"  approximation. 

The  situation  here  is  not  unlike  that  for  electrons.  The  wavefunction 
for  electrons  is  a  Bloch  function  consisting  of  a  microscopic  part 
describing  the  charge  in  each  unit  cell  and  an  envelope  function.  Strictly, 
both  parts  must  be  considered  in  ensuring  continuity  of  current  at  a 
boundary,  but  it  is  often  a  good  approximation  to  restrict  the  matching 
procedure  to  the  envelope  function  only.  Our  "acoustic  approximation"  is 
a  somewhat  analogous  procedure  which  throws  all  the  emphasis  on  the 
macroscopic  property. 

3.5.  Hybrid  Optical  Modes 

In  order  to  satisfy  these  boundary  conditions  for  p-polarized  modes 
(i.e.,  modes  with  polarization  in  the  plane  of  incidence)  it  is  necessary  to 
take  a  linear  combination  of  the  bulk  solutions  having  identical  frequencies 
and  identical  in-plane  wavevector.  (For  the  s-polarized  TO  mode  this 
procedure  is  unnecessary.  It  is  only  with  p-polarized  modes  that  problems 
exist.)  In  a  single  layer  -a/2  <  z  <  a/2  with  the  x-axis  parallel  to  the 
direction  of  propagation,  the  resultant  p-polarized  hybrids  have  relative 
ionic  displacements  of  the  form: 

r  1 

u  =k  A  cosk.  z+c^cosk^z-c  cosk  z  e  * 

X  X  I  LT  Tp  pj 


(6) 
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Uz  =  ikj^A| 


k 

sink,  2-c^ — 2L_sink^z  +  c  — 5_sink  z 
L  ^k,k^  T  pk  k  P 

L  T  L  p 


i(k  x-fiX) 
le  * 


where  the  subscripts  L,  T  and  P  refer  respectively  to  the  LO,  TO  and  P 
components.  A  second  form  with  sine  replacing  cosine  also  exists.  In  the 
surrounding  material  I  z  I  >  a/2,  the  solutions  in  most  cases  consist  of 
components  which  exponentially  decay  away  from  the  layer.  For  given 
frequency  and  in-plane  wave  vector  the  wave  vector  kt,  kj  and  kp  are  all 
determined.  The  boundary  conditions  select  the  allowed  frequencies  and 
quantify  the  fractional  amplitudes  cj  and  cp.  Standard  quantization 
procedure  determines  the  amplitude  A,  and  expresses  u  in  terms  of 
annihilation  and  creation  operators.  (Strictly,  hybrids  of  the  same  parity 
should  be  subjected  to  a  rigorous  orthogonalization  procedure  but  in  many 
cases  of  interest  the  departure  from  orthogonality  is  small  and  can  be 
neglected.  Of  course,  opposite  parity  modes  are  always  orthogonal).  In 
non-polar  material,  Cp  =  0.  In  polar  material  at  in-plane  wavevectors  of 
order  106  cm-1,  typically  involved  in  the  interaction  with  electrons,  the 
unretarded  limit  can  be  taken  for  the  polariton,  which  then  turns  out  to  be 
an  interface  mode  with  kp  =  -  ikx.  At  frequencies  above  ©to  in  polar 
material  the  TO  component  is  also  an  interface  mode  with  kj  =  -  ioj.  For 
frequencies  significantly  above  (Oto,  oct  is  typically  so  large  that  its  precise 
magnitude  cannot  be  predicted  accurately  by  a  long-wavelength  model  with 
quadratic  dispersion,  and  it  is  then  necessary  to  obtain  it  from  the  complex 
vibrational  band  structure  predicted  by  lattice  dynamics. 

3.6  Interaction  with  Electrons 

Discussion  here  will  be  limited  to  the  polar  interaction.  The  latter  is 
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Double  modulation-doped  AIGaAs/  InGaAs  heterostructure  with  a  graded 
composition  in  the  quantum  well 


affected  by  hybridization  in  several  ways.  First,  the  Frohlich  interaction 
with  LO  modes  is  reduced  by  confmement  (which  puts  a  lower  limit  on  the 
total  wavevector)  a  reduction  which  increases  with  decreasing  well-width. 
Second,  the  interaction  with  IP  (interface  polariton)  modes  is  also  modified 
by  confmement  in  some  cases  through  the  coherent  association  with  the  LO 
component.  When  this  occurs  this  interaction  also  diminishes  with 
decreasing  well-width.  However,  for  strongly  polar  materials  in  which  the 
LO  and  TO  bands  do  not  overlap  (e.g.,  AlAs)  there  exist  frequencies  for 
which  the  LO  component  must  be  replaced  by  a  mode  drawn  from  the 
LO/LA  branch  of  the  complex  vibrational  band  structure.  In  this  case  the 
IP  mode  would  act  approximately  as  an  unhybridized  mode,  and  its 
interaction  with  electrons  would  then  increase  with  diminishing  well-width. 
(Such  a  behaviour  is  predicted  in  the  case  of  AlAs/GaAs  from  lattice 
dynamics  [[9]).  It  turns  out  that  in  the  case  of  a  slab  with  stress-free 
surfaces,  hybridization  is  weak,  and  the  IP  mode  is  expected  to  interact 
with  electrons  approximately  in  the  familiar  unhybridized  way.  But  for 
cases  where  the  mechanical  mismatch  effectively  reduces  the  relative  ionic 
displacement  to  zero  at  the  interface,  and  where  LO/IP  hybridization  can 
take  place  strongly,  the  polar  interaction  with  well  modes  will  under  some 
circumstances  diminish  in  strength  with  diminishing  well-width.  Finally, 
the  polar  interaction  will  be  affected  by  the  more  complex  dispersion 
which  occurs  as  a  result  of  LO/IP  hybridization.  This  is  particularly  the 
case  in  superlattices,  where  the  dispersion  becomes  highly  anisotropic,  and 
consequently,  so  does  the  interaction  with  electrons. 

The  interaction  of  interface  polaritons  with  electrons  takes  place  (in  the 
usual  Coulomb  gauge)  via  a  vector  potential  A  -  specifically  as  -e  A  • 
p/m*,  where  p  is  the  momentum  operator  and  m*  is  the  effective  mass  of 


the  electron.  Polaritons,  being  transverse  EM  modes,  do  not  have  a  scalar 
potential  (t),  and  so  the  IP  interaction  is  not  like  that  for  LO  modes. 
However,  in  many  cases  it  has  been  assumed  without  proof  that  in  the 
unretarded  limit  it  is  permissable  to  calculate  the  interaction  as  if  the  IP 
mode  had  a  scalar  potential  associated  with  it.  This  matter  has  been 
clarified  recently  [7].  In  the  unretarded  limit  for  interface  polaritions  it  is 
possible  to  associate  a  potential  A  with  the  mode  such  that  A  =  -^VA/e. 
A  unitary  transformation  generated  by  the  operator  U  =  exp  (-LA)  then 
converts  the  Hamiltonian,  to  first  order,  into  one  which  contains  an 
operator  which  behaves  like  a  scalar  potential  ())  plus  a  self-energy  term. 
For  real  transitions  involving  matrix  elements  on  the  energy  shell  the  use 
of  ep  gives  exactly  the  same  results  as  e  A  .  p/m*,  but  for  virtual  and 
higher  order  transitions  this  would  not  be  true.  Thus,  either  A  or  p  can  be 
used  in  the  unretarded  limit  provided  only  scattering  rates  are  required. 

The  interaction  with  holes  involves  the  deformation  potential  associated 
with  all  three  components  of  the  hybrid,  and  this  will  be  true  also  for 
electrons  in  L-valleys.  While  deformation  potential  scattering  by  LO  and 
TO  modes  is  familiar  enough,  that  by  interface  polaritons  is  novel. 

3.7.  Results 

The  theory  has  been  applied  to  the  cases  of  a  slab  with  stress-free 
surfaces,  and  to  a  quantum  well  and  a  superlattice  with  rigid  interfaces.  In 
the  case  of  the  non-polar  slab  the  theory  mixes  LO  and  TO  guided  waves 
and  predicts  the  existence  of  a  surface  wave  which  is  the  optical  mode 
analogue  of  the  Rayleigh  wave  (101.  In  the  case  of  a  polar  slab, 
hybridization  turns  out  to  be  weak  [11]  and  the  modes  can  be  described  to  a 
good  approximation  purely  by  taking  EM  boundary  conditions,  as  was 
done  by  Fuchs  and  Kliewer  long  ago  [12j. 


In  the  case  of  a  polar  quantum  well  with  u  =  0  at  the  interfaces, 
hybridization  is  strong  at  all  frequencies  provided,  there  is  no  LO/TO  gap 
at  the  zone  boundary[3j.  Such  is  the  case  with  GaAs,  and  Fig.  3.1  shows 
the  dispersion,  with  the  influence  of  the  IP  mode  evident  in  producing  anti¬ 
crossings  and  mode  conversions.  In  particular,  the  LOl  mode  with  kL  = 

:t/a  converts  to  an  L03-like  mode  with  Icl  =  2.86  Ji/a,  while  the  L03  moves 
toward  LOS.  Fig.  3.2  shows  the  intrasubband  scattering  rates  at  threshold 
in  an  infinitely  deep  quantum  well  associated  with  LO-like  hybrids. 

Although  the  IP  component  of  the  hybrid  enhances  the  interaction 
significnatly  form  the  bare  LO  contribution  it  gets  weaker  with  narrowing 
wells. 

Approximately  pure  IP  modes  exist  at  and  around  special  points  in  the 
disperion.  These  points  are  where  the  IP  dispersion  crosses  the  dispersion 
of  an  LO  mode  of  symmetry  opposite  to  that  of  the  IP,  and  are  thus 
identified  by  special  pairings  of  kt,  kx  .  These  IP-like  hybrids  will  interact 
with  electrons  somewhat  differently  from  the  interaction  associated  with 
LO-like  hybrids.  Being  limited  to  small  regions  of  kt,  kx  space  they  will 
not  be  as  generally  available  as  LO-like  modes  for  scattering  purposes  and 
so  their  contribution  to  the  overall  scattering  rate  is  expected  to  depend 
more  sensitively  on  the  specific  structure.  This  adds  to  the  difficulty  of 
estimating  scattering  rate  in  the  presence  of  hybridization. 

The  superlattice  dispersion  of  GaAs  hybrids  when  tanh  Oja  =  1  is  given 
by 

sink,  a(2rs^d  +  (k^ /k‘ )sinh  k  asinhk  b)  +  2s(k  /k,  ) 

L  X  L  X  X  X  L 

sinhk  b  +  rsinhk  acosk  (a +  b)-cos  k,  a(cosh  k  asinhk  b  +  rsinhk  acoshk  b)l  =  0 

X  xz  L  XX  XX  J 
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Figure  .3.2.  Intrasubband  scattering  rates  at  threshold  for  LO-Iike  modes  in  a  GaAs  infinitely  deep 
quantum  well. 
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Figure  3,4.  The  change  of  frequency  of  the  LOl  mode  from 
on-axis  to  in-plane  propagation  as  a  function  of  m  the  number 
of  monolayers  in  a  m/m  GaAs/AlAs  superlattice.  Open  squares; 
theory:  closed  squares  experiment  (4}. 
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Figure  3.5.  Patterns  of  potential  for  the  first  three  symmetric 
modes  in  a  GaAs  quantum  well  of  a  model  GaAs/AlAs  superlattice. 
The  barrier  width  is  5.66A  and  the  in-plane  wavevector  equals  that 
at  the  threshold  of  emission  (2.52  x  lO^cm-l).  (a)  Q  s  kz(a+bO  =  0; 
(b)  Q  =  T. 
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Figure  3.6,  Intrasubband  scattering  rates  at  threshold  for  GaAs  modes  in 
the  model  Gaas/AlAs  superlattice  with  b  =  5.66A;  contributions  of  the  first 
three  symmetrical  modes  are  shown,  (a)  O  =  0;  (b)  O  =  t. 
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Fieure  3.7.  Intrasubband  scattering  rates  at  threshold  for  well  and  barrier 
modes  for  the  same  system  as  Figs.  5  and  6;  (a)  Q  =  0;  (b)  Q  =  x.  Note 
that  the  sum  ignores  the  difference  in  threshold  energy  of  well  and  barrier 
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In  spite  of  the  above  agreement  with  experiment,  it  is  too  early  to  be 
entirely  confident  that  a  continuum  approach  to  the  problem  of  optical 
vibrations  in  layered  material  is  generally  adequate.  Indeed,  we  have 
already  seen  that  in  some  circumstances  modes  are  called  for  whose 
amplitudes  vary  significantly  over  a  distance  of  order  of  the  unit  cell 
dimension.  Such  components  of  a  hybrid  cannot  be  described  adequately 
without  help  from  lattice  dynamics.  Experimental  data  is  still  sparse  for 
layered  material  in  which  the  elastic  boundary  condition  is  not  as  simple  as 
u  =  0,  such  as,  for  example,  AlxGai-xAs/GaAs  with  x  <~  0.3,  though  this 
system  introduces  the  added  complication  of  a  two-mode  spectrum.  For 
the  polar  slab  with  stress-free  surfaces  hybridization  turns  out  to  be  too 
weak  for  the  approach  to  be  distinguishable  from  a  conventional  dielectric 
continuum  model.  Which  is  support  of  a  kind,  since  there  is  a  lot  of 
experimental  evidence  for  the  latter  model  in  this  particular  case.  A 
prediction  for  the  non-polar  slab  is  the  existence  of  optical  Rayleigh  waves, 
and  this  should  be  amenable  to  experimental  test. 

The  importance  of  a  continuum  model  resides  in  its  usefulness  for 
transport  theory.  Its  predictions  concerning  the  scattering  of  electrons  in 
real  systems  have  still  to  be  worked  out.  Unfortunately  the  paradigmatic 
system,  ALAs/GaAs,  contains  the  problem  of  describing  dispersion  in  ALAs, 
which  is  small.  How  small  determines  the  degree  of  LO/IP  hybridization 
and  this,  in  turn,  crucially  determines  the  strength  of  fringing  fields  in  the 
GaAs  which  can  scatter  electrons.  Hitherto,  it  has  been  assumed  that  no 
LO/IP  hybridization  occurs  in  AlAs  [15],  but  this  assumption  needs 
examination.  It  leads  to  the  conclusion  that  scattering  in  narrow  wells  is 
dominated  by  the  unhybridized  AlAs  IP  mode,  and  there  is  experimental 
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support  for  this  [16].  Nevertheless,  the  whole  question  of  hybridization  in 
materials  with  zone-edge  gaps  in  the  optical  mode  spectrum  needs  to  be 
investigated. 

The  results  of  the  scattering  rate  in  a  superlattice  with  narrow  opaque 
barriers  presented  here  shows  that  IP -like  hybrids  belonging  to  the  well 
make  a  significant  contribution.  Although  hybridization  produces  highly 
modified  properties  it  is  still  possible  to  identify  LO-like  and  IP-like 
hybrids,  which  are  responsible  for  highly  anisotropic  scattering  rates. 

Finally,  it  should  be  recalled  that  i'>  deformation  potential  interaction 
with  holes  or  with  L-valley  elec  ins  will  involve  the  interface  polariton 
component  of  a  hybrid  as  well  as  the  LO  and  TO  parts.  A  comprehensive 
description  of  the  interaction  of  hybrid  optical  modes  with  holes  does  not 
yet  exist,  but  is  certainly  needed. 

3.9  References 

[1]  X.  Ziarmi,  P.N.  Butcher  and  I.  Dharssi  J.  Phys.  Condens.  Matt.  4  L77 
(1992). 

[2]  KJ.  Nash  Phys.  Rev.  14^7723  (1992). 

[3]  B.K.  Ridley  Proc.  SPIE  Symp.  on  'Compound  Semiconductor  Physics 
and  Devices  -  High  Speed  Electronics  and  Optoelectronics’,  Somerset, 
NJ  1675  492  (March  1992). 

Phys.  Rev.  ML  4592  (1993). 

[4]  M.  Haines  and  G.  Scamarcio  NATO  ARW  "Phonons  in 
Nanostructures"  St.  Filiu,  Spain,  (1992). 

[5]  M.  Babiker  and  B.K.  Ridley  B43  9096  (1991);  B15  (1993)  (reply  to 
comment  by  Knipp  and  Reinecke). 

[6]  P.A.  Knipp  and  T.L.  Reinecke  Phys.  Rev.  B 15  (1993)  (comment  on 
above  paper). 

[7]  M.  Babiker,  N.C.  Constantinou  and  B.K.  Ridley,  submitted  for 
publication 

[8]  H.  Akera  and  T.  Ando  Phys.  Rev.  B40  2914  (1989). 

[9]  H.  Rucker,  E.  Molinari  and  P.  Lugli  Phys.  Rev.  B45  6747  (1992). 

32 


» 


» 


» 


» 


» 


» 


I 


» 


I 


» 


» 


•  •  • 


High-current  lattice-strained  lno.59Gao.4iAs/lno,52Alo.48As  modulation-doped 
field-effect  transistors  grown  by  molecular  beam  epitaxy 

I  Q  * 


» 


[10]  B.K.  Ridley  Phys.  Rev.  B44  9002  (1991). 

[11]  B.K.  Ridley,  O.  Al-Dossary,  N.C.  Constantinou  and  M.  Babiker 
submitted  for  publication 

[12]  R.  Fuchs  and  K.L.  Kliewer  Phys.  Rev.  AMO  2076  (1965). 

[13]  N.C.  Constantinou,  O.  Al-Dossary  and  B.K.  Ridley  Solid  St.  Commun 
M  191  (1993). 

[14]  T.  Tsuchiya,  H.  Akera  and  T.  Ando  Phys.  Rev.  B22  6025  (1980). 

[15]  O.  Al-Dossary,  M.  Babiker  and  N.C.  Constantinou  Semicond.  Sci. 
Technol.  2  B91  (1992). 

[16]  K.T.  Tsen,  K.R.  Wald,  T.  Ruf,  P.Y.  Yu  and  H.  Morkoc  Phys  Rev, 
Lett.  £7  2557  (1991). 

[17]  K.  Huang  and  B-F  Zhu  Phys.  Rev.  B38  13,  377  (1988). 

[18]  B.K.  Ridley  Phys.  Rev.  B29  5282  (1989). 


I 


» 


» 


33 


» 


•  •  • 


IV. 

Double  modulation-doped  AlGaAs/ InGaAs  heterostructure  with  a  graded 
composition  in  the  quantum  well 
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A  new  double  modulation-aoped  AlOaAs/InGaAs  heterostructure  with  a  graded  composition 
in  the  InGaAs  quantum  well  (QW)  has  been  designed  and  electncally  characterized.  .An 
InGaAs  QW  with  a  rectangular-like  potential  profile  in  the  presence  of  the  two  dimensional 
electron  gas  is  obtained  by  a  two-step  grading  of  the  In  composition,  which  results  in  a  broad 
and  symmetric  electron  distribution  profile  even  under  various  voltages.  The  Hall  measurement 
shows  a  very  high  electron  mobility  of  7230  em'/V  s  and  an  electron  sheet  density  of  4.1  x  lO'"/ 
cm"  at  room  temperature.  To  our  knowledge,  this  is  the  highest  mobility  ever  reported  so  far  for 
the  double  modulation-doped  .AIq  jGa^ -As/Ing -Gao  jAs  field-effect  transistor  structure. 


Pseudomorphic  AlGaAs/InGa.As/GaAs  heterostruc¬ 
tures  have  been  extensively  investigated  in  recent  years  to 
improve  the  performance  of  modulation-doped  field  effect 
transistors  (MODFETs)  as  well  as  opto-electronic  de¬ 
vices.'  The  InGaAs,  when  grown  within  a  certain  critical 
layer  thickness  (CLT)  on  Ga.As.  can  provide  a  larger  con¬ 
duction  band  discontinuity  at  the  .AlGaAs/InGa.As  inter¬ 
face  than  that  of  AlGaAs/GaAs.  which  leads  to  an  in¬ 
crease  m  the  electron  sheet  density  in  the  quantum  well 
(QW;.  The  electron  sheet  density  also  increases  with  the 
In  amount  in  the  InGaAs  QW.  though  the  thickness  of  the 
well  IS  limited  by  a  CLT  capable  of  accommodating  a 
misfit-dislocation-free  interface."  The  increase  of  the  elec¬ 
tron  sheet  density  (.Vj  in  the  QW'  is  crucial  to  the  perfor¬ 
mance  of  MODFETs  including  the  gain  cutoff  frequency 
for  high-speed  operation,'  noise  figure  in  high  frequency, 
low-noise  operation.'*  and  current-dnving  capability  for 
high-power  operation.'  In  order  to  increase  .V,  without 
increasing  the  In  concentration,  several  MODFET  struc¬ 
tures,  such  as  a  double  modulation-doped  field  effect  tran¬ 
sistor.”'  a  doped  quantum  well  FET.'  and  a  double  quan¬ 
tum  well  MODFET."*  were  proposed.  These  approaches, 
however,  often  cause  the  degradation  of  the  electron  trans¬ 
port  pr  ernes.  The  mobilities  of  the  electron  in  the  doped' 
and  doi  quantum  well  structure"  were  quite  low  be¬ 
cause  o.  ne  ionized  impurity  scattenng  in  the  QW  in  the 
former  structure,  and  the  interference  of  each  QW  with  the 
donor  supplying  layer  between  the  QW's  in  the  latter  case 
In  double  modulation-doped  structures,  it  is  very  difficult 
to  obtain  high  mobility  charactenstics  due  to  the  surface 
roughness  at  the  interface  below  the  QW'  channel'®  and  the 
Si  segregation  effect."  Si  atoms  in  the  donor  supplying 
layer  below  the  QW'  move  nonuniformly  to  the  vicinity  of, 
or  even  into,  the  QW  during  the  molecular  beam  epitaxy 
(  MBE ,  growth.  Therefore,  one  of  the  most  important  con¬ 
cerns  in  MODFETs  is  to  increase  the  .V,  without  sacrific¬ 
ing  the  electron  mobility  for  practical  high-frequency  ap¬ 
plications  In  this  letter,  we  describe  a  new  design  of  double 
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modulation-doped  heterostructure  and  characterize  the 
electrical  properties  of  the  MBE-grown  wafers. 

A  typical  double  modulation-doped  field  effect  transis¬ 
tor  (D-MODFET)  structure  is  shown  in  Fig.  1.  It  consists 
of  an  Ino-Gaoj.As  QW,  the  In  composition  of  which  is 
uniform  over  the  QW'.  two  GaAs/Alg  iGao  -As  spacers  ( .A 
and  B).  and  two  .Alo -.Ga^ -  As  barriers  including  Si  atomic 
planar  doping  (.APD)  layers.  In  panicular,  we  put  a  20  .A. 
undoped  Ga.As  layer  between  the  .Al,,  ;Gao  -.As  spacer  and 
Ino  -Gay  .-.As  QW.  which  greatly  helps  reduce  the  surface 
roughness  at  heterointerfaces,  especially  with  high 
.Al  (jOGt)  and  In  (2092)  concentrations,  and  also  allows 
for  the  change  in  the  growth  temperature  from  580  °C  for 
AlGaAs  to  500  °C  for  InGa.As.  "  The  necessity  of  the  20  A 
of  the  Ga,As  layer  between  the  pseudomorphic  InGa.As 
channel  and  AlGa.As  spacer  suggests  that  there  may  be 
serious  interface  roughness  at  the  heterointerface,  even 
though  the  theoretical  room-temperature  mobility  is  dom¬ 
inantly  determined  by  the  optical  phonon  scattenng  rather 
than  the  remote  ionized-impunty.  or  interface  roughness, 
scattenng. ' ' 

To  better  understand  why  the  interface  problems  fre¬ 
quently  take  place  in  the  MODFET  structure,  we  calcu¬ 
lated  the  electron  distribution  in  the  channel  by  solving  the 
one-dimensional  Poisson's  equation  and  Schrodinger  equa¬ 
tion  consistently  and  quantum-mechanically  along  the 
depth  direction.''*  The  simulation  produces  conduction 
band  energy,  electron  distribution  profile,  and  subband  en¬ 
ergy  levels  at  different  gate  bias  voltages.  Figure  2(a)  pre¬ 
sents  the  conduction  band  and  electron  profile  of  the  con¬ 
ventional  D-MODFET  corresponding  to  the  structure  of 
Fig  1  at  zero  gate  voltage.  Two  asymmetnc  peaks  of  the 
distribution  are  located  very  close  to  each  interface,  which 
indicates  that  the  electron  transpon  can  be  very  sensitive 
to  interface  properties  The  two  peaks  will  become  closer  to 
each  interface  and  w  ill  be  separated  from  each  other  fur¬ 
ther  if  the  QW  thickness  becomes  larger  Moreover,  assum¬ 
ing  that  the  electron  profile  can  be  regarded  as  the  super¬ 
position  of  two  distributions,  each  with  its  own  peak  (one 
closer  to  the  top  interface  and  the  other  to  the  bottom 
one),  they  will  probably  interfere  together  along  the  irans- 
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ized  with  *  network  analyzer  and  a  pair  of  microwave 
probes  calibrated  by  the  short-open-load-through  tech¬ 
nique  on  an  impedance  standard  substrate.  The  current 
gain  cutoff  frequency /,  extrapolated  at  —  6  dB/octave 
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QW  .  Ino  iAi  Lniform-Composiuon  Scracture 
.  In^Gai-^Ai  .  GracJcO  Structure 

FIG  1.  Schemauc  cross  section  of'  conventional  and  graded  double 
modulation-doped  structures  ;  leu  QW  convenuonai  siruciure.  nght  QW 
graded  structure; 

port  path  where  the  aspect  ratio  of  the  typical  source-drain 
distance  (1  am)  to  the  QW  thickness  (100  A)  is  more 
than  100.  The  interference  of  two  maior  peaks  seems  to  be 
another  reason  that  the  mobility  of  D-MODFETs  reponed 
so  far  has  not  been  as  good  as  that  of  the  single-doped 
MODFET''  (7  000  em'-W  s).  Most  of  the  best  reponed 
mobility  for  typical  Air,  ;.Gao -.As/ In/,  ^Ga^  jAs  D- 
MODFETs  are  scattered  around  6  000  cm‘/V' s  with  .V,  of 
3.5— i.Ov  10‘'/cm'  at  room  temperature. 

To  improve  the  interface-related  mobility  characteris¬ 
tics  of  the  D-MODFET.  we  designed  a  new  Alo  jGaorAs/ 
In,);Ga<i  jAs  D-MODFET  structure  by  linearly  grading 
the  in  composition  in  the  channel.  W'e  anempted  to  make 
the  conduction  band  in  the  channel  flatter  than  that  of  the 
conventional  D-MODFET  by  decreasing  the  In  composi¬ 
tion  at  the  top  and  the  bottom  interfaces,  and  increasing 
the  In  composition  at  the  center  of  the  channel.  This  com¬ 
positional  grading  modifies  the  potential  shape  of  the  QW, 
resulting  in  a  rectangular-like  QW,  which  pushes  both  elec¬ 
tron  peaks  near  the  interfaces  toward  the  QW  center.  As  a 
result,  a  broad  and  relatively  symmetric  electron  distribu¬ 
tion  can  be  obtained,  as  shown  in  Fig.  2(b).  If  we  try  to 
cause  a  more  sophisticated  In  profiling  such  as  a  nonlinear 
or  a  multistep  grading,  more  symmetric  electron  distribu¬ 
tions  farther  from  the  interface  can  be  obtained  under  any 
operating  gate  voltage  and  donor  density.  It  is  worthwhile 
to  point  out  that  the  energy  difference  between  the  first  and 
the  second  subbands  of  the  graded  structure  [Fig.  2(b)] 
case  is  larger  by  1 1  meV  than  that  of  the  conventional 
D-MODFET  [Fig.  2(a)].  which  implies  that  there  is  less 
probability  of  the  iniersubband  scattering. 

Based  on  the  above  band  gap  grading  concept,  we  grew 
a  conventional  D-MODFET  structure  (sample  No.  4129) 
with  uniform  In  concentration  of  20^  in  the  QW.  and  a 
compositionally  graded  one  (sample  No.  4130)  by  MBE. 
In  the  graded  structure,  the  In  concentration  is  15%  at 
bottom  interface,  linearly  increased  up  to  22%  through  the 
center,  and  decreased  down  to  15%  at  the  top  interface. 
The  average  In  concentration  is  18.5%.  .All  other  struc¬ 
tural  and  growth  paramelers  c.scepi  in  the  QW  siruciurc 


FIG  2.  Conduction  band  and  electron  distribution  profiles  at  f',=0  for 
the  D-MODFET  (a)  conventional  structure  with  uniform  20%  In  con¬ 
centration,  (b)  graded  structure  with  linear  gradient  of  lS%-22%-t5% 
In  concentration,  (c)  graded  structure  with  linear  gradient  of  5%-25%- 
5%  In  conccauation. 

were  identical  in  both  structures.  In  both  cases,  we  used 
the  well-calibrated  Si  APD  technique  and  the  90  A 
Alg  jGao.7As/lO  A  GaAs  superlattice  buffer  layers. 

For  accurate  Hall  measurements,  we  chose  3  samples 
from  different  positions  in  a  2  in.  wafer  for  each  structure 
and  made  the  van  der  Pauw  clover  leaf  patterns.  We  re¬ 
peated  Hall  measurements  after  each  chemical  etching  of 
the  GaAs  cap  layer  in  a  small  step  (stripping  Hall  mea¬ 
surement).  The  mobility  and  the  sheet  density  of  the  con¬ 
ventional  D-MODFET  were  measured  to  be  6520  cm‘/V  s 
and  4. lx  10'’/cm’  at  room  temperature,  respectively.  On 
the  other  hand,  those  of  the  graded  structure  were  the 
e.xtremely  high  mobility  of  7230  em’/'V  s  and  the  sheet 
density  of  4.1  x  lO'Vcm’  at  room  temperature.  The  77  K 
mobility  and  the  sheet  den.sity  are  also  as  high  as  21  630 
em'/V  s  and  3  9  ■:  lO'Vcm’,  respectively  To  our  knowl- 
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The  In-plane  effective  mass  In  strained-layer  quantum  wells 
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A.  -  K 

4129 

20 

20 

20 
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22 

15 
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4  I 

3  9 

4195 

23 

20 

20 

20/  yo 
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7  260 

2  2 

20 

4199 

5 

2^ 

< 

0/30 

5. '90 

16  050 

2  5 
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edge,  this  high  mobility  with  such  a  high  sheet  density  of 
4.1  X  10'‘/cm"  is  the  best  ever  reponed  for  the  .\lGa.A.s/ 
InGaAs  D-MODFET.  This  mobility  over  7(300  cm'/V  s  at 
room  temperature  is  quite  sig.iificant  because  with  such  a 
high  electron  density,  some  parts  of  parallel  conduction 
may  still  take  place  in  the  Alp^Gao-As  barrier  layers, 
where  the  electron  mobility  is  quite  poor,''*  Table  I  sum¬ 
marizes  the  QW  structures  and  Hall  measurement  data  at 
3(30  and  77  K. 

To  verify  the  effectiveness  of  the  compositional  grading 
on  the  mobility,  we  tned  another  MBE  growth  run.  We 
intentionally  degraded  the  mobility  by  using  higher  growth 
temperature  to  enhance  the  Si  segregation."  and  using 
thinner  .AIq  ^Ga^  -.As  spacer  ( 30  .A  instead  of  50  .A ) .  We 
grew  two  kinds  of  D-MODFET  structure:  one  (sample 
No.  4198)  is  the  conventional  structure  with  a  uniform  In 
composition  of  209^.  and  the  other  (sample  No.  4199)  is 
the  graded  one  with  a  steeper  gradient  in  In  concentration. 
The  In  composition  was  linearly  changed  from  5%  at  both 
interfaces  to  25%  at  the  center,  as  summanzed  in  Table  I. 
In  the  latter  case,  we  did  not  use  the  20  .A  GaAs  spacer 
because  the  In  concentration  at  the  heterointerface  is  only 
5%.  According  to  the  simulation  shown  in  Fig.  2(c).  the 
electron  density  is  strongly  concentrated  around  the  QW' 
center  with  the  full  width  at  half  maximum  less  than  60  A. 
Although  the  mobility  was  degraded  on  both  structures  as 
expected,  the  graded  structure  showed  much  less  degraded 
mobility  of  5390  cm'/V  s.  compared  to  that  of  the  conven¬ 
tional  structure.  3030  cm'/V  s.  This  seems  to  suggest  that 
the  grading  of  the  composition  is  indeed  efficient  in  con¬ 
trolling  the  electron  profile  and  providing  the  relative 
mobility-insensitivity  to  the  interface  conditions. 

In  conclusion,  the  potential  shape  of  the  InGaAs  quan¬ 
tum  well  in  AlGaAs/InGaAs  MODFET  structures  can  be 
controlled  by  the  band  gap  grading  of  the  well.  The  elec¬ 
tron  distnbution  profile  in  the  rectangular-like  InGaAs 
QW  is  broad  and  symmetric.  Hall  measurement  shows  that 


the  very  high  mobility  and  sheet  density  can  be  obtained  at 
the  same  time  in  the  graded  structure,  which  can  improve 
the  performance  of  .MODFETs, 
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fonduction  band  and  valence  bands  can  be  decoupled,  and 
when  the  spin-orbii  splitting  is  large  the  split-offhand  can  be 
decoupled  froni  the  rest.  As  a  result  the  8x8  Hamiltonian 
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Progress  in  High  Frequency  Heterojunction  Field  Effect  Transistors 


Lester  F.  Eastman 

School  of  Electrical  Engineering  and  National  Nanofabrication 
Facility,  Cornell  University.Phillips  Hall, 

Ithaca,  New  York,  14853-5401 

Abstract.  The  concepts,  technology  and  the  present  experimental  limits  of 
performance  of  heterojuncdon  modulation  doped  field  effect  transistors  are 
covered.  Both  lattice-matched  and  strained,  pseudomorphic  quantum  well 
channels  on  GaAs  substrates,  as  well  as  lattice-matched  quantum  well 
channels  on  InP  substrates  are  included.  Power  gain  frequency  limits  to  450 
GHz,  and  noise  figures  as  low  as  .3-.9  db  at  60  GHz  are  presented  for  300K 
operation  of  the  latter  devices. 


Microwave  and  millimeter  wave  field  effect  transistors  have  reached 
impressive  performance  levels  by  the  use  of  heterojunctions.  In  these 
structures,  a  lower  band-gap  layer  holds  the  electrons,  while  a  higher  band- 
gap  layer  is  doped  with  donor  atoms.  As  well  as  lattice-matched 
combinations  on  GaAs  and  InP  substrates,  thin  strained  layers  have  been  used 
to  hold  the  electrons.  Very  short  gates  have  been  realized  using  electron- 
beam  lithography.  This  presentation  covers  the  materials  used,  the 
fabrication  methods,  the  materials  and  device  assessment  methods,  and  the 
present  performance  limits. 

2.  Materials 

The  initial  results  used  A1  jGa  7As/GaAs,  where  an  undoped  .5  -  1.0  ^m 
GaAs  buffer  layer  was  grown  on  the  GaAs  semi-insulating  substrate.  The 
electron  sheet  density  in  the  Ga.As  is  limited  to  .3  -  .9  x  10* -/cm-  in  such  a 
structure,  due  to  the  limited  conduction-band  potential  step  of  .24V.  At  low 
temperatures  there  is  a  high  electron  mobility,  ranging  from  90,000  cm^/V-s 
up  to  220,000  cm^/V-s  at  77K.  The  latter  result  is  only  possible  with  low 
electron  sheet  density  and  with  a  substantial  (-  200  A  )  undoped  region  in  the 
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Ai.GaAs  near  the  heterojunction.  Increasing  the  fraction  of  aimrinuni  does 
not  allow  higher  electron  sheet  density  because  the  energy-  of  the  donors  in 
the  Al.GaAs  does  not  rise  more  than  .16  eV  above  the  GaAs  conduction  band 
edge. 

In  order  to  increase  the  donor  energy  above  the  conduction  band  edge. 
Indium  can  be  added  to  a  thin,  pseudomorphic  channel  layer  holding  the 
electrons.  With  In  ^5Ga  in  this  channel,  electron  sheet  density  values  to 

1.6  -  1.7  X  10^^/cm-  are  achievable,  and  with  In  75AS  the  sheet  density 

can  reach  Z.4  -  2.5  x  10'“/cm-.  In  this  latter  case,  the  conduction  band 
potential  step  is  .44V.  Mobility  of  electrons  is  lowered  a  modest  amount  at 
300K,  but  is  severely  lowered  at  77K,  dropping  from  90,0QJ^cm-/V-s  for  no  In, 
to  1S,000  cm~/V-s  for  25%  Indium. 

Using  InP  substrates,  Al  ^gln  52As/Ga  47ln  33AS  structures  have  been  used  to 
achieve  improved  results.  In  this  case  the  conduction  band  potential  step  at 
the  heterojunction  ;s  .52  eV  and  the  electron  sheet  density  can  be  at  least  3  x 
’0^-/om-.  In  addition,  electrons  have  only  two  thirds  as  much  effective 
mass,  which  raises  the  room  temperature  mobility  by  about  50%,  The  77K 
electron  mobility  is  limited  by  alloy  scattering,  and  to  date  has  not  surpassed 
the  60,000-70,000  cm-/V-s  region. 

In  order  to  achieve  high  performance  in  transistors  with  short  gates,  electrons 
must  be  confined  by  a  potential  barrier  from  beneath  as  vvell  as  on  top. 
Instead  of  thick  layers  of  the  barrier  allov,  a  superlattice  such  as 
A1  3Ga  pAs/GaAs/Al  jGa  jAs  etc,  is  used. 


On  top  of  the  doped  barrier  layer,  which  usually  contains  Aluminum,  there  is 
normally  placed  a  cap  layer  having  no  Aluminum,  This  layer  is  doped,  as 
needed,  to  limi.  the  depletion  of  the  desired  electron  sheet  density  along  the 
channel.  The  ohmic  contact  alloying  can  more  easily  penetrate  such  cap 
layers  having  no  Aluminum  due  to  the  absence  of  Aluminum  oxide. 

3.  Transistor  Structures 

The  modulation-doped  field  effect  transistors  (MODFET's)  are  fabricated 
using  electron-beam  lithography.  High  performance  is  achievable  with  short 
gates,  with  -  .15  4m  gate  lengths.  Such  gates  are  made  with  a  cross  section 
having  a  mushroom  shape,  or  T-shape,  with  a  short  "footprint"  but  with  an 
enlarged  metal  conductor  region  above  the  semiconductor.  Such  structures 
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can  be  fabricated  with  S  200  I^/mm  resistance  along  their  length.  Without  the 
mushroom  shape,  such  gates  would  have  8-10  times  as  high  a  resistance.  It 
would  require  six  gate  fingers,  each  with  one  sixth  the  length,  to  achieve  the 
same  performance  as  two  fingers  with  mushroom  shape. 

The  layout  of  the  MODFET  pattern  should  be  designed  to  have  low  gate-to- 
drain  feedback  capacitance,  in  order  to  yield  a  high  value  of  power  gain.  This 
allows  the  unity-power-gain  frequency,  ffnax'  ^-5 ' 2.0  times  as  high  as 
the  unity-current-gain  frequency  fj. 


On  GaAs  substrates,  the  pseudomorphic  InyGa^.yAs  channel  has  interrelated 

limits  on  thickness  and  Indium  fraction.  Above  y  s  ,22  the  surface  of  the 
In.GaAs  is  often  rough,  lowering  yield.  At  y  =  .25  the  maximum  thickness  of 
this  layer  is  limited  to  -  100  A,  since  dislocations  form  for  thicker  layers.  As 
shown  in  Figure  1,  this  causes 
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For  quanrum  well  thickness 
=  0,  this  is  simply  a  GaAs 
channel  MODFET.  The  gf^ 
and  f  j  are  both  about  50% 
higher  for  the  optimum. 

In  25Ga  -j.^ks  quantum  well, 
compared  with  a  GaAs 
quantum  well.  Figure  2  shows 
the  current  gain  versus 
f.'equency  for  a  .MODFET  with 
an  optimum  In  25Ga  -^.As 
channel  and  .15  qm  x  150  pm 
mushroom  gate.  The  room- 
temperature  ft  shown  is  the 
present  record  value  for 
pseudomorphic  InGaAs/GaAs. 
(Nguyen,  et  al,  1989)  The  room- 
temperature  of  a  MODFET 
with  the  same  quantum  well 
and  .15  x  6-0  nm  gate  was 
250  GHz.  With  lower  f^  (100  GHz), 

fmax  been  as  high  as  350  GHz  on 
such  devices  (Lester  et  al.  1988). 
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Fig.  1.  Mutual  transconductance 
of  In  25Ga  75  As /GaAs  MODFETs 
versus  strained-laver  thickness. 
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Fjg.  2  ^Current  gain  performance  of  In.25Ga ^gAs/GaAs  pseudomorphic 
0.15xl50pm  MODFET  versus  frequency. 


Figure  3.  Conduction  band  potential  profile  of  an  A1.43ln.52As/Ga.48ln,52As/ 
Ga.47ln  53As/InP  MODFET  with  superlattice  buffer  layer. 


Progress  in  hi^li  frequencv  heiernniiurion  held  eli'eci  transisiors  62.’ 

Devices  ofjh^  type  have  \-ielded  1.5  db  r.sise  dgure  at  60  GHz  for  .15  jim  gates 
at  room  temperature. 

MODFETs  with  A1  43ln  52AS  barriers,  and  Ga  ^yln  53AS  channels  lattice- 
matched  to  InP  have  yielded  the  highest  frequency  performance  for  a  give 
gate  length.  Figure  3  shows  a  scale  drawing  of  the  conduction  band  profile  for 
such  a  device.  The  conduction  band  potential  step  is  .52  V  at  the  interface, 
allowing  electron  sheet  density  up  to  -  3  .x  10^-/cm-  as  well  as  -  50%  higher 
mobility  values  than  for  quantum  well  channels  on  GaAs.  In  addition, 
electron  average  transit  velocity  is  at  least  one  third  higher.  There  are  still 
problems  related  to  gate  leakage,  channel  breakdown,  electron  traps  in  the 
AlInAs,  and  reliability  with  these  devices  on  InP,  but  dedicated  research  to 
improve  these  areas  has  just  begun.  These  devices  have  yielded  .S-.9  db  noise 
figures  at  60  GHz  for  .15  )im  gates  at  room  temperature.  The  best,  room 
temperature  values  for  f^  (Mishra)  and  (Chao  and  Smith  (a))  are  250  and 
450  GHz,  respectively  for  these  devices.  In  comparing  devices  on  Ga.As  and 
InP  substrates,  the  a.xpected  frequency  response  as  a  function  of  gate  length  is 
shown  in  Figure  4.  The  quantitative  values  for  the  ma.ximum  electron  sheet 
density,  and  average  electron  transit  velocity  are  given  in  Table  1. 


Fig.  4.  Expected  frequency  response  of  optimum  MODFETs  on  GaAs  and  InP 
substrates  versus  gate  length. 
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TABLE  1 

Channel 

Electron  Sheet  Density 

Average  Electron 

Transit  Velocity 

GaAs 

.8-.9x10''2/cm2 

1.2x10^  cm/s 

'".25^3.75^® 

2.4-2.5x10''2/cm2 

1.8x10^  cm/s 

2®.47''^.53As 

3-3.2x10''2/cm2 

2.4x10^  cm/s 

A.  Conclusions 

It  is  concluded  that  MODFET's  have  the  highest  frequency  resoonse 
and  the  lowest  noise  figures  of  any  transistors.  Those  MODFET's  on 
InP  substrates  have  the  highest  f(  and  f-nax  lowest  noise 

figure,  while  those  MODFET’s  on  GaAs  substrates  have  the  highest 
power  (Chao  and  Smith  (b))  and  efficiency  (Chao  and  Smith  (b)) 
(.66W/mm  at  41%  efficiency  at  60  GHz  and  .36  W/mm  at  23% 
efficiency  at  94  GHz). 
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High-current  driving  capability  is  realized  in  submicron-gate  lattice-strained  Ing  59Gao4|As/ 

(Aa/a  =  4x  10~^)  modulation-doped  field-efTect  transistors.  Full-channel 
drain  current  in  excess  of  1.10  and  1.90  A/mm  are  obtained  at  80  K  for  the  singly  doped 
and  doubly  doped  structure,  respectively.  By  using  the  double  modulation  technique 
and  a  buried  p-buffer  layer,  excellent  high-frequency  performance  and  flat  transconductance 
characteristics  were  realized  over  a  very  broad  range  of  gate  and  drain  bias  voltages. 

These  devices  are  very  suitable  for  large-signal  or  power  device  operations.  The  high  quality 
Ing.j2M).48'A^  buffer  layer  eliminates  the  hysteresis  and  current  instability  (or  the  kink 
effect)  in  submicron-gate  devices  at  both  300  and  80  K. 


Excellent  high-frequency  performance  was  previously 
reported  in  the  lattice-matched  ItiojsGao  47AS/ 
In^  57AI0.41AS  modulation-doped  field  effect  transistors 
(MODFETs).'"  However,  most  of  these  devices  utilize  a 
single  electron  supplying  layer,  and  the  current  driving 
capability  is  limited.  Lattice-strained  lOxGai.^^As/ 
InojjMx'W'A^  0.33)  modulation-doped  beterostruc- 
tures  promise  higher  electron  mobility  and  larger  conduc¬ 
tion  band-edge  discontinuity  than  the  lattice-matched 
structures.^’  Higher  electron  mobilities  and  velocities  are 
expected  in  the  strained  structures  due  to  the  better  elec¬ 
tron  confinement  and  reduced  scattering.*’’  Larger 
conduction-band-edge  discontinuity  in  the  strained  struc¬ 
tures  also  results  in  a  better  electron  confinement  and  a 
larger  sheet  electron  concentration.  Reduced  scattering*  in 
the  channel  further  raises  the  electron  mobility  and  effec¬ 
tive  velocity.  In  this  letter  we  report  the  realization  of  high- 
frequency,  high-current  MODFETs  by  utilizing  a  lattice- 
strained  quantum  well  chaimel  and  doubly  doped  wide-gap 
electron  supplying  layers.  A  full  channel  drain  current  over 
1  A/mm  is  demonstrated  in  both  the  strained  singly  doped 
(SD)  and  doubly  doped  (DD)  InGaAs/InALAs 
MODFET  structures  with  6%  excess  indium  in  the  chan¬ 
nel.  The  large  current  density  together  with  the  higher 
thermal  conductivity  of  InP  than  GaAs  makes  these  de¬ 
vices  potential  candidates  for  high-speed  large-signal  appli¬ 
cations  such  as  millimeter-wave  power  amplification  and 
laser  drivers  for  multigigabit  optical  link. 

The  strained  SD  and  DD  heterostructures  were  grown 
by  molecular  beam  epitaxy  on  the  semi-insulating  InP  sub¬ 
strates.  Schematic  cross-sectional  layer  structures  and  the 
conduction-band-edge  diagrams  are  shown  in  Figs.  1(a), 
1(b),  and  1(c),  respectively.  For  the  SD  structure,  both 
the  capping  and  electron  supplying  layer  have  a  silicon 
doping  level  of  3  X  10'*/cm*.  For  the  DD  structure,  the 
silicon  doping  level  of  the  capping,  upper,  and  lower  elec¬ 
tron  supplying  layer  are  3x10'*,  3xl0'*,  and  3xl0'*/ 


cm*,  respectively;  and  the  beryllium  doping  level  in  the 
buried  p-buffer  layer  is  SxlO'Vcm*.  The  Inoj,Gao4,As 
channel  growth  was  performed  with  20%  lower  Ga  flux  to 
achieve  the  6%  higher  indium  mole  fraction.  The  potential 
barrier  of  the  p-buffer  layer  in  the  DD  structure  further 
improves  the  earner  confinement  and  isolates  the  two- 
dimensional  electron  gas  (2  DEG)  from  the  substrate  de¬ 
fects  and  impurities, 

MODFETs  were  fabricated  with  a  recess-gate  process. 
A  planar  isolation  was  realized  by  boron  ion  implanution 
at  33  keV  with  a  dose  of  1 X  10‘*/cm*.  Breakdown  voltages 
exceeding  200  V  were  measured  between  two  active  regions 
separated  by  a  15  pm  spacing,  justifying  the  effectiveness  of 
the  isolation  process.  From  the  four-point  probe  measure¬ 
ments  of  the  annealed  Ni/AuGe/Ag/Au  ohmic  contacts, 
contact  resistances  of  0.14  and  0.13  11  mm,  specific  contact 
resistivities  of  1.79x10“*  and  1.89x10"*  fl  mm%  and 
sheet  resistivities  of  1 10  and  92  fl/0  were  measured  at  300 
K  for  the  SD  and  DD  structure,  respectively.  The  sheet 
resistivity  was  measured  with  the  capping  layers.  At  80  K. 
the  sheet  resistivity  dropped  below  50  fl/O  for  both  struc¬ 
tures  as  the  result  of  the  increase  in  the  carrier  mobilities 
and  2  DEG  densities.  A  recess  etch  was  performed  to  re¬ 
move  the  n  *  -InGaAs  capping  layer  before  depositing  the 
electron  beam  defined  T-shaped'®  Ti/Pd/Au  gates  on  the 
200  A  undoped  InAJAs  layer.  This  is  essential  to  reduce 
the  gate  diode  leakage  current.  The  measured  dc  end-to- 
end  gate  resistances  were  smaller  than  200  fl/mm  for  gate 
footprints  of  various  submicron  dimensions. 

Figure  2  shows  a  typical  drain  current-voltage  charac¬ 
teristics  of  a  0.36x30  pm*  gate  device  for  the  DD  struc¬ 
ture  at  300  K.  Similar  characteristics  are  obtained  in  the 
SD  structure  with  a  lower  current  level.  At  a  bias  point  of 
K,.  =  2.3V  and  Fp  s  0  V,  the  devices  are  able  to  deliver  a 
current  density  of  0.9  (3(X)  K)  and  1.07  (80  K.)  A/mm  for 
the  SD  structure  and  1.3  (  300  K)  and  1.8  (80  K)  A/mm 
for  the  DD  structure.  There  is  a  20%  increase  of  drain 
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FIG.  2.  Cuirent'voluge  ciutacteracics  of  a  0.36x50  /im*  gate  doubly 
doped  device  at  300  K. 
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current  for  both  structures  at  77  K.  The  full  channel  drain 
current  is  over  1.10  A/mm  for  the  SD  structure  and  1.90 
A/mm  for  the  DD  structure  at  80  K.  The  high  current 
level  also  indicates  the  effective  activation  of  silicon  doping 
through  the  high  quality  epitaxial  growth  and  efficient 
charge  transfer  from  the  electron  supplying  layer  to  the 
quantum  well  channel  Hysteresis  and  drain  current  anom- 
aJities  (or  the  kink  effects)  in  the  current-voltage  curves,  as 
commonly  observed  in  the  short-gate  InP-based  devices, 
are  not  present  at  both  300  and  80  K.  We  attribute  this  to 
the  low  impurity  and  defect  densities  in  the  high  quality 
InAlAs  layer.  With  an  appropriate  heat  sink  for  power 
consumption,  the  devices  can  sustain  a  drain  voltage  up  to 
3.5  V. 

The  peak  e.xtrinsic  dc  transconductance  is  511 
(300  K)  and  645  (  80  K)  mS/mm  for  the  SD  structure  and 
545  (  300  K)  and  646  (80  K)  mS/mm  for  the  DD  struc¬ 
ture  at  =  1.0  V  as  shown  b  Fig.  3.  The  current  density 
under  the  high  g„  bias  condition  is  still  over  0.50  and  0.65 
A/mm  at  both  300  and  80  K  for  the  SD  and  DD  structure, 
respectively.  The  improved  transconductance  and 
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FIG.  I.  Epitaxial  layer  sinictura  of  the  lactice-ttnined  (a)  singly  and 
(b)  doubly  doped  In^,,CaB.4|At/In«.5]Ala.4|As  MODFETl.  (e) 
conduciiofl-band-cdge  diagranu  for  both  structures. 
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FIG.  3.  de  transconduetance  vs  gate  voluge  tor  a  036x  30  ftm  gate 
singly  and  doubly  doped  device  at  300  and  80  K.  (Solid  fines)  sing  y 
doped  device;  (broken  lines),  doubly  doped  device. 
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fig.  4.  Typical  revene-biaacd  current-volUge  eharaeteraties  of  a 
0.36x30  lun*  gaie-to-dnin  Schonky  diode  at  3C0  K. 

increased  drain  current  at  80  K  are  the  results  of  better 
electron  confinement  and  mobility  enhancement  from  the 
lattice-strained  heterostructures.  The  high  g„  is  maintained 
over  a  broad  range  of  from  0.8  V  and  up.  At  =  1.0 
gm  >6  greater  than  S20  mS/mm  in  both  structures  for  a 
full  range  of  from  —  1. .55  to  —  2.60  V  at  80  K.  This 
bias  insensitivity  of  g„  is  very  desirable  to  minimize  the 
harmonic  distortions  for  large-signal  operations. 

Figure  4  shows  the  reverse-biased  gate-to-drain 
Schottky  diode  characteristics  of  a  0.36x50  nm~  gate  de¬ 
vice  at  300  K.  A  breakdown  of  —  6  V  is  measured.  Al¬ 
though  not  quite  comparable  with  the  similar  GaAs/ 
AlGaAs  devices,  the  gates  perform  reasonably  well  among 
the  usually  leaky  InGaAs/TnAlAs  devices  especially  under 
the  present  condition  of  high  carrier  densities.  This  indi¬ 
cates  the  effectiveness  of  the  200  A  undoped  InAlAs  bar¬ 
rier  layer.  The  pinch-off  voltages  for  both  structures  are 
deep.  This  makes  high  drain  voltage  and  very  negative  gate 
voltage  bias  conditions  impractical  due  to  the  gate  diode 
leakage  current  consideration.  (  For  example,  at  =  3.2 
V  and  F’p  =  —  3  V,  a  6.2  V  reverse  voltage  is  applied  to 
the  gate  ^ode.)  To  operate  the  devices  in  the  medium  and 
ultimately  high  power  regime,  further  improvements  in  the 
gate  characteristics  can  be  achieved  by  employing  a  thicker 
InAlAs  barrier  layer  or  increasing  A1  mole  fraction  of  the 
InAlAs  barrier  layer. 

The  microwave  performance  of  devices  was  character¬ 


ized  with  a  network  analyzer  and  a  pair  of  microwave 
probes  calibrated  by  the  short-open-load-through  tech¬ 
nique  on  an  impedance  standard  substrate.  The  current 
gain  cutoff  frequency  /„  extrapolated  at  —  6  dB/octave 
from  Ayi*  is  Sl:fc2  and  42:t;l  GHz  for  the  SD  and  DO 
structure,  respectively.  The  maximum  stable  gain  is  1 1  ±  1 
and  1 1.64: 1  dB  at  26  GHz  for  the  SD  and  DD  structure, 
respectively.  The  lower/,  ( -I-  C^J)  in  the  DD 
structure  is  due  to  larger  parasitic  capacitances 
(Cp  +  Cpj)  despite  the  higher  g^  Both  the  current  and 
power  gains  vary  very  little  over  a  broad  gate  and  drain 
bias  range.  Qymbined  with  the  simultaneous  high  current 
density  and  improved  gate  breakdown  voltage,  this  struc¬ 
ture  is  very  suitable  for  the  high-speed  digital  and 
millimeter-wave  applications. 

In  summary,  we  have  fabricated  the  lattice-strained 
lno.]9Ga()^4,As/InoL52Alo.4«As  MODFETs  with  a  high  cur¬ 
rent  driving  capability  and  kink-free  characteristics.  The 
increase  of  conduction-band-edge  discontinuity  by  employ¬ 
ing  a  strained  structure  is  effective  in  enhancing  2  DEG 
densities.  The  broad  bias  range  of  gate  and  drain  voltages 
which  achieve  high  transconductance  and  high  current  si¬ 
multaneously  is  very  attractive  for  large  signal  operations. 
With  further  optimization  in  gate  characteristics,  they  are 
potential  candidates  for  power  device  applications. 

This  work  was  supported  in  pan  by  GE,  ARO,  and 
AT&T. 
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The  in-plane  effective  mass  in  strained-layer  quantum  wells 

B.  K.  Ridley^' 

Sch(H>l  of  Electrical  Engineering.  Cornell  University.  Phillips  Hall.  Ithaca.  .\ew  York  14853-5401 
( Received  28  March  1990;  accepted  for  publication  3  July  1990) 

The  problem  of  calculating  the  valence-band  structure  of  strained-layer  quantum  wells  in  the 
effective-mass  approximation  is  reviewed.  Using  the  spherical  approximation  and  exploiting 
the  simplicity  of  the  infinitely  deep  well  model  we  show  that  the  in-plane  effective  mass  is 
determined  by  two  factors — a  splitting  contribution  which  is  dominant  at  large  strains,  and  a 
quantum  confinement  contribution.  A  model  for  finite-depth  wells  is  presented  which  gives 
analytic  expressions  for  the  zone-center  in-plane  mass  and  associated  nonparabolicity  factor, 
and  it  is  applied  to  the  system  In.Ga,  ,  As/GaAs.  The  model  allows  the  computation  of 
valtiKe-band  structure  using  no  more  than  a  pocket  calculator.  It  is  shown  to  give  results  in 
tsonable  agreement  with  experiment. 


I,  INTRODUCTION 

It  has  been  known  for  over  30  years  that  uniaxial  or 
biaxial  strains  can  split  the  heavy-  and  light-hole  valence 
bands  in  semiconductors  at  the  zone  center. '  An  interesting 
result  of  this  splitting  is  the  phenomenon  of  mass  reversal  in 
which  the  heavy-hole  band  acquires  a  light  mass  in,  say,  the 
XV  plane  but  retains  its  heavy  mass  in  thez  direction,  whereas 
.he  light-hole  band  acquires  a  heavy  mass  in  the  xy  plane  but 
retains  its  light  mass  in  the  z  direction,  the  xyz  directions 
being  determined  by  the  orientation  of  the  strains.  The  po¬ 
tential  of  this  structure  to  exhibit  negative  differential  resis¬ 
tance  (NDR)  was  realized  long  ago,’  and  recently  its  ad¬ 
vantages  for  lasers  regarding  the  provision  of  a  low  threshold 
current,'  reduced  intervalence-band  absorption,  and  re¬ 
duced  Auger  recombination"'  has  led  to  a  major  activity  in 
the  investigation  of  strained  layers.'  The  latter  are  obtained 
by  growing  an  epilayer  on  a  lattice-mismatched  substrate, 
and  these  layers  are  stable  and  dislocation-free  provided  the 
layer  thickness  does  not  exceed  a  critical  value,  which  is  typi¬ 
cally  50-1000  A  (the  bigger  the  strain,  the  smaller  the  criti¬ 
cal  thickness).  Interest  has  therefore  focused  on  strained- 
layer  quantum-well  structures  in  which  the  effects  of 
quantum  confinement  are  added  to  those  associated  with 
strain. 

Theoretical  work  in  this  area  is  based  on  the  effective- 
mass  formalism  of  Luttinger  and  Kohn,"  '  a  formalism 
w  hich  works  in  the  case  of  quantum  confinement  only  pro¬ 
vided  that  the  periodic  part  of  the  BIcKh  function  in  the  well 
is  the  same  as  that  in  the  barrier.  Where  this  is  not  the  case  a 
full-blown  band-structure  calculation  cannot  be  avoided. 
Even  when  the  effective-mass  formalism  holds  well,  keeping 
track  of  the  altered  mixing  of  heavy-  and  light-hole  states 
brought  about  by  confinement  warrants  a  numerical  treat¬ 
ment.  The  problem  has  been  discussed  by  Altarclli''  and 
Mar/in.'  When  the  growth  of  the  layer  is  along  a  (III) 
direction  the  complexity  is  compounded  by  the  appearance 
of  a  substantial  electric  field  via  the  piezoelectric  effect  '  In 
many  cases  in  practice  the  effective-mass  approximation  is 
expected  to  give  reasonable  results,  and  we  will  restrict  our 
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attention  to  these  cases  in  what  follows.  For  purposes  of  de¬ 
vice  modeling  it  is  desirable  to  minimize  the  amount  of  nu¬ 
merical  computation,  limiting  the  latter  if  possible  to  what 
can  be  done  on  a  pocket  calculator.  To  this  end  there  are 
several  simplifying  assumptions  which  can  be  adopted,  and 
these  will  be  briefly  reviewed  in  the  next  section.  Our  aim 
will  be  to  evolve  out  of  these  a  simple  model  which  contains 
the  essential  features  of  the  physics,  and  we  will  apply  this  to 
theln.Ga,  .As.^GaAs  system. 

11,  BASIC  THEORY 

The  basic  problem  in  the  context  of  the  effective-mass 
formalism  and  k*p  theory  is  to  obtain  the  allowed  electron 
energy  states  taking  into  account  the  coupling  between  two 
conduction-band  states  (spin  up-spin  down),  and  six  va¬ 
lence-band  states  (three spin  up-spin  down),  but  neglecting 
changes  brought  about  by  more  remote  bands.  (This  is  the 
first  approximation.)  The  effect  of  strain  via  deformation 
potentials  can  be  incorporated  into  the  Luttinger-Kohn 
Hamiltonian,  as  can  piezoelectric  fields,  but  we  will  limit 
attention  to  simple  biaxial  or  uniaxial  strains  along  (100) 
directions  so  that  piezoelectric  fields  do  not  exist.  Piezoelec¬ 
tricity  is  connected  with  the  lack  of  a  centre  of  inversion 
symmetry  in  the  crystal  lattice,  and  this  produces  an  asym¬ 
metry  in  the  energy  structure,  but  this  effect  is  small  and  we 
will  neglect  it.  Each  eigenfunction  is  a  linear  combination  of 
the  eight  basic  wave  functions.'’  which  are  themselves  linear 
combinations  of  P){  x).  r).  z> )  orbitals  chosen  to  proper¬ 
ly  quantize  total  angular  momentum  along  the  z direction  in 
the  presence  of  the  spin-orbit  interaction.  In  bulk  material 
the  Hamiltonian  is  an  8  ■  8  matrix  with  the  allowed  energies 
being  given  by  the  solutions  of  an  8  ■  8  secular  determinant. 
When  an  interface  is  present,  continuity  of  wave  function 
and  current  leads  toa  lb  ■  lb  determinant  to  solve  (32  ■  32 
if  no  inversion  symmetry  ) . 

While  this  presents  no  difficulty  in  principle  many  ofthe 
inherent  details  which  expand  the  need  for  computational 
power  are  often  of  minimal  interest  for  example,  the  aniso¬ 
tropic.  as  distinct  from  the  nonparabolic,  features  of  hand 
structures.  Moreover,  the  interaction  of  bands  vviiiely  sepa¬ 
rated  in  energy  is  weak  and  this  can  form  the  basis  of  a  signif¬ 
icant  simplification.  Thus  when  the  band  gap  is  large  the 
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conduction  band  and  valence  bands  can  be  decoupled,  and 
when  the  spin-orbit  splitting  is  large  the  split-offband  can  be 
decoupled  from  the  rest.  As  a  result  the  8  x  8  Hamiltonian 
matrix  splits  into  twfo  2x2  matrices  and  one  4x4  matrix. 
This  is  a  pow'erful  simplification  which  is  often  Justifiable  in 
practice.  The  2x2  matrices  describe  bands  (conduction 
band  and  split-offband)  which  are  spherical  and  parabolic 
to  a  good  approximation. 

The  solution  of  the  4x4  secular  determinant  in  bulk 
material  gives  the  energ>  of  the  light-hole/heavy-hole  bands 
as  a  function  of  k  and  its  components  A:,,  A,,  k,  along  the 
cube  axes  as  follows: 

E  =£•,  +  (fr/2m){-  r,k'  ±  [4}'3A:"-t-  12(r;  -  f:) 
\(k\k  I  +  k  Ik  I  +  k]k  I ) 

+  2r.{k\  -^-k:  -2k:)K-  +  K*]'  -},  (1) 

where  are  the  Luttinger  parameters,’  is  the 

shift  in  energy  due  to  the  hydrostatic  component  of  the 
stress,  and  K  is  the  equivalent  wavevector  associated  with 
the  splitting  strain  energy  5,  viz., 

K-  =  2mS/fi\  (2) 

In  the  case  of  a  biaxial  stress, 


^  T  (cii  -C12)  .  (c,,  -i-2c,, ) 

E,=2a - e,,  S  =  b - e,. 

Cm  Cm 


where  a  and  b  are  valence-band  deformation  potentials,  c, , 
and  C| ,  are  elastic  constants,  c,  is  the  strain  associated  with 
the  difference  between  lattice  constants  of  well  and  bar¬ 
rier  (the  latter  assumed  to  be  unstrained),  and  a,,  is  the  lat¬ 
tice  constant  in  the  well,  and  we  assume  that  compressive 
strains  are  positive.  In  Eq.  ( 1 )  the  upper  sign  refers  to  the 
heavy-hole  (HH)  band,  the  lower  to  the  light-hole  (LH) 
band. 

Each  of  the  four  eigenfunctions  is  of  the  form 

vk  =  (.4  5  M,!)  +  C -  1)  -h  £>  1^,  -  P  (4) 

They  can  be  most  easily  represented  as  column  matrices,’" 
thus 


H,{k)  = 


L,  (k)  = 


R, 

0 

-  E ' 

e*L  Z/j(k)  = 

-M- 

^  L 

0 

-L 

-M 

Ri 

0 

0 

L,(k)  =  ^ 

-A/- 

L 

where 

R,=P  Q+E-E,.  /?,=/*-(- C 

P^^r,k\  Q  =  +  kl  -2kl)  +6, 

2m  2m 
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I - —ik^)k,, 

2m 


M  =  -^  [r-v-MA  ;  A;)  -/y.2v3A.A  1. 

2m 

The  character  of  the  state  becomes  pure  heavy  hole  or 
pure  light  hole  when  A,  and  A,  are  zero.  In  this  case  Z-  and  M, 
the  off-diagonal  terms  in  the  Hamiltonian,  are  zero  and 


— (Ki  4  2r,  lA;  +5. 
Zm 


Thus,  R ,  is  finite  and  /? ,  is  zero  when  the  upper  sign  is  taken, 
which  corresponds  to  the  heavy-hole  state  with  effective 
mass  ratio  equal  to  (y,  —  2y, )  '.  The  lower  sign  gives  the 

pure  light-hole  state  with  effective-mass  ratio  equal  to 
(y,  +  2y, )  '.  The  two  states  are  split  by  an  energy  2b. 

Since  the  deformation  constants  a  and  b  are  negative,  com¬ 
pressive  strains  make  <5  negative,  which  means  that  the  HH 
band  lies  above  the  LH  band  (Fig.  1.) 

The  character  of  the  bands  becomes  increasingly  mixed 
with  increasing  A,, ,  [  the  wave-vector  component  in  the  (x,y) 
plane)  with  the  effective  mass  in  the  (jc,y)  plane  showing 
mass  reversal.  The  extent  of  this  mass  reversal  varies  with 
direction  in  the  plane,  the  magnitude  of  the  anisotropy  being 
determined  by  the  difference  between  the  Luttinger  param¬ 
eters  y,  and  y, .  In  many  cases  this  difference  is  small  and  it  is 
often  convenient  and  justified  to  take  a  spherical  average  and 
henceforth  ignore  such  anisotropy.  It  is  straightforward  to 
show  that 

<A;AJ-|-A:’A.^4-A;’A-:)  =  {A*.  (8) 


FIG  1,  Valence-hand  slruclurc  (schematic)  for  (a)  compressional  stress 
(<^''0).  ( h)  tenstoiial  stress  (/>  0),  Nfite  the  possibility  of  degeneracy  in 
(he  lalter  case. 
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Thus,  in  the  spherical  approximation,  which  we  will  adopt 
from  now  on, 

£=£,+—[  -y,k-  ±(4rk^ 

2m 

+  ly(k'  -lk,)K-  + 

(9) 

where 

r  =  (2^3  +  3}^  )/5  (10) 

and  we  put  y,  =-  =  y  elsewhere. 

Introducing  an  interface  means  that  we  have  to  combine 
forward  and  backward  traveling  waves  of  the  same  parity  on 
one  side,  and  fit  to  evanescent  functions  on  the  other  to  ob¬ 
tain  bound  states  conserving  amplitude  and  current  flow.  If 
the  potential  is  asymmetric,  such  as  in  the  case  of  a  single 
heterojunction,  we  must  consider  both  spins  separately, 
since  spin  degeneracy  is  not  maintained  when  the  system 
lacks  a  center  of  inversion  symmetry.  Here  we  will  limit  our 
discussion  to  the  symmetric  quantum-well  case,  so  we  need 
consider  only  one  spin. 

In  choosing  the  combination  of  forward  a  backward 
traveling  waves  in  the  quantum  well  it  is  not  merely  a  matter 
of  reversing  the  sign  of  k.,  because  more  than  one  quantiza¬ 
tion  of  angular  momentum  is  involved  in  general  as  a  conse¬ 
quence  of  the  mixed  light-hole/heavy-hole  character  of  the 
waves.  A  specular  reflection  is  equivalent  to  a  rotation  in  the 
(xy)  plane  followed  by  a  space  inversion.  Rotation  affects  a 
state  through  the  rotation  operator 

R,=e  (11) 

where  J,  is  the  component  of  the  total  angular  momentum 
along  the  z  direction,  which  in  our  case  takes  on  the  values 
j,  —  —  j’.  This  causes  the  sign  to  oscillate  according  to 

1,  -i-  1,  —  1,  -t-  1.  Thus  suitable  WJ’  ve  functions  for  the 
quantum  well  are 


//,  (k  k. )  = 


R,  cos  k.z 

—  iL  ■  sin  k^z 

—  M  ■  cos  k.z 

0 


i,(k,*,) 


—  L  cos  k.z 
iR;  sin  k.z 
0 

—  A/'  sin  k.z 


(12) 


These  two  wave  functions  are  degenerate  at  a  particular  en¬ 
ergy  E  for  wave  vectors  k.  =  k,,  (say)  for  the  HH  band  and 
k.  =  k,  for  the  LH  band,  so  that,  for  k  =  0, 

( +  2y)k  j  -  E'  =  (y,  ~  2y)k  +  E\  ( 13) 

The  envelope  function  inside  the  well  is  then  a  linear  combi¬ 
nation  of  the  form 

'*',.,11  =  ■'*^1  (k  k„)  -+■  BH^  (k,|k,  ) 

+  CL,{y,k„)  +DL,(k,k,).  (14) 

Note  that  in  the  presence  of  strain  k,,  or  k,  may  be  imagi¬ 
nary  even  though  they  are  in  the  well.  In  the  barrier  we  take  a 
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similar  linear  combination,  but,  for  bound  states,  and  k, 
are  both  imaginary  and  V„ ,  the  barrier  height,  is  added  to  the 
energy.  The  allowed  energies  are  obtained  by  solving  the 
8x8  determinant. 

For  k||  =  0  the  situation  simplifies  since  H,  (0,k,  )  =  0 
and  L,  {0,k„ )  =  0,  and  we  obtain  from  the  boundary  condi¬ 
tions  at  z  =  ±  L/2  the  familiar  conditions 

tan  k,L  /2-F 
or 

cot.,Z./2=  -F,  (15) 

F=  [»«■.(  -  EQ)/m'gEQ\''\ 

where  ky  =  k„  or  kj  ,ml.  and  m'g  are  the  z-direction  effec¬ 
tive  masses  in  the  well  and  barrier,  respectively,  and  E^  is  the 
quantum  confinement  energy,  which  is  the  kinetic  energy  of 
the  hole  measured  from  the  bottom  of  the  quantum  well.  But 
when  k,|  ^0  we  are  back  to  our  8x8  determinant. 


III.  SIMPLE  MODEL 

There  is  one  case  where  an  analytic  solution  can  be 
found,  namely,  when  the  well  is  infinitely  deep.  In  this  case 
the  wavefunction  vanishes  at  z  =  L  /2  and  there  is  no  current 
to  conserve.  As  a  result  the  8x8  determinant  reduces  to  a 
4x4  determinant.  Working  this  out  gives  the  following 
equation: 

[E-  +  E(.P„+P,,)+  P„P,  -  Q„Q,  -MM‘] 

Xsin  k„L  sin  k,  L  —  [2('f^/2m)^y^^k^^k„ki^ 

X  ( 1  —  cos cos  k,  L)  =  0,  (16) 

where  Pi,,Pl  and  Qh,Ql  are  given  in  Eq.  (16)  with 
ky  =  k„,ki^.  (In  the  spherical  approximation  we  replace 
by  —  y-)  When  k||  =  0  the  solution  is 

k„L  =  n7r  or  k,L  =  niT,  (17) 

where  n  is  a  nonzero  integer.  However,  we  note  that  k  „  ( and 
k,  )  deviates  from  the  familiar  condition  of  Eq.  (17)  when 
k,;  ytO.  The  quantized  wave  function  in  the  confinement  di¬ 
rection  is  thus  dependent  on  k,, .  Reference  to  Eq.  (9)  shows 
that  the  change  of  energy  with  k,  then  arises  from  terms 
directly  dependent  on  k^.  plus  a  contribution  from  the  vari¬ 
ation  of  the  confinement  energy  (i.e,  k^)  on  k„.  Both  of 
these  contributions  combine  to  determine  the  in-plane  effec¬ 
tive  mass. 

When  k||  is  small,  specifically  when 
k-<k]  (2rk:’-A’),  (18) 

Eq.  (16)  becomes,  in  the  spherical  approximation, 

2k„k,  (1  -  5„  )(1  ~  S,  )sin  k„L  sin  k,  L 

-H  3k  ’  ( I  —  cos  k„L  cos  k,  L)  ~0.  (19) 

where  we  have  introduced  the  dimensionless  quantities 
S„.S, .  viz. 

S„  K'/2rk'„.  S,=K'/2yk).  (20) 

These  strain  ratios  measure  the  strength  of  the  strain  relative 
to  the  strength  of  the  confinement.  Putting  k„  =  k,,,, 
-i-  Ak„,  Ak„  4k„„,  where  k„,,f,  =  hit.  we  can  determine 
the  variation  of  the  heavy-hole  confinement  wave  vector 


•  • 


•  • 
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with  ,k,  being  given  by  Eq.  (13),  We  can  also  perform  a 
similar  calculation  for  the  light-hole  band.  We  find  that  Ak„ 
is  simply  proportional  to  k  ’ .  Consequently,  if  we  expand  Eq. 
(9)  to  first  order  in  k  '  an  in-plane  effective  mass  for  the  HH 
band  can  be  defined  for  the  wth  subband  as  follows: 


(r.  -  '///)') 

3(;',  -  2y)  [cos  A,  +  (  -  1 ) 
(1  -  )  ( 1  ~  S,  )ki  L  sin  A,  L 


m-  (HH) 


(21) 


with  A,  (and  hence  S,  )  determined  by  Eq.  (13)  when 
A„  =  mr/L.  and 

r/„  -  (2 +5„)/(l  -S,,).  (22) 

Note  that  there  are  two  contributions  on  the  right-hand  side 
of  Eq.  (21 )  which  we  may  term  the  splitting  component  and 
the  confinement  component.  Similarly,  for  the  LH  band. 


D'l  +  1),  y) 

3(1';  4  2)')  [cos  A„/.  4-  (  —  1)" 


»;•  (LH) 


(1  ■S',/)(l  -  i>;  )k,,L  s\n  ki/L 


(23) 

with  k,i  (and  hence  S,,)  determined  by  Eq.  (13)  when 
A,  >177-/ L.  and 

V,  -  (2  +  Si  )/(  1  -  S,  ).  (24) 

When  the  strain  is  zero  these  expressions  reduce  to  those  of 
Nedorezov  '  ‘ 

Equations  (21  )-(24)  are  valid  provided  .S„  or  .V,  -  1 
When  Sii  -  1  it  turns  out  ( see  Eq.  (  13 )  |  that  A,,  =  A, .  and 
S,  =  1  also.  This  situation  can  happen  only  for  positive 
strain  energies,  corresponding  to  lattice  tension,  when  the 
l.H  band  is  above  the  HH  band  and  they  cross  in  the  confine¬ 
ment  (7)  direction  at  A„  =  k,  =  htt/L,  independent  of  A 
according  to  Eq.  (19).  Solving  Eq.  (9)  in  this  case  for  small 
A  gives 


E  =  {fr/2m)  [  -  (A  i,  -i-  A  ’ )  f  2A„A  )'\3]. 


(25) 


It  should  be  noted  that  special  cases  cKcur  for  all  LH  and  HH 
subband  degeneracies.  i.e.,  when  A„=«„!7/L, 
A,  =  7t/L  correspond  to  the  same  energy.  Whenever 
these  occur,  Eqs.  (21 )  and  (23)  are  no  longer  valid.  A  de¬ 
scription  in  these  cases  is  straightforward  to  obtain,  but  we 
w  ill  not  pursue  the  topic  of  subband  degeneracies  further. 

The  big  drawback  to  the  model  as  it  stands  is  its  depend¬ 
ence  on  the  assumption  of  an  infinitely  deep  well,  w  hich  is 
know  n  to  lead  to  a  very  poor  estimate  of  subband  energy  in 
real  cases  It  turns  out,  however,  that  the  in-plane  mass  is 
very  iiisensitive  to  well  width  in  the  infinitely  deep  well  miKl- 
el.  the  variations  m  the  splitting  and  confinement  compo¬ 
nents  lending  to  cancel  one  another  out  This  encourages  the 
assumption  that  the  mass  is  not  sensitive  to  the  details  of  the 
■onfinement  and  hence  we  can  use  the  in(iiiile-vvell  moelel  to 
give  us  an  estimate  of  the  /one-center  mass  In  order  to  do 
this  we  adopt  an  "equivalent  infinitely  ileep  well"  (LIW) 
model  In  this  we  obtain,  say  for  the  HH  confinement.  A„  by 
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solving  the  correct  boundary  conditions  with  A  =0  via  Eq. 
(15),  and  hence  define  S,,.  We  then  take  the  confinement 
contribution  to  the  in-plane  effective  mass  to  be  the  same  as 
in  the  case  of  an  infinitely  deep  well  of  width  L  '  such  that 
k„l.*  =  nrr.  This  means  that  we  can  use  Eq.  (21)  with 
L  —  L'  and  k,L'  determined  from  Eq.  (13)  with 
A„Z.  ■  =  riTT.  Effectively  we  are  lumping  all  the  effects  due  to 
the  portions  of  the  8x8  determinant  which  were  zero  in  the 
derivation  of  Eq.  ( 16)  into  the  effective  well  width  L".  This 
priK'edure  uses  the  correct  confinement  energy  and  hence 
preserves  the  correct  ratio  of  strain  to  confinement  strengths 
and  uses  the  infinite-well  model  to  estimate  the  confinement 
contribution  to  the  effective  mass. 

It  is  also  useful  to  determine  the  degree  of  nonparaboli- 
city  and  this,  too,  can  be  done  in  the  El  W  mixiel  ( though  the 
algebra  is  simiewhat  tedious).  We  define  a  nonparaboheity 
factor  a  as  follows: 

a'k  '/2m'  --  E  (1 


aE 


(26) 


where  E  is  the  in-plane  kinetic  energy.  To  obtain  o  we  ex¬ 
pand  A„  and  A,  as  follows: 


-  ok  ^  -h  {b/k'fi,^)k^ 


and 


k  i  —  k  1 ^  ck  ((J  /k  1  ..  )k  ^  (27) 

in  Eqs.  ( 9 )  and  (  16).  Here  o.h.c.  and  d  are  coefficients  to  be 
determined  (not  to  be  confused  with  deformation  potentials, 
etc.,  used  earlier!).  We  obtain  for  a  HH  band: 


a,  r- 


[4  S,!  4  40(1  -  5",,)] 


2(5(  1  -  S,,)' 


/m'\-(ri  -2y)'b 
\  m  )  E„„ 


where 


(28) 


(29) 

(.^0) 

(31) 

(32) 


3[cos  A;  ,,2,  ■  4  (  -  1 )"  '  ] 

(1  -  )  ( 1  —  5,  )A  ,  ,,2.  ■  sin  A,  ,,2  " 

,  ,  ,  [r,  -  2^1  -  2S„)|  ,  , 

A,  ,,2.  ■  .  - - - ITT-. 

E,,.,  -  (^'/2w)(;',  2y)k'„,,. 

ft  ^  (<■  -  /cot  A,  „L  '  )/g. 

In  the  expression  for  b.  we  have 
e  -  a  k  ,,2,  ■  ■  (  1  S/i )  (  1  S,  ) 


a[k  ,  ,,L  ‘  (  I  5",  )  (4o  1  *  ) 

r  tt  -  (  1  .V„)(4c  I  +  V;  I  ]  6t(  77-  c.  (33) 

/  2o(7i.tA,  ,,2.  ■(  I  S,i )  I  1  .S',  ). 

.g  4A,.,2.  "(1  .S„)(l  .V,  ) 

SmceA,  IS  chosen  to  giv  c  the  same  energy  asA„  the  factors 
and  c  are  relaleit.  viz  , 

<  i/„  ■  )/,  )  I  U) 

Lor  a  LH  band  we  interchange  .V,,  ami  .S',  ,  o  and  i ,  A  ,,  ,nul 
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TABLE 

-  t.S.Vv 

I.  Parameters 

+  0.45.r 

for  GaAs  and 

InAs. 

EjiOO  K)  =  1.43 

i 

a 

b 

r, , 

r,.  a. 

4 

(eV) 

(  UV  '  dyn  cm  )( A) 

GaAs 

7.0 

--  2.0 

-  17 

1  4.5 

11.88 

5.38  5  65.55 

InAs 

4-  4.7- 

-  1..V‘ 

-  1.7 

0..55 

8..5.5 

4.53  6.0584 

s5 

l'\ 

V: 

/  • 

r 

m'ff/m 

m,*  /m  m,‘  /m 

GaAs 

6.85 

2.10 

2.90 

2  61 

0.613’' 

0.0824"  0067 

InAs 

20.5 

8. 57 

9.29 

8.4,5 

0.379’' 

0.0261"  0.027 

‘  Derives  from  H,, 

,  a  -  t  eV  and  A/',  :^E 

-  S.,:  a< 

-  1:2. 

'  Derived  from  the  Luttinger  parameters  in  the  spherical  approximation 


k, and  change  the  sign  of  }'and  the  strain  energy  S  in  Eqs. 
(28)  -  (33). 

These  expressions  yield  the  correct  result  of  vanishing 
nonparabolicity  when  either  f  approaches  zero  or  the  strain 
splitting  energy  is  large  and  negative  (compression).  Taken 
with  the  expressions  for  zone-center  mass  [Eqs.  (21 )  and 
( 23 )  ]  they  provide  a  complete  analytical  description  of  the 
band  structure  near  the  zone  center,  which  is  applicable  to 
arbitrary  well  depths  and  widths  within  the  EIW  approxi¬ 
mation  and  is  exact  for  infinitely  deep  wells. 

It  is  time  to  summarize  the  discussion  of  this  section.  A 
study  in  the  case  of  an  infinitely  deep  well  shows  that  the  in¬ 
plane  effective  mass  near  the  zone  center  is  determined  by- 
two  components,  one  depending  on  the  splitting  of  the  HH 
and  LH  bands  brought  about  by  strain  but  modified  by  quan¬ 
tum  confinement,  the  other  depending  on  how  the  confine¬ 


ment  changes  with  in-plane  motion.  These  components 
quantify  the  degree  of  HH-LH  mixing  out  from  the  zone 
center.  In  the  case  of  a  well  of  finite  depth  the  confinement 
energies  can  be  determined  simply  only  at  the  zone  center. 
However,  provided  the  splitting  and  confinement  compo¬ 
nents  of  the  mass  are  close  to  those  for  an  infinitely  deep 
well — which  the  independence  of  mass  on  well  width  sug¬ 
gests  is  true — an  EIW  model  can  be  used  to  determine  the 
zone-center  mass  and  how  it  varies  away  from  the  zone  cen¬ 
ter. 


IV.  APPLICATION  TO  In^Ga,  As/GaAs  QUANTUM 
WELLS 

We  have  applied  these  ideas  to  the  system 
In,  Ga,  ,  As/GaAs.  The  parameters'-  used  in  the  calcula¬ 
tion  are  given  in  Table  I.  The  extrapolations  between  GaAs 
and  InAs  were  all  assumed  to  be  linear  except  for  the  band 
gap.  The  distribution  of  band-edge  discontinuities  arising 
from  alloying  was  taken  to  be  the  same  as  obtained  from  the 
GaAs  dilalational  deformation  constants,  viz.,  11%  for  the 
conduction  band.  Figure  2  shows  the  variation  of  subband 
energies  at  A  =0  and  overall  band  gap  (El  to  HHl )  with 
indium  compo.silion  up  to  40%  assuming  a  50-A  well.  Note 
that  the  LH  band  drops  below  the  GaAs  valence-band  edge 
implying  that  light  holes  are  unconfined.  The  system  is  type  I 
for  heavy  holes  but  type  II  for  light  holes.  (This  will  not  be 
the  case  if  the  interaction  with  the  split-off  band  is  strong. ) 

The  dependence  of  subband  structure  on  well  width  for 
X  =  0.3  is  shown  in  Fig.  3.  For  Z.  =  50  A  there  is  only  one 
electron  subband  and  two  HH  subbands.  Figure  4  depicts 
the  variation  of  in-plane  mass  with  composition  and  with 


FIG  2  F‘nerg>  levels  as  a  function 
of  indium  cimiposilion  (  x)  for  u 
5{)-A  In,Ga.  ,  As/GaAs  quan¬ 
tum  well  (a)  Suhband  levels  and 
band  edges,  fb)  energv  gap  G-l- 
HHl ) 
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FIG.  3.  Subband  structure  as  a  function  of  well 
width  for  x  =  0.3.  The  discontinuities  in  the 
conduction-band  (A£', )  and  valence-band 
( A£,  )  edges  are  also  shown. 


well-width  calculated  using  Eq.  (21)  with  the  equivalent 
infinitely  deep  well  approximation.  This  shows  that  the 
zone-center  mass  remains  remarkably  constant  with  varying 
well  width  for  the  odd-numbered  subbands.  Figure  5  shows 
the  dependence  of  the  non-parabolicity  factor  on  x  and  L. 

It  is  interesting  to  compare  these  results  with  the  experi¬ 


mental  results  obtained  from  a  study  of  magnetolumines¬ 
cence  by  Jones  et  al.'^  The  comparison  is  shown  in  Table  II. 
We  have  used  a  7:2  ratio  for  the  band  offsets.  The  energy 
gaps  differ  by  3.5%,  1.7%,  and  1.7%,  respectively,  for  the 
three  samples,  which  is  a  very  reasonable  agreement.  Apart 
from  the  case  of  the  sample  with  25%  In  the  masses  exhibit 


KK i  4  Varialii.n  of  /onf.cfnltT  cffccio i-  ntas-,  wilh  ( a  )  conipo.'itM.n  for  a  So. A  well  aiul  ( fr  t  well  wiillfi  for  \  O  f 
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FIG.  5.  Dependence  of  the  nonparabolicity  factor  a  on  (a)  composition  for  a  50-A  well  and  (b)  well  width  for  x  = 


0.3. 


TABLE  II.  Comparison  of  EIW  predictions  with  experiment. 


Specimen 

In 

% 

L 

(A) 

( meV ) 

ar(eV)  ' 

Expt. 

Theor. 

Expt. 

(  ±  8%) 

Theor. 

Expt. 

(  ±  20%) 

Theor. 

AK242 

15 

125 

1388.1 

1340 

0.084 

53.3 

22.4 

BC042 

20 

80 

1326  2 

1304 

0.086 

0.084 

29.8 

18  3 

TO  195 

25 

90 

1288.3 

1266 

0.094 

18.1 

15,7 

excellent  agreement.  Comparison  with  the  experimental 
nonparabolicity  factors  shows  only  rough  agreement — bet¬ 
ter  for  higher  In  contents — but  the  same  trend  is  exhibited. 
Overall,  the  comparison  suggests  that  the  EIW  model  is  a 
good  one. 

V.  CONCLUSIONS 

We  conclude  that  a  relatively  simple  analytic  model  for 
the  valence-band  structure  can  provide  a  reasonably  sound 
basis  for  calculating  parameters  of  direct  relevance  to  quan¬ 
tum-well  lasers  and  other  devices.  The  model  is  applicable 
for  arbitrary  strains  and  arbitrary  quantum  confinement  ex¬ 
cept  where  “accidental"  light-hole/heavy-hole  degeneracies 
occur. 
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We  idapt  '.he  standard  theory  of  the  free  acoustic  vibrations  of  an  inAnice  plate  of  homogeneous  iso¬ 
tropic  elastic  matter  to  the  corresponding  case  of  optical  vibrations.  Treating  nonpolar  maienai  .irst.  we 
show  that  the  effect  of  the  tree  surface  is  to  couple  LO  and  TO  modes,  and  we  demonstrate  the  e.xistence 
of  the  optical  analog  of  Rayleigh  waves.  Interface  and  guided  modes  are  both  present,  and  their  respec¬ 
tive  mode  patterns  are  derived.  In  polar  materials  the  coupling  between  LO  and  TO  is  different  because 
of  the  frequency  splitting  due  to  the  ionic  fields,  but  surface  modes  are  still  present.  Tais  result  contra¬ 
dicts  the  conclusion  of  the  hydrodynamic  model  that  surface  modes  do  not  e.'tist.  The  polar  character 
also  allows  the  e-tistence  of  sunace  polantons.  It  is  shown  that  the  standard  description  of  these  modes, 
which  neglects  the  elastic  properties  of  the  matenal.  ts  physically  invalid.  The  effect  of  the  free  sunace  is 
to  couple  sunace  polantons  and  LO  modes,  and  a  descnption  is  given  of  the  mode  patterns  that  may 
occur.  General  e.xpressions  for  energy  flux  are  given,  and  boundary  conditions  for  the  general  .'ase  are 
suggested.  This  treatment  goes  some  way  towards  reconciling  ihe  vanous  theoretical  models  of  phoncti 
confinement  that  have  been  advanced  recently. 


» 


» 


» 
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I.  INTRODUCTION 

■A  description  of  the  free  acoustic  vibrations  of  an 
infinite  piate  of  homogeneous  isotopic  elastic  matter  has 
been  available  for  over  a  century,*  but  a  corresponding 
description  for  optical  vibrations  has  not  been  given,  to 
the  author’s  knowledge.  Here  we  report  such  a  descrip¬ 
tion  for  nonpolar  material,  and  subsequently  for  polar 
matenal.  The  motivation  was  to  establish  a  consistent 
account  of  longitudinally  polanzed  (LO)  and  transversely 
polarized  (TO)  optical  modes  in  a  thin  layer  treated  as  an 
isotopic  elastic  and  dielectric  continuum  to  act  as  a  basis 
for  calculating  the  electron-phonon  and  hole-phonon 
scattering  rates  in  layered  semiconductors.  Here  we  limit 
attention  to  the  free-standing  plate  since  this  is  the  sim¬ 
plest  system  that  illustrates  the  basic  physics.  The  latter 
emerges  in  the  form  of  an  unavoidable  coupling  of  LO 
and  TO  modes  brought  about  by  the  presence  of  a  sur¬ 
face.  Such  a  coupling  between  acoustic  modes  (LA  and 
TAl  is  well  known,  but  its  significance  for  optical  modes 
has  not  been  widely  appreciated  hitherto. 

Models  of  confined  optical  modes  which  are  currently 
in  the  literature  achieve  diffenng  levels  of  sophistication. 
The  earliest  and  simplest  was  the  dielectric-continuum 
(DC  model  in  which  only  electrical  boundary  conditions 
were  used  to  determine  LO  mode  patterns.*'^  Such  a 
model,  clearly,  could  not  descnbe  confinement  in  nonpo¬ 
lar  matenal.  and  it  soon  came  into  senous  conflict  with 
the  predictions  of  linear  chain  models'*  through  its  lack  of 
boundary  conditions,  which  referred  to  mechanized  sta¬ 
bility.  A  more  successful,  but  still  overly  simple,  model 
described  LO  confinement  using  hydrodynamic  (HD) 
boundary  conditions.^  in  which  it  was  argued  that  since 
an  LO  mode  was  charactenzed  by  zero  electric  displace¬ 
ment  It  contained  no  electromagnetic  energy,  and  so  only 
purely  mechanical  conditions  'nad  to  be  satisfied  at  a 
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boundary.  The  HD  model  showed  much  closer  agree¬ 
ment  with  linear  chain  models,  though  it  was  soon  point¬ 
ed  out  that  HD  boundary  conditions  were  not  stnctly 
consistent  with  the  actual  ionic  motion  at  the  interface.* 
-Nevertheless,  the  HD  model  displayed  the  twin  virtues  of 
conserving  energy  and  establishing  mechanical  stability 
through  its  expansion  of  the  Bora-Huang  model  to  in¬ 
clude  dispersion.  On  the  other  side  of  the  balance  sheet. 

It  displayed  the  perceived  defect  of  allowing  tangential 
electric  fields,  and  with  them,  the  scalar  potential,  to  be 
discontinuous  at  the  interface. 

As  three-dimensional  models  of  the  lattice 
dynamics — usually  referred  to  as  microscopic  (M) 
models — became  more  sophisticated,  two  aspects 
emerged.  One  was  the  primacy  of  mechanical  conditions, 
and  the  other  was  the  added  complexity  introduced  by 
the  elastic  anisotropy  of  the  crystals  considered  (mostly 
CaAs  and  AlAs).  Huang  and  Zhu  attempted  to  recon¬ 
cile  the  DC  and  iV  models  by  ad  hoc  addition  of  scalar 
potentials  so  both  electric  and  mechanical  boundary  con¬ 
ditions  were  satisfied,  and  by  a  reinterpretation  of  off-axis 
modes.  Bechstedt  and  Gerecke*  attempted  a  similar 
reconciliation.  Both  pairs  of  workers  obtained  a  hybridi¬ 
zation  of  LO  guided  modes  with  Fuchs-KIicwer^  (FK) 
surface  polaritons.  The  additions  conceived  by  Huang 
and  Zhu.  besides  being  ad  hoc,  destroyed  the  orthogonali¬ 
ty  of  the  modes,  but  this  has  been  rectified  recently  by 
Haupt  and  'Wendler'®  who  have  gone  on  to  calculate 
electron-phonon  scattering  rates.  Some  unsatisfactory 
elements  in  all  of  these  approaches  have  been  pointed  out 
recently,*' 

A  continuum  theory  of  confined  optical  modes  should 
be  able  to  stand  on  its  own  and  be  internally  self- 
cor.sisteni.  U  should  also  be  applicable  to  nonpolar  and 
polar  materials.  Only  when  these  criteria  are  satisfied 
can  a  judgment  be  made  concerning  the  validity  ot  apply- 
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ing  it  to  a  real  system  composed  of  atoms.  The  HD  mod¬ 
el  sacisdes  these  criteria  by  limiting  itself  strictly  to  LO 
modes  and  seeing  any  discontinuity  of  the  scalar  potential 
that  may  arise  as  resolvable  only  on  an  atomic  scale.  The 
models  mentioned  above  stan  from  the  basis  of  the  DC 
model,  which  is  clearly  not  self-consistent.  In  what  fol¬ 
lows  we  develop  a  continuum  model  that  goes  beyond  the 
HD  model  in  that  it  treats  TO  modes  and  the  LO-TO  in¬ 
teraction  and  which  rediscovers  the  HD  model  as  a  spe¬ 
cial  case.  It  bases  itself  on  the  firmly  established  theory 
of  elasticity. 

Our  approach  is  straightforward.  We  first  identify  the 
dispersion  of  optical  modes  with  the  macroscopic  elastic 
properties  of  the  matenal.  This  immediately  allows  us  to 
take  over  the  theory  established  for  acoustic  modes  in 
nonpolar  matenal  and  apply  it  to  optical  modes.  In  do¬ 
ing  so  we  discover  the  optical-mode  analog  of  Rayleigh 
waves.  The  transition  to  polar  material  is  effected  simply 
by  taking  into  account  the  splitting  of  LO  and  TO  fre¬ 
quencies  and  introducing  a  new  mode  —  the  surface  po- 
lariton.  Throughout,  we  assume  that  the  medium  is 
elastically  isotropic,  which  allows  us  to  keep  LO  and  TO 
modes  distinct  and  allows  us  to  onentate  Canesian  a.xes 
irrespective  of  crystallographic  direction.  We  take  the  x 
a.xis  parallel  to  the  direction  of  propagation  in  the  plane 
of  the  plate  and  the  r  a.xis  to  be  perpendicular  to  the 
plane  of  the  plate.  Modes  with  displacements  in  the 
plane  of  incidence  i.rc)  are  the  LO  mode  and  the  p- 
polarized  TO  mode  lor  just  p-TO);  the  mode  with  dis- 
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FIC.  1.  Coordinate  system  and  mode  designation.  Depicted 
are  the  propagation  and  polarization  directions  of  the  plane- 
wave  components  that  together  with  their  reflections  make  up 
waves  with  a  standing  component  in  the  c-direction.  which 
propagate  along  the  x  a.sis. 


placement  at  nght  angles  to  the  plane  of  incidence  (i.e.. 
along  the  y  direction)  is  the  r-polarized  TO  mode  (or  r- 
TO)  (Fig.  1). 

II.  OPTIC.\L  PHONONS 

The  dispersion  relation  for  long-wavelength  optical 
phonons  in  nonpolar  matenal  can  be  wntten  as  follows: 

(tu-  — a>5'!uik)=  — iif(k)uik.  .  (1) 

where 


Ck,k, 

Cx.kx 


Ckxk, 

Ak’-rBikp-^k'] 

Ck.k. 

•  7 


Ck,k. 

Ck.k. 

7  • 

Ak;-Bik;^k-) 


(2) 


is  a  3X3  matnx.  u(k)  is  the  relative  displacement,  and  Og 
is  the  frequency  at  k  =0.  We  take  if(k)  to  be  identical  to 
the  matna  for  acoustic  modes,  from  which  it  follows  that 

/I  =l’£  =C|, /p  ,  B  =Vf  =  C^x/p  , 

(3) 

C  =iC^2^C^^)/p  . 

where  U[_,  u-p  are  the  velocities  of  LA  and  TA  modes,  c,,. 
etc.  are  the  elastic  constants,  and  p  is  the  mass  density. 
Thus  the  dispersion  relations  for  LO  and  TO  modes  are 

(iy'=u!j  —  ulk~  LO 

=  — (sandp),  (4) 

where  k-  =  k;  ~  k'- k:  for  all  directions  provided  the 
matenal  is  elastically  isotropic,  viz.. 

C||— •C.;-"2Ct4— 0  (5) 


equivalent  strains  are  given  by 


5.1 


dy 


3u_, 

3y 

dz 

3w, 

dx 

du.. 

ay 

(6) 


We  will  assume  this  to  be  the  case. 

Dispersion  arises  as  a  consequence  of  the  elastic 
stresses  produced  hv  a  traveling  optical  wave  The 


The  minus  signs  arise  as  a  consequence  of  the  oui-ot- 
phase  vibration  of  the  two  atoms  in  the  unit  cell.  The 
strcs.ses  are  given  by  the  usual  set  of  equations: 
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r,  =C|i5| -rC|j(S; -r^j  )  , 

(7) 

r,=2c«5,  . 
r<  =  ;c«s,  . 
r,=:c«5, . 

Equations  (6i  and  (7)  allow  us  to  calculate  dilatational 
and  shear  stresses  producni  by  a  travelling  optical  wave. 

The  assumption  of  isotropy  allows  us  to  decompose 
any  displacement  u  into  a  longitudinally  polarized  part 
and  a  transversely  polarized  part  U7- 


FIG.  2.  Schematic  dispersion  relationship  of  (a)  nonpolar  and 
lb)  polar  material. 


ci}'  —  ct}Q  —  ulllc‘-Tkl)  =  cj^  —  uiik;  —  kf)  .  (14) 


U  =  Ur— Ut-  (8) 

with  U£  and  Up  defined  by'* 

VXUr  =0  . 

V-u-=0  .  (9) 

Of  course,  Uj_  and  Uj-  must  have  the  same  time  depen¬ 
dence.  i.e,  the  same  frequency. 

III.  nonpolar  M.ATERLAL 


It  is  convenient  to  define  a  “plate”  velocity  u  and  a  pa¬ 
rameter  s  such  that 

cjl—cj'=u‘k~  and  s—u'/uf  ,  (15) 

from  which 

kr^iys  —  Dk^,  k^  =  i.s  —  \)k^,  (16) 

where  y  =iif /ur  =1:14/0,,. 

As  indicated  in  the  .Appendi.x.  the  boundary  conditions 
can  be  satisfied  provided  the  following  equation  is  true: 


The  boundary  condition  which  must  be  satisfied  when 
the  surface  is  free  is  that  the  dilational  stress  perpendicu¬ 
lar  to  the  surface  and  shear  stress  across  the  surface  van¬ 
ish.  The  situation  is  simplest  for  r-poianzed  TO  modes. 
A  solution  is 

,  (10) 
which  satisfies  V-u=0.  It  must  also  satisfy 


(qr  —pt)-%ink-pL  sinx^L 

-2qrp;[cosiI:, -\v)L -1]=0  ,  (17) 

where 

p  —  yis~l),  q^ly'^s  —  l, 

_  (18) 

r=^lVys  -  1.  t  =S  — 2  . 


T,^Q,  74=0,  75=0 


(II) 


at  r  =0  and  z—L.  For  this  wave  the  stresses  73  and 
are  zero  everywhere.  To  obtain  Ti  =0  we  must  have 

ik  X 

u^—e  cQsk.z  ,  /c,L=n?r  ,  (12) 

where  n  is  an  integer.  In  this  case  no  mixing  with  other 
modes  is  required. 

The  surface  couples  LO  and  p-polarized  TO  modes. 
Thus  we  take 


u^—k^e  '  (  Ae  ~  -r Be  ‘-  ) 

—  Kj-e  '  ( Ce  '  -rDe  ^  )  , 

,  tk  X  ik.;  tk,: 

u.=ki_e  ‘  Ae  •  -Be  *  ) 

,  iJcppiS  _  _  -k,Z 

—  k,e  '  iCe  •  -De  '  )  . 


113) 


This  admits  of  two  waves.  One  is  a  guided  mode  consist¬ 
ing  of  coupled  phased-matched  LO  and  TO  waves,  the 
other  is  a  surface  mode  consisting  of  coupled  evanescent 
LO  and  TO  waves.  The  latter  is  the  optical  analog  of  the 
Rayleigh  wave  in  acoustics.  Some  mode  patterns  are  de¬ 
scribed  below. 


A.  Guided  modes 

These  waves  have  phase-matched  components  with  the 
following  wave  vectors: 

k[_L=n[_:r,  kj-L=n-rrr,  ,  O')) 

where  n^,  rij-,  and  m  are  ail  integers.  The  displacement 
in  the  z  direction  can  be  descnbed  by  sine  or  cosine.  The 
sine  solution  is 


We  distinguish  the  z  components  of  the  LO  and  TO  wave 
vectors  k^_  and  kj.  respectively.  Equation  i9)  is  satisfied 
tor  the  individual  LO  and  TO  components,  and  it  is  clear 
that  the  condition  7^=0  holds  everywhere.  The  wave- 
vector  component  along  the  surface  musi  be  common  to 
both  components.  They  must  also  have  the  same  fre¬ 
quency,  from  which  '  Fig  2) 


u^=2.Ie  '  \k^zosk^z  —kj-^zoskyz 

1  9 


u.—liAe  ‘  ki^i\nki_z  —k^—'iinki-z  \  , 

9  1 


and  the  cosine  solution  is 
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FIG.  3.  Relation  between  and  k,  for  guided  modes  in 
nonpolar  material  with  L  is  the  width  of  the  plate  and 

2m  int_  ^ki_L/~,  iir^kfL/ir). 


u,=2i/4e'*'^  \krSink,  2  —kr—sink-rz 

X  [  ‘  t  ‘ 

u.  =  lAe''''  \  k,^zo%kiZ  —k^—ZQ%kjZ 
I  ^ 


(21) 


LO  and  TO  modes,  which  arises  because  their  frequencies 
lie  close  together.  The  frequency  of  the  TO  mode  is  the 
higher  as  a  consequence  of  the  smaller  velocity  of  TA 
modes.  One  consequence  of  this  is  that  the  smallest  TO 
vector  component  in  the  :  direction  that  is  allowed  is  one 
where  =  1  and  /i7-  =  3,  to  satisfy  Eq.  (19).  Confinement 
eliminates  n7-=l  and 


B.  Surface  modes 


No  surface-wave  solutions  exists  for  s-polarized  TO 
modes,  but  Eq.  (17)  allows  both  k[_  and  k-p  to  be  imagi¬ 
nary.  For  brevity  we  quote  mode  patterns  only  for 
L  —  00  and  L  —>0,  since  simple  expressions  are  obtained 
in  these  limits. 

For  L  —  00 ,  the  antisymmetric  solution  is,  with 


«,=2Ac 


k^e  “^^^coshoilx  — L/2) 


«.  =  2^e  ' 


‘cosfaaT-^r  — 1./2) 


—  iaie  “‘■^^sinhai_(z —L /2) 


(22) 


-rk^e  °’'^"*sinhcric  — L /2) 
9 


Figure  3  shows  the  relationship  between  ki_  and  k,  when 
y=4-.  The  spectrum  can  be  seen  as  consisting  of 
branches  defined  by  the  integer  m.  Figure  a  depicts  the 
corresponding  dispersion.  All  branches  are  associated 
with  a  velocity  {actor  j  » 3,  corresponding  to  a  velocity  v 
close  to  0^.  Note  that  the  choice  of  y  (  =c«/cii  •  to  be 
one-third  in  these  examples  is  reasonable  in  view  of  the 
actual  ratios  of  elastic  constants  found  among  semicon¬ 
ductors. 

These  guided  waves  are  a  profound  hybridization  of 


and  the  svmmetric  solution  is 


u,=2^e 


tt.  =2.^e 


-a,L/Z 


kjB  *  ‘sinhor  (r  —  L /2) 


—ia-j-^e  ’sinha7-(r  — 1./2) 

9 

ioje  “^^"^cosha^^lz  —  L/2) 

—  Icj^e  "*’’^^cosha7-(z  — L/2) 
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(23) 
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FIG-  ^  Dispersion  for  guided  modes  in  nonpolar  maienal 
with  y  =  y  and  £j^/t‘r-  =  20.  (G  LO  modes;  A  TO  modes.l 

I  fl .  ^  kj_L/rrt 


In  this  limit  Eq.  (17)  reduces  to 

9r-rpr=0,  (24) 

which  leads  to  a  cubic  equation  for  s,  the  velocity  factor, 

viz., 

s^-8j--!-8(3-2r)s-l6(I-r)=0  .  (25) 

This  equation  is  just  that  found  in  the  theory  of  Rayleigh 
12 

waves. 

With  r  ~  T  this  reduces  to 

(j-4)(3s--l2s  +  8)=0  .  (26) 

The  solutions  are,  therefore. 


1 

1-r-L 

,  i  =2  i 

1 

V'S 

Only  the  last  is  consistent  with  having  k^  and  kj  imagi¬ 
nary.  We  thus  obtain  the  relationship  between  a^,  Ot- 
and  A;,,  viz.. 
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FIG.  5.  Mode  patterns  for  nj-  =  3  in  nonpolar  materi¬ 

al.  (Dashed  line,  LO;  dotted  line  TO;  continuous  line.  LO  — 
TO  .) 


ar  = 


1/2 


1/2 


(28) 


Note  that  the  LO  component  drops  off  from  the  surface 
more  rapidly  than  does  the  TO  component. 

For  L  —0  only  the  antisymmetric  mode  survives,  the 
mode  pattern  being  that  of  Eq.  (22)  with 
exp(  —  a  j-L  /2 )  =  1 .  In  this  case  we  find 

r=2(l-2r),  (29) 

from  which 

,  ar  =  i4y-l)‘^’-lc,  .  (30) 

This  completes  our  description  of  the  allowed  modes  of 
a  nonpolar  plate.  A  depiction  of  mode  patterns  is  given 
in  Figs.  S  and  6. 


IV.  POLAR  MATERIAL 

When  the  two  atoms  in  the  unit  cell  are  oppositely 
charged  their  oscillations  are  accompanied  by  long-range 
electric  fields  that  modify  the  elastic  restoring  forces  and 
alter  the  frequencies  of  LO  and  TO  modes.  The  disper¬ 
sion  relation  for  long-waveiength  optical  modes  becomes 
modified  along  the  following  lines. 

(&)-— £j^')u(fc)  =  [  — H (A:i-r(t)^  — tu^)/(A:)iu(k)  ,  (31) 


FIG.  6.  Mode  patterns  for  surface  waves  in  nonpolar  materi¬ 
al  ^designations  as  in  Fig  5). 


where  tuj-  are  the  LO  and  TO  zone-center  frequen¬ 
cies.  and 


k  fc 


kjky  k^k, 

k;  k^.k. 
k.ky  k: 


(32) 


The  LO  frequency  is  shifted  upwards  relative  to  oiq  of 
Sec,  III.  The  main  effect  is  this  shift  in  frequency.  For 
simplicity  we  will  continue  to  regard  the  material  as 
elastically  isotropic  with  effective  elastic  constants 
modified  by  the  polar  fields.  Consequently  we  can  em¬ 
ploy  the  analysis  of  the  previous  section,  merely  changing 
the  definition  of  certain  terms. 

Thus  we  maintain  the  distinction  between  irrotational 
and  divergenceless  modes  according  to  Eq.  (9).  It  is  con¬ 
venient  to  modify  the  definition  of  the  plate  velocity  u  as 
follows; 


01^  —  =  (33) 

from  which  it  follows  that  Eq.  (16)  becomes 

=(ys  —  1  )fc*  ,  kf  =  is—\)k'  —  kl,  (34) 


where 


(35) 


The  secii'ar  equation.  Eq.  ipi,  remains  valid,  but  now 
Eq.  (IS)  becomes 
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p=T(i— 2),  9  =  27[(i -l)fc;-k5l‘'’/fc,  , 

_  ,  ,  (36) 

r^2'^rs—l,  t  ={s —D—ik^/k;)  . 

The  splitting  of  the  LO  and  TO  frequencies  profoundly 
modifies  the  ability  of  the  surface  to  couple  these  modes 
together.  Over  the  frequency  range  coupling  of 

guided  LO  inodes  can  only  be  accomplished  with  rapidly 
varying  TO  evanescent  modes,  and  for  only 

short-wavelength  LO  modes  can  be  involved.  Because  of 
these  consequences,  a  description  of  coupled  modes  tends 
to  stretch  any  continuum  model  to  its  limits  of  validity 
and  perhaps  beyond  it.  With  the  caveat  we  will  proceed 
first  to  describe  the  guided  modes. 

.\.  Guided  modes 

It  may  be  noted  that  our  results  for  r-polarized  TO 
modes  remain  unaffected  by  the  ionic  nature  of  the  ma¬ 
terial.  and  so  we  turn  immediately  to  the  coupling  be¬ 
tween  guided  LO  and  evanescent  TO  modes.  Taking 
Ic^asijlcg  and  assuming  for  simplicity  that  kg  is  large,  we 
obtain. 


=  2.*le  liJc^siTikfX 


—  — cost k.  Lie  i  , 


u,  =  2.-le'  '  ik,  cosk,  r 


—  k,y[e  cos(k,L)e  ^“e  ’"]  j 

with  kj_L  an  integer  (Fig.  7).  In  this  coupling 

the  contribution  of  the  TO  mode  to  the  dilational  stress  is 
negligible,  but  the  presence  of  the  TO  mode  is  still  vital 
to  cancel  out  the  shear  stress.  If  shear  is  neglected,  no 
TO  mode  need  be  involved,  and  this  is  e.sactly  the  impli¬ 
cit  assumption  underlying  the  HD  model,  whose  bound¬ 
ary  condition  7-u=0  at  the  surface  is  e.tactly  the  same  as 
our  condition  7^=0  when  y  =0. 

It  1$  interesting  to  observe  that  the  mode  pattern  of  Eq. 


(37)  satisfies  the  electromagnetic  boundary  conditions. 
The  electric  and  magnetic  fields  associated  with  the  TO 
mode  are  vanishingly  small,  so  they  do  not  enter.  The 
LO  mode  has  zero  magnetic  field  and  zero  electric  dis¬ 
placement  D.  Thus  O.  =0.  and  since  k[_'=nir,  the 
tangential  electric  field  vanishes  at  each  surface.  The 
continuity  of  D.  and  £,  is  then  satisfied  by  having  no 
fields  in  the  vacuum. 

Tne  coupling  between  guided  LO  and  guided  TO  for 
loSiiif  gives  the  mode  patterns  already  given  in  Eqs.  (2) 
and  (21)  but  now  with  (Ot  replacing  C0[_  in  Eq.  (33), 
kj-L  k^L  —n^TT  a.nd  — .e_r  =2m. 

kv  =  yk5 -rk,'( j'f  —  1)  ,  kr  =  'j  — Dk"  (38) 


p  =  -/.:k5^(s-2)k  =  ]/k^  c  =  :;,Vi-i  , 
r=il/k^)V'  ykl^k'iys  -  ,  t=\s—l)  . 


B.  Surface  modes 

There  are  still  no  surface  s-polarized  TO  modes,  but  in 
general,  surface  waves  consisting  of  mixed  LO-p- 
poiarized  TO  modes  do  e.tist.  and  their  mode  patterns  for 
L  —  =0  are  described  by  Eqs.  22)  and  (23)  with  p,  q, 
df  —:k,.  aj-  =  ikj-,  given  by  Eqs.  (3a)-!36).  Specific 
magnitudes  must  satisfy  Eq.  (2-  .  which  now  leads  to  a 
quanic  for  the  velocity  factor  s- 

r*-2(-t-u)s"-i3(3-2r)-12J  -a-r* 

-[  16(l-r)-^Sa(3-2v’--iG->  -4<j==0  ,  (40) 


a  -kl/k;  . 

Note  that  when  a— 0  Eq.  (40)  reduces  to  the  nonpolar 
equation  [Eq.  (25)].  In  the  limit  a  —  » ,  the  solution  is 


independent  of  y,  from  which 

al=(l—2y)k-  ,  af=k^-k-  .  (42) 

In  this  limit  solutions  e.'tist  provided  y  E  t.  Waves  propa¬ 
gate  with  surface  phase  velocity  equal  to  The  an¬ 

tisymmetric  solution  is 


u,  ==2/4k,e'*‘^  e  “‘•^‘*^cosha(_(r  — L/2) 


—  2-r^e  *’^'coshk,(z  — L/2)  , 

kg 


u.~2iAai_e'‘‘''  — e  "’'■^^^sinhoj^tr  —  L /2) 


FIG.  7.  Mode  patterns  for  '1^  =  1  in  polar  matenal  (with 
k,=ki_)  (designations  as  in  Fig.  5).  .Mode  patterns  for  surface 
waves  are  approximately  those  of  Fig.  6  in  the  absence  of  the 
TO  component. 


-1-2  — e  ’  'sinhkglr  — L /2) 


and  the  svmmetric  solution  is 


=  n„(u))  -■ 


/*f(  I  — a  ) 
■  Mfitu)- 


where  F  is  the  steady-sute  power  density  input. 

The  time  constants  characterizing  the  response  are  tj-. 


are  strong  enough  to  establish  a  common  electron-hole 
temperature  there  is  no  need  to  consider  well  capture 
explicitly — it  occurs  as  rapidly  as  the  thermalization  pro¬ 
cess. 

Energy  relaxation  of  the  strongly  coupled  electron-hole 
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The  existence  of  this  optical-mode  analog  to  Rayleigh 
waves  in  polar  material,  as  well  as  its  appearance  in  non¬ 
polar  material,  depends  upon  the  shear  modulus  being 
nonzero.  If  y=0.  Eq.  (42)  shows  that  ai=k;.  In  this 
case  the  LO  component  has  the  propeny  7Xu=0  as  well 
as  VXu=0.  Such  a  mode  must  have  zero  amplitude, 
which  is  consistent  with  the  prediction  of  the  HD  mod¬ 
el'^  that  no  surface  modes  are  possible  in  a  free-standing 
plate.  However,  this  conclusion  derived  from  the  HD 
model  appears  to  be  incorrect  for  real  elastic  solids,  even 
though  Eqs.  !43)  and  (44)  show  that  the  amplitudes  of  the 
TO  components  are  vanishingly  small.  We  note  that  in 
this  limit  la —  *  i.  the  solutions  obtained  above  are  valid 
for  ail  plate  thicitnesses  of  interest.  However,  we  also 
note  that  the  tangential  component  of  the  eiectnc  iieid. 
which  is  proportional  to  u,,  is  not  zero  at  the  surface, 
and  so  conventional  eiectnc  boundary  conditions  are 
violated. 


C.  Surface  polaritoos 

One  well-known  consequence  of  polarity  is  that  it  al¬ 
lows  electromagnetic  waves  to  directly  couple  with  TO 
modes  to  form  sunace  poiantons.  as  described  by  Fuchs 
and  Kliewer.’  We  will  refer  to  these  as  FK  modes.  Their 
dispersion  is  depicted  in  Fig.  8.  For  large  fc,  the  mode 
patterns  are  of  the  following  forms:  antisymmetnc, 

u,  =  2  .-le  '  ‘e  "  coshA,  ( r  -  £  /2 )  , 

(45) 

u.  =  —2iAe'^’'e  *'^^^^c,sinh/c,(r  — £ /2) 
and  symmetric. 


FIG  S  Fiichs-Kliewer  Jispersion. 


u,  = -T-2/.-le'*'%  *'^''^c,sinhlc,(z  — £/2)  , 
u.  =  — 2/le"*‘%  *‘^‘*k,cosh^,(z  —  £/2)  . 


(46) 


When  It,  is  large  compared  with  the  wave  vector  of  light 
at  the  same  frequency,  these  solutions,  share  the  same 
frequency  that  lies  between  and  and  is  given  by 


oi'r  , 


(47) 


where  e,,  are  the  static  and  high-frequency  permittivi¬ 
ties  and  e.j  is  the  permittivity  of  the  vacuum.  For  simpli¬ 
city  in  our  discussion  we  will  assume  that  this  condition 
is  met. 

Two  comments  can  be  made.  First,  it  would  appear 
from  Eqs.  (45)  and  (46)  that  we  are  dealing  with  a  null 
mode  since  both  V-u  and  7Xu  are  zero.  This,  however, 
is  not  really  so  since  the  decay  constant  in  hyperbolic 
terms  is  only  approximately  equal  to  the  wave  vector  It, . 
In  fact  the  modes  are  TO  modes  with  finite,  but  in  the 
present  approximation,  vanishingly  small  rotation,  which 
describes  the  rate  of  change  of  magnetic  field.  The  solu¬ 
tions  quoted  are  those  obtained  neglecting  retardation 
(i.e..  assuming  the  velocity  of  light  to  be  infinite).  They 
match  smoothly  with  an  evanescent  electromagnetic 
wave  in  the  vacuum. 

The  second  comment  is  that  these  solutions  violate  the 
condition  that  the  elastic  stress  vanish  at  the  sunace. 
Fuchs-Kliewer  modes,  conventionally  descnbed.  are 
therefore  unphysical.  This  rather  important  defect  can 
be  rectified  by  coupling  with  an  LO  guided  mode  of  the 
same  frequency,  as  was  done  in  Sec.  IV  A.  For  £  —  oe 
only  cosine  waves  couple,  and  the  corresponding  mode 

«,  =  2.-lc'*''*[/I:,sinA:i,r 

•^ik^ie  •  —coski^Le  '  e  )]  , 


FIG.  9.  .Mode  patterns  for  hybndized  FK  modes  in  polar 
material  (with  KL—-irr.  'Dashed  line.  LO:  dotted 

line.  FK;  continuous  line.  LO  —  FK.;  .Note  that  this  combina¬ 
tion  of  FK  and  LO  is  possible  only  if  the  two  modes  have  the 
same  frequency.  In  general  only  certain  values  of  k,  arc  al¬ 
lowed. 
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with  <[_L  ='iT.  satisfying 

“  u}~  “ '  )  (*^9^ 

This  mode  pattern  Fig.  i9'  satisfies  all  elastic  and  elec¬ 
tromagnetic  boundary  conditions.  Note  that  when  n  is 
an  odd  integer,  corresponding  to  an  asymmetric  guided 
mode,  the  coupling  is  :o  the  asym.-netric  FK  ode.  and 
when  n  is  an  even  integer  both  components  are  sym- 
metr.c.  In  other  woros.  like  parity  hybndizes  with  like 
parity. 

.\  further  discussion  of  hybndization  of  FK  and  LO 
modes  taking  into  account  FK  dispersion  will  be  reported 
soon  in  the  future. 


the  medium  cannot  be  ignored.  The  important  parame¬ 
ters  are  the  wave  vector  k,j.  Eq.  i35).  and  the  dimension 
of  the  unit  cell  j,.  So  for  a  continuum  theory  to  be  valid 
we  e.spec: 

'.kajQi  <  1  .  (50) 

For  Ga.As,  =  1. 

12)  The  relation  between  HD  boundary  conditions  and 
those  arising  from  the  linear-chain  model  was  illuminated 
by  .Akero  and  .Ando.’  When  the  force  constants  on  either 
side  of  the  interface  are  the  same,  and  when  the  interface 
IS  taken  to  lie  midway  between  atom  A  and  atom  3.  the 
condition  for  mechanical  siaoiiity  simply  entails  the  con¬ 
tinuity  of  the  displacements  a  ,  and  uj,  viz.. 


V,  DISCUSSION 


Our  e.stension  of  classical  elasticity  theory  to  include 
optical  modes  has  allowed  us  to  described  the  mi.xing  of 
LO  and  TO  modes  at  a  free  surface  and  to  obtain  the 
confined  phonon  spectrum.  In  panicular.  it  has  allowed 
us  to  demonstrate  the  e.xistence  of  surface  optical  modes, 
which  are  analogous  to  Rayleigh  waves.  Hitherto  only 
HD  theory  has  predicted  interlace  modes,  but  only  at 
internal  surfaces  and  only  for  LO  modes.  DC  theory  and 
the  more  sophisticated  versions  mentioned  in  the  Intro¬ 
duction  do  not  predict  interface  mooes  other  than  FK 
poiantons.  so  our  conclusions  here  are  quite  new.  Decid¬ 
ing  the  ontology  of  surface  waves  is  of  some  imponance 
for  the  electron-phonon  interaction,  since  under  some  cir¬ 
cumstances  surface  modes  can  interact  strongly  with  elec¬ 
trons.'*  For  the  same  reason  it  is  important  to  ascertain 
the  strength  of  the  interaction  with  FK  polaritons.  It  has 
been  argued  elsewhere  '  that  such  an  interaction  has  the 
nature  of  a  magnetic  ^  A  pi  rather  than  of  an  electrical 
ied)  interaction,  but  clearly  an  even  greater  reappraisal 
of  the  FK  interaction  is  now  necessary  as  a  result  of  our 
demonstration  of  hybndization  of  FK  and  LO  modes. 

The  foregoing  treatment  has  included  nothing  but  the 
simplest  cases  and  has  merely  sketched  a  rough  outline  of 
the  topic.  An  extension  to  the  quantum-well  situation  is 
underway,  where  continuity  of  stress  replaces  the  vanish¬ 
ing  of  stress,  and  continuity  of  energy  flow  is  entailed.  It 
is  clear  that  a  self-consistent  picture  of  confined  modes 
can  emerge  from  this  program,  so  it  is  of  interest  to  ask 
how  valid  such  a  picture  can  be  as  a  description  of  optical 
waves  in  real  matter.  There  are  at  least  four  problems, 
and  there  may  be  more.  These  have  to  do  with  the  fol¬ 
lowing:  (1)  The  spatial  vanation  of  the  mode  pattern,  (2) 
the  relation  of  relative  displacement  i  ui  to  ionic  displace¬ 
ment.  (3)  the  elastic  anisotropy  of  real  crystals,  and  (41 
The  electncal  boundary  conditions.  We  wiil  comment  on 
each  in  turn. 

M)  It  has  already  been  remarked  in  Sec.  IV  that  the 
polarity  of  the  medium  typically  forced  one  of  the  com¬ 
ponents  of  a  LO-TO  hybndization  to  be  rapidly  varying, 
the  more  so  the  stronger  the  poianty.  i.e..  the  more 
disparate  the  LO  and  TO  zone-center  frequencies.  When 
the  spatial  variation  of  ;he  mode  function  becomes 
significant  over  a  primitive  unii  cell,  the  atomic  nature  if 
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(51) 


(Note  there  are  no  shear  stresses  here.;  The  parameter  j, 
is  now  the  interatomic  spacing,  assumed  to  be  identical  in 
both  matenals.  and  it  is  also  assumed  that  u  ^  and  Ug  are 
iescr.bable  by  slowly  varying  envelope  functions. 
Opticai-mode  relative  displacement  a  for  long  wave¬ 
lengths  IS  related  to  u  ^  and  u,  os  follows: 

u.  = - - - '  u,  = - - -  ,  (52; 

*  I-i.V/^/.Vfa;  ^  l-'.Vfa/.V/.,  i 

where  \f  ^  .Mg  are  the  atomic  masses.  Substitution  into 
Eq.  (51)  yields  the  connection  rule 


(53) 


where 


T 


r 


4(3?  ~R  ■')  T-'i?  -R  ■') 
R-R-'  ^:R-R~') 


(54) 


Here  r  =  ^  r  “ ' ,  r,  =  ( .Vf  ,  /.Vf ,  i  ‘ R  =  % ,  /r„  j.  This 

connection  rule  satisfies  continuity  of  energy  Sow.  The 
corresponding  connection  rule  for  the  continuum  theory 
is 


(55) 


which  agrees  with  microscopic  theory  only  when  R  =1. 
I.e..  the  mass  ratios  are  common.  In  reality,  the  interface 
is  never  as  precisely  located  as  conventional  microscopic 
theory  assumes,  and  the  position  and  properties  of  the  in¬ 
terface  aifect  the  olf-diagonal  elements  in  Eq.  (54).  Con- 
sequentlv,  boundary  conditions  relating  relative  displace- 


tively.  The  new  factors.  /.  r  and  p  are  all  dimensionless. 
The  first,  /.  may  be  termed  the  injection  hot-phonon  fac¬ 
tor,  it  is  a  measure  of  the  energy  injected  into  the  phonon 
system  and  hence  of  the  increase  of  phonon  number  m  a 
.mode.  The  second  factor,  r.  is  the  recombination-heating 


The  departure  from  ideal  behavior  is  often  described  in 
terms  of  a  nonlinear  gain  coefficient  such  that  g^go/i  1 
The  parameter  is  found  from  the  ratio  of  the 
damping  rate  y  to  the  square  of  the  resonance  frequency, 
f\.  This  ratio,  usually  named  K.  is  proportional  to  e,j  at 
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ment  and  siastic  stresses  in  a  continuum  theory  may  not 
be  a  major  source  of  error. 

1 3'  Elastic  anisotropy  adds  signtncant  complexity.  In 
directions  other  than  the  m.-’jor  crystallographic  ones 
there  is  no  ciean  distinction  between  LO  and  TO.  and  the 
degree  of  surface-induced  hybridization  depends  upon 
airection  foiio'.ving  the  directional  dependence  of  disper¬ 
sion.  However,  orientationai  etfects  of  this  sort  are  not 
likeiy  to  aud  anyt.-.ing  quaiitativeiy  new  to  the  picture  ob¬ 
tained  on  the  basis  of  an  isotropic  modei.  In  any  case,  a 
straightforward  appiication  to  the  problem  o.'  the 
eiectron-phonon  interaction  will  require  some  angular 
averaging.  .Anisotropy  is  therefore  not  likely  to  lead  to 
signincant  probiems. 

•i)  There  is  no  question  but  that  a  transversly  polar¬ 
ized  wave  of  eiectromagnetic  character  must  satisfy  elec¬ 
tromagnetic  boundary  conditions.  The  same  is  true  of 
static  fields.  The  question  however  is  open  concerning 
the  applicability  of  electr.c  boundary  conditions  to  LO 
Seids.  We  noticed  that  the  surface  optical  modes  de- 
scr.bed  in  Sec.  IV  B  did  not  satisfy  these  conditions.  Nor 
do  linear-chain  models  in  cenain  situations.  Indeed,  it 
was  this  problem  that  brought  the  DC  modei  into  clear 
conflict  with  microscopic  modeis.  The  same  question 
anses  in  the  conte.xt  of  piasma  waves.  The  special  prop¬ 
erty  of  LO  waves  is  that  they  have  no  eiectncal  energy 
because  the  permittivity  is  zero.  Electrical  neids  under 
this  circumstance  cannot  e.xist  in  the  absence  of  ionic  po¬ 
larization.  .At  the  sur'ace  of  an  ionic  solid  the  eiectnc 
heid  must  therefore  drop  to  zero  over  a  distance  of  the 
order  of  atomic  dimensions.  If  this  is  true,  there  is  as 
much  need  for  a  continuum  theory  to  address  the  issue  as 
It  would  be  for  it  to  address  the  problem  of  continuity  of 
interatomic  forces,  or  indeed  the  continuity  of  matter.  .A 
surface  or  an  intenace  is  already  a  discontinuity  of 
matter  pius  its  propenies.  including  polanzation  and  as¬ 
sociated  fields,  which  suggests  that  it  is  unnecessary  to 
aad  nonmechanical  boundary  conditions.  If  this  view  is 
adopted,  it  becomes  clear  that  the  DC  modei.  insofar  as  ii 
refers  to  the  confinement  of  LO  modes,  is  fundamentally 
flawed,  and  discontinuities  of  LO  fields  are  to  be  expect¬ 
ed.  On  the  other  hand,  if  this  is  denied,  then  one  of  the 
things  that  follows  is  that  the  surface  modes  of  Sec.  IV  B 
cannot  survive.  This  implies  th;  t  surface  modes  exist  in 
nonpolar  matenal  but  not  in  polar  matenal  irrespective 
of  the  strength  of  poiar.ty.  provided  the  latter  is  nonzero. 
We  appear  to  have  exchanged  one  type  of  discontinuity 
:or  another.  Note  one  final  comment.  The  necessity  for 
any  boundary  condition  must  be  justified  by  arguments 
based  on  the  physics  of  the  situation.  Hitheno  there  ap¬ 
pears  :o  be  no  such  argument  specifically  directed  at 
=  LO  modes  for  ensuring  the  continuity  of  the 
tangential  component  of  the  eiectnc  fieid. 

It  turns  out.  however,  that  both  elastic  and  electromag¬ 
netic  boundary  conditions  can  be  satisfied  by  invoking  a 
-.cneme  of  triple  hyoridization  involving  LO.  TO.  and  FK 
modes  and  disiinguisning  carefully  between  scalar  and 
vector  potenuals. 

'•''e  conclude  with  some  comments  on  energy  flow. 
General  espressions  for  the  energy  flow  m  the  c-direction 
ma;.  oe  written  as  .'oilows. 


Soiol 

S.~ - — ^(u;V-u-u.V  u‘)  LO  , 

Zi 

pcovi 

= - [u;(VXu)„-u,(VXu*)J  TO.  '5o) 

where  p  is  the  reduced  density.  The  solutions  we  have 
derived  all  obey  5.  =0.  In  the  general  case  a  basic 
boundary  condition  must  be  the  continuity  of  S..  Simi¬ 
larly  we  must  have  continuity  of  the  stress  components. 
These  considerations  suggest  that  if  the  c  axis  is  perpen¬ 
dicular  to  the  interface,  the  general  boundary  conditions 
to  be  satisfied  for  LO  and  TO  modes  are 

5.-.  Tj-.  r,  continuous  .  (57) 

In  general,  we  must  add  the  usual  electromagnetic 
boundary  conditions. 
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APPENDIX 


The  displacement  of  Eq.  13)  must  satisfy  the  following 
boundarv  conditions  at  c  =  0  and  L: 


1  ou. 

OU. 

'5  <^11  j 

T  —  — ./• 

ax 

3u. 

or 

‘  5  ^*4 

3z  dx 

=0  , 


This  entails  the  vanishing  of  the  determinant 

-q  -q 


=0  . 


(Al) 


(A2) 


{A3) 


P  P 
r  —r  [ 

|p/t.  PSl  -R/t  -RSt 

i/'/t.  --Rr 

where 

,  q=2ykj./k^ 

r-=Zki_/k^  .  i  =1  kf/k;:—  \  , 

;,'  =  <?  '  •  Sr  =  '^ 

This  leads  to  the  secular  equation 

ipr  ~pt)'%\nkfL  iink^L  —  Zarpticosik^  —k-j-'^L  —  1  ]=0 


(A4) 


with  the  constraint  on  the  frequency,  namely, 


I 


» 


» 


» 


» 


(A5' 


£0*  ■“  ct>5  f  ^jr  *  ^ 


The  amplitudes  are  as  follows: 

B - - - ^ 

^  gri/u-gL--gr  -pi(fL-SL^  0 

(. - ^ - -  ^ 

A  q 

^  _  orilfr-fi  ^-pnf.^  -gi_ )  0 

D - - - -  ^ 

9 


(A6) 


(A7) 


(A8) 


£^-qr{fj.-gT)-rpt{fT^gT--gL'>  ■  <A10) 

Sometimes  it  is  useful  to  express  the  secular  equation 
(AS)  as  follows: 

iqr  Sind i^cosd T-^pt  cos0|_sine;-) 

Xiqr  cosdi_sin6-p^pt  sind^cQsdj-)=0  .  (All) 

where 


(A9) 


e^=ki_L/2  and  dr  =  k-rL/2  . 
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ABSTRACT 

Hvbrid  optical  modes,  consisting  of  a  unique  linear  combination  of  LO,  TO  and  IP  (interface  polaritonl  modes 
satisfying  both  elastic  and  electromagnetic  boundary  conditions,  are  described  in  a  continuum  approximation. 
Mode  patterns  and  mode  dispersion  are  obtained  which  agree  with  microscopic  calculations.  Inter-  and 
intrasubband  electron  scattering  rates  are  obtained  which  also  agree  with  those  derived  on  the  basis  of 
microscopic  models. 

1.  iNTROOUCnON 

Continuum  theories^--  of  confined  LO  modes  hitherto  fail  to  satisfy  both  elastic  and  electromagnetic  boundary 
conditions  and  do  not,  in  general,  describe  properly  the  mode  patterns  and  mode  dispersion  obtained  in 
microscopic  calculadons.-^’’*-^  In  strucrures  where  wells  and  barriers  are  many  unit  cells  thick,  conbnuum  theory 
has  an  important  role  to  play  in  iaciiitar.ng  the  calculation  of  electron  scattering  rates  and  of  phonon  decay 
rates,  and,  not  least,  in  providing  a  conceptual  framework.  Here  we  describe  how  hybrid  opncal  modes 
consisting  of  unique  linear  com.binations  of  -C,  TO  and  IP  (interface  poiariton)  modes  solve  the  boundary- 
condition  problem,  and  provide  a  descr.ption  which  agrees  remarkably  well  with  the  results  of  microscopic 
theory. 

2.  LO,  TO  AND  IP  HYBRIDS 

.Assuming  a  quadratic  dispersion  for  LO  and  TO  modes  and  electromagnetic  dispersion  for  IP  modes  in  an 
elastically  isotropic  medium,  we  take  the  relative  ionic  displacement  u  to  be  a  linear  combination  of  LO,  TO 
and  IP  components  which  share  the  same  frequency  and  the  same  in-plane  wavevector,  which  makes  the 
combination  unique: 

u  =  uIjO  +  irpo “IP 

with  an  associated  electric  field  E  given  by 

E  =  -Po(uLO  +  sujp),  2. 

where  Po  =  e*/eoVo  is  the  space-charge  density,  e*  is  the  effective  ionic  charge,  Eq  free-space 

permittivity  and  Vg  is  the  volume  of  a  primitive  unit  c^ll.  have,jassumed  that  the  TO  mode  has 
vanishingly  small  electric  field.  The  parameter  s  =  (u-  -  <“-fo)/(WL“o  “  *^TO^  ^  usual  fraction  of  the  LO 
space-charge  density  which  is  associated  with  an  IP  mode  of  frequency  o),  and  colq,  toyo  zone- 

centre  frequencies.  The  displacement  components  obey  longitudinal  and  transverse  condidons:  d  x  ulO  = 
A.uTO  =  0>d.uip  =  0. 

For  elastically  infinitely  rigid  interfaces  (roughly  corresponding  to  the  AlAs/GaAs  interface)  the  boundary 
conditions  are  u  =  0,  Ex  continuous,  D^  continuous  at  z  =  ia/2,  where  x  is  along  the  plane,  z  is  perpendicular  to  it, 
and  a  is  the  well-width.  The  hybrids  whose  frequencies  lie  near  to  u>LO  (ss  distinct  from  eojO^ 
which  interact  most  strongly  with  electrons.  Mode  patterns  for  these  turn  out  to  be  given  by 


2  ■  $Pi^  Voi  1675  Quirxum  Well  tnd  Superimtee  Physics  IV  (1992) 


0^1  $4-03 36-0/92/S-*  00 


ux  « ikx  ~  -  arsu'i' tr*  -  *p 


(a) 


3. 


Odd 


Uz  =  kL  ~  [cos  kL2  -  aT<k^/kLkT)cosh  kTZ  -  ap(kx/kL)cosh  kxzl 


(b) 


or 

ux  =  kx  ~  [cos  kLZ  -  bT  cosh  Icp  -  bp  cosh  kxzl  (a) 

Even 

uz  »  ikL  ~  [sin  kLz  -t-  bT(k“/kLkT)sinh  k-fz  +  bp(kx/kL)sinh  kxz] 

(b) 


where 


<a^  =  (i>Lo  -  '‘x  ^  “  ®TO  ~  ~  (c^/eiuHioKk^  -  k^)  • 

Here  VL,  vj  are  velocities  describing  the  dispersion,  c  is  the  velocity  of  light,  e((o)  »  CmCo^  -  - 

0)-^),  Ho  is  ihe  permeability,  and  kp  is  the  evanescent  wave-vector  of  the  IP  mode,  which  in  etp.  3  an  4  has 
bew  taken  to  be  equal  to  kx,  assuming  the  non-retarded  limit.  The  parameters  ax,  ap,  bx,  bp  are  given  by 

(1  -p|tanh(kxa/2)lsin(kLa/2)  sin(kLa/2)  _ I _ 

cosh<kxa/2)  '  ~  cosh(kxa/2)  '  ~  s(tanh(kxa/2)+r) 


5. 


1 -p2Coth(kxa/2)cos(kLa/2)  co5(kLa/2>  1 

sinh(kxa/2)  '  *^"'^sinh(kxa/2)  '  P2*s(coth(kxa/2)+r) 

In  these  equations  r  *  e(a))/eB,  where  cb  is  the  permittivity  in  the  barrier,  assumed  for  simplicity  to  be 
frequency  independent,  and  it  is  further  assumed  that  kx  is  a  large  evanescent  wavevector  such  that 
tanh(kxa/2)  «  1.  The  real  waveveaor  of  the  LO  component,  kL,  is  determined  from  the  dispersion  relations, 
which  are 

Odd  modes:  cot(kLa/2)  »  (kx/kLXpi  +  (kx/kxKl  -  P'jtanh(kxa/2)J  (a) 

6. 

Bvennrxxies:  Un(kLa/2)  »  -  (kx/kL)(p2  +  (kx/kxXl  -  P2COth(kxa/2))l  (b) 

In  the  barriers  u  s  0  an  there  are  only  evanescent  electric  fields. 

Mode  patterns  for  small  in-pIane  wavevector,  kx,  are  depicted  in  Rg.  1  for  the  four  highest  frequency  modes. 
These  are  closely  like  patterns  obtained  in  microscopic  theory.  Fig.  2  depicts  the  dispersion  for  a  lOOA  well. 
Odd  modes  are  affected  by  the  IP  dispersion,  even  modes  are  not.  In  particular,  the  LOl  mode  rapidly  evolves 
into  a  mode  near  the  original  L03  frequency.  This  ’disappearance*  of  the  LOl  mode  is  characteristic  of  the 
effect  of  hybridization  and  was  first  pointed  out  by  Huang  and  2hu^  in  their  microscopic  theory.  It  is 
interesting  that  this  property  is  also  obtained  in  a  continuum  theory. 
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Fig-  2  Dispersion  plots  of  ItL  versus  kx,  solutions  of  cq.  6,  for  a  =  100 A. 
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3. 


THE  ELECTRON-HYBRIDON  INTERACTION 


These  hybrids  have  mainly  mechanical  energy,  as  have  phonons,  but  the  interaction  with  electrons  is  partly 
electric  and  partly  magnetic,  the  LO  component  contnbuting  a  scalar  potential,  o,  the  IP  component  a  vector 
potential  A.  Thus 

0  =  Po  sin  k^z  (a) 


Odd  7. 

Ax  *  ^  ap  sin  kxZ  (b) 

Az  =  apcosh  kxZ  (c) 

0  = -po  Ae*^^xx~“t)  coskLZ  (a) 

Even  8. 

Ax  =  bp  cosh  kxZ  (b) 

Az  =  -  .^giCkxXHDt)  j (c) 


It  is  interesting  that  in  general  neither  0  nor  A  is  continuous  at  the  interface,  but  what  matters  is  the  interaction 
with  an  electron,  H  =  -e0  +  (e/m)  A  .  p,and  this  turns  out  to  be  continuous  provided  that  the  electron  travels 
coherently  with  the  ware,  i.e.  with  a  group  velocity  equal  to  ta/kx-  This  condition  generalizes  the  usual 
electrostatic  argument  for  o  to  be  continuous  to  the  time-dependent  situation.  It  should  be  noted,  in  this  context, 
that  in  microscopic  theories  hitherto,  electric  fields  have  been  always  associated  with  a  scalar  potential,  an 
association  which  ignores  the  time-dependent  properties  of  the  system. 

The  hybrid  field  can  be  quantized  to  give  "hybridons"  and  the  scattering  rates  of  intrasubband  processes  can  be 
computed  using  the  Fermi  Golden  Rule.  Fig.  3  shows  the  intrasubband  rate  for  an  electron  in  the  lowest  subband 
of  an  infinitely  deep  well  in  GaAs  when  the  electron  energy  =  ha).  The  interaction  is  solely  with  GaAs 
hybridons.  Comparison  is  made  with  the  predictions  of  the  hydrodynamic  model,  the  dielectric  continuum 
model,  the  bulk-phonon  model  an  the  Huang-Zhu  model.^  The  agreement  with  the  latter  model,  which  is 
based  on  a  scheme  of  assigning  envelope  functions  to  the  results  of  microscopic  theory,  at  small  well-widths,  is 
striking.  Fig.  4  shows  the  rates  for  the  intersufaband  process  of  an  electron  an  the  bottom  of  subband  2  scattering 
into  subband  1  via  the  errussion  of  a  (converted)  LOl  hybridon.  (Addition  of  contributions  from  LOS  etc.  would 
close  the  gap  between  our  result  and  that  using  Huang  and  Zhu  modes.) 

4.  CONCLUSIONS 

The  continuum  model  of  hybridization  successfully  predicts  confined-mode  patterns  and  electron  scattering 
rates.  It  shows  that  hybridization  causes  the  rates  to  lessen  with  decreasing  well-width.  Previous  models 
showed  this  reduction  to  be  true  only  for  LO  mode  interaction  but  not  for  the  IP  interaction.  Hybridization 
forces  the  IP  interaction  into  the  LO  mold.  It  is  only  in  a  material  like  AlAs  where  LO  dispersion  is  too  weak  to 
allow  hybridization  over  all  the  frequency  range  between  (olq  ‘^TO  modes  can  retain  their 

characteristic  behaviour. 
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Fig.  3  Intrasubband  scattering  rate  for  an  electron  with  energy  hcs  interacting  with  L02,  in  an  infinitely  deep 
CaAs  well.  (HD  =  hydrod}mamic  model,  DC  =  dielectric  continuum  model,  HZ  =  Huang-Zhu  model). 


Fig.  4  Intersubband  scattering  rate  for  an  electron  at  the  bottom  of  subband  2  scattering  to  subband  1  emitting  a 

"converter"  LOl  hybridon. 
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A  simple  analytical  model  is  obtained  to  descrioe  the  effect  of  carrier  heating  on  the  frequency 
response  of  a  quantum  well  laser.  The  principal  factors  are  taken  to  be  injection  hearing, 
recombination  heanng,  and  hot  phonons.  The  model  is  applied  to  the  GaAs/GalnAs  stramed 
layer  system  and  is  shown  to  qualitatively  account  for  many  of  the  nonideal  features  observed. 
The  nonlinear  effects  cannot  be  described  satisfactorily  by  a  single  phenomenological  “gain 
suppression”  factor.  However,  at  low  drives  the  conventional  gain  suppression  factor  can  be 
e.tpressed  in  terms  of  the  phonon  lifetime  and  the  temperature- relaxation  time.  The  response  is 
mediated  by  several  time  constants  which,  in  our  e.xample.  combine  to  give  an  effective  time 
constant  of  about  10  ps.  The  modulation  frequency  response  becomes  seriously  impaired  when 
the  differential  gain  is  lowered  by  a  factor  of  2  and  the  time  constants  describing  scattering  and 
phonon  lifetime  are  increased  by  a  factor  of  2- 


I.  INTRODUCTION 

The  effect  of  carrier  heating  on  the  performance  of 
injection  lasers  has  been  considered  by  a  number  of  au¬ 
thors.'"*  The  source  of  heating  has  oeen  considered  to  be 
the  total  power  dissipated.’  .Auger  processes.'  recombina¬ 
tion.'*  and  free-camer  absorption."’  *  and  recent  work  has 
related  earner  heating  to  so-called  nonlinear  gain  suppres¬ 
sion  directly^  and  in  combination  with  hole  burning.*  That 
earner  heating  occurs  in  lasers  is  well  established^  and 
femtosecond  studies  suggest  that  it  alone  can  ex,  'n  gain 
suppression.*  Since  much  of  the  foregoing  work  ha*  ...tn 
on  bulk  lasers,  it  is  of  interest  to  study  carrier  heating  in 
quantum  wells  and  in  particular  to  elucidate  its  effect  on 
the  modulation  frequency  response,  and  that  is  the  topic  of 
this  anicle. 

Our  interest  will  be  focused  on  high-speed  (and  there¬ 
fore  short )  quantum  well  strained  layer  lasers  in  which  the 
effective  masses  of  the  electrons  and  holes  in  their  lowest 
subbands  are  similar  in  magnitude.^  The  main  source  of 
external  heating  in  the  quantum  well  case  is  that  provided 
by  energetic  carriers  injected  at  the  top  of  the  well.  These 
relax  via  intercamer  collisions  and  via  optical  phonon 
emission,  thus  heating  both  comer  and  phonon  gases.  Of 
the  heating  meciianisms  invoked  in  previous  work  we  re¬ 
tain  only  recombination  heating  and  regard  the  contribu¬ 
tion  from  Auger  processes  and  free-carrier  absorption  as 
being  negligible.  Comprehensive  work  on  energy  relaxation 
in  bulk  and  quantum  well  material,  recently  reviewed.* 
shows  clear  evidence  that  at  high  carrier  densities  the  role 
of  hot  phonons  cannot  be  neglected.  We  have  therefore 
assumed  that  the  most  important  factors  to  be  taken  into 
account  are  injection  heating,  recombination  hearing,  and 
hot  phonons,  and  we  explore  how  they  influence  the  per¬ 
formance  of  a  laser.  As  far  as  we  are  aware  such  an  analysis 
has  not  been  carried  out  previously. 

The  system  to  be  desenbed  is  a  complex  one  consisting 
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of  several  populations  (electrons,  holes,  photons,  phonons) 
dynamically  coupled  together.  Detailed  modeling  can  only 
be  carried  out  numerically.  Our  purpose  in  this  article  is 
not  to  attempt  a  numencal  analysis  but  instead  it  is  to 
construct  an  analytic  model  which  embodies  the  funda¬ 
mental  physical  processes  involved  and  which  illustrates 
the  factors  of  importance.  Such  a  model  is  a  preliminary 
(in  our  view  a  necessary  preiimmary)  to  a  computer  sim¬ 
ulation. 

This  article  is  organized  as  follows:  In  Sec.  II.  we  dis¬ 
cuss  energy  relaxation  in  a  degenerate  gas.  A  number  of 
investigations  for  GaAs  show  that  at  densities  of  order  I0‘* 
cm'*  a  time  constant  characterizing  electron-electron  scat- 
tenng  and  plasmon  emission  is  of  order  60  fs  or  less.* 
Given  the  assumed  near  equality  of  electron  and  hole 
masses,  that  allows  us  to  adopt  the  simplification  of  a  com¬ 
mon  carrier  temperature  established  virtually  instanta¬ 
neously.  We  simplify  the  electron  (hole) -phonon  interac¬ 
tion  in  a  quantum  well  by  ignonng  hybridization  of  the 
LO.  TO.  and  IP  (interface  polariton)  modes  and  take  the 
rate  to  be  proportional  to  the  density  of  subband  states. 
This  allows  us  to  write  down  simple  coupled  equations  for 
the  carrier  temperature  and  the  phonon  number,  taking  the 
degenerate  limit  and  we  bnefly  explore  the  linear  response 
of  the  system.  This  leads  to  the  identification  of  a  new  time 
constant  which  characterizes  the  temperature  relaxation  in 
a  degenerate  gas.  In  Sec.  III.  the  equations  are  combined 
with  the  standard  ones  for  photon  and  carrier  densities, 
and  the  injected  energy  is  related  to  the  current  density. 
The  linear  response  of  this  coupled  system  of  carriers,  pho¬ 
nons.  and  photons  is  analyzed  and  the  resultant  response 
function  is  compared  with  the  ideal  (no  carrier  heating) 
case.  A  distinction  is  discovered  between  the  frequency  of 
maximum  response  and  the  frequency  of  relaxation  oscil¬ 
lations.  and  the  inadequacy  of  describing  the  nonlinear  ef¬ 
fects  which  are  predicted  by  a  single  gain  suppression  fac¬ 
tor  IS  pointed  out.  In  Sec.  IV,  the  theory  is  applied  to  the 
GaAs/Gao.7lnojAs  strained  layer  laser  and  the  response 
function  is  found  for  a  number  of  situations  in  all  of  which 
the  current  density  and  photon  flux  is  taken  to  be  high  in 
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order  to  illustrate  the  noniinearines.  Qualitative  features  of 
nonideal  behavior  observed  in  practice  are  well  reproduced 
by  the  theory.  However,  because  of  the  simplistic  nature  of 
the  model  no  attempt  is  made  to  obtain  a  quantitative  fit. 
The  article  concludes  in  Sec.  V  with  a  brief  discussion. 

Our  aim  throughout  is  to  present  a  simple  analytic 
model  in  order  to  illustrate  the  processes  involved.  We 
conclude  that  carrier  heating  itself  does  not  seriously  im¬ 
pair  the  frequency  bandwidth,  but  that  it  complicates  the 
response  and  forces  a  reappraisal  of  how  nonlinear  effects 
should  be  interpreted. 

_ _ 1 

— 

It 


IL  ENERGY  RELAXATION 

The  rate  of  increase  in  energy  density,  U,  can  be  ex¬ 
pressed  in  terms  of  the  rate  of  change  of  occupation  prob¬ 
ability,  /{£),  of  a  state  with  energy  E  as  follows: 

dU  r  df{E) 

where  S(E)  is  the  density  of  states  function.  In  cases 
where  the  occupation  probability  changes  solely  by  the 
emission  and  absorption  of  optical  phonons  of  energy  ffu 
Eq.  { 1 )  can  be  written  in  the  following  form: 


J EN{E)dE^  ^(.V-(£,£’){n(e))/(£')[l-/(£)]-[n(a))-(-l]/(£)[l-/(£')]} 


x5(£’-£-hi?a»)Af(£')d£')->-(3/-(£,£”){(n(iu) 
X  5(£"  -  £-ffai  )^(£”  )d£"). 


where  M^{ES')  is  the  squared  matrix  element  for  the 
transition,  n(<u)  is  the  phonon  occupation  number,  and  the 
second  integral  is  over  allowed  final  states.  This  integral  is 
particularly  simple  when  the  matrix  element  is  indepen¬ 
dent  of  phonon  wave  vector,  for  then: 

X  [  1  -/(£)  ]  -  [n(6))  + 1  ]/(£)  ( 1  -/(£-&))  ]  1 

-/t(a>)/(£)[l-/(£^fiai)Jlj.  (3) 

where  T{E±1mJE)  is  the  scattering  time  constant  for  the 
process.  The  two  types  of  electron-phonon  interaction  of 
importance  are  the  polar  and  the  nonpolar  (in  a  polar 
material).  To  a  good  approximation  the  nonpolar 
interaction — relevant  for  holes  and  electrons  in  L — is  me¬ 
diated  by  a  matrix  element  which  is  independent  of  wave 
vector,  but  the  polar  interaction  is  not  so  simple  in  general. 
However,  in  a  two-dimensional  (2D)  system  the  depen¬ 
dence  of  the  polar  interaction  on  wave  vector  is  weakened 
by  confinement,  and  where  confinement  is  strong  Eq.  (3)  is 
approximately  valid  even  for  the  polar  case.  A  further  sim¬ 
plification  which  can  be  made  for  a  2D  system  in  which 
only  the  lowest  energy  subband  is  occupied  is  that  the 
scattering  time,  which  is  inversely  proportional  to  the  sin¬ 
gle  spin  density  of  sutes,  can  be  taken  to  be  energy- 
independent  if  nonparabolicity  is  ignored.  In  what  follows, 
both  of  these  simplifying  features  will  be  adopted. 

With  t(£)  =T  and  Vf  F)  =  iV  plus  assumption  that 
the  distribution  function  is  of  Fermi-Dirac  form  charac¬ 
terized  by  an  electron  temperature  T„  the  integrals  in  Eq. 


> \)f(E")[  1  -/(£)  ]  -n(£^)/(£)  [  1  -/(£”)  1 1 


(2) 


r 


(3)  can  be  evaluated  analytically  in  the  degenerate  limit 
VIZ.:  exp(£fAar,)>l,  where  £/-is  the  Fermi  energy  mea¬ 
sured  from  the  subband  edges.  Since  degeneracy,  or  near 
degeneracy,  is  a  prerequisite  for  laser  action  we  have  taken 
the  degenerate  limit  in  order  to  obtain  simple  analytical 
results  which  illustrate  hot-carrier  effects.  Thus,  for  a  2D 
degenerate  electron  gas  the  rate  of  increase  in  energy  den¬ 
sity  reduces  to 

dU  fia 

-j-= - {rt,(<y)  — rt(a/)  JiVAy,  (4) 

dt  T 


where  n,(w)  is  the  Bose-Einstein  occupation  number 
when  the  phonon  temperature  equals  the  electron  tem¬ 
perature.  The  rate  is  zero  when  £,=  T^^Tl,  where  is 
the  lattice  temperature,  as  it  should  at  thermodynamic 
equilibrium.  The  factor  N'fko  is  the  effective  density  of  elec¬ 
trons  which  can  interact  with  the  phonons,  which  repre¬ 
sents  a  fraction  fko/Ef  of  the  totaL  In  a  nondegenerate 
system  the  corresponding  equation  would  have  the  total 
carrier  density  replacing  Nfiiit.  Thus,  degeneracy  reduces 
the  energy  relaxation  rate  by  the  fraction  fia/Ef. 

In  a  nondegenerate  gas.  energy  relaxation  translates 
linearly  to  temperature  relaxation,  but  this  is  not  the  case 
when  the  gas  is  degenerate.  In  2D  the  mean  energy  (£)  is 
given  for  exp{Ef/kgT,)>l  by 


17^  (kaTA 


(5) 


[Note  that  in  the  derivations  of  Eqs.  (4)  and  (S),  identical 
methods  were  used  in  evaluating  the  integrals,  and  these 
equations  are  therefore  mutually  consistent.]  With  U 
=n{E)  and  Ea=^n/N,  where  n  is  the  electron  density,  Eq. 
(4)  becomes 


» 


» 


» 


I 


» 


» 


> 


i 

i 


► 


2580  J.  Appl.  Phy».,  Vol.  72.  No.  7.  1  October  1992 


L  F.  Lesur  and  B.  K.  Ridley  2580 


^4-y  krjT,  S  [«,{«)  -/i(ffl)  lJV*a. 


(6) 

The  temperature  relaxation  rate  then  depends  upon  the 
earner  density  relaxanon  rate  as  well  as  the  energy  relax¬ 
ation  rate.  Indeed,  if  the  latter  is  zero  the  temperature  rises 
as  the  earner  density  falls,  and  this  effect  was  termed  re¬ 
combination  heatmg  by  Bimberg  and  Mycielski.’'‘°  On  the 
other  hand,  if  the  earner  density  remains  constant  the  tem¬ 
perature  rela.xauon  rate  is  given  by 


dT,  ftbi 


;cs— [«,(&») 


where  the  new  time  constant  is  given  by 
rkaT.r 


_* _ 

'  r" 


l-ftm 


(7) 


(S) 


This  time  constant  would  enter  a  descripdon  of  the  tem¬ 
perature  fall  off  after,  say,  the  opdcai  injection  of  hot  elec¬ 
trons,  as  observed  by  time-resoived  photoluminescence, 
provided  no  hot  phonon  effects  were  present. 

To  calculate  the  effect  of  hot  phonons  we  adopt  the 
simplifymg  assumption  that  a  phonon  temperature  can  be 
defined  for  the  modes  interacting  most  strongly  with  the 
electrons.  We  then  assume  that  the  rate  of  change  of  oc¬ 
cupation  number  can  be  expressed  as  follows: 


dn((ii)  dU  V  [n(o))—n,(u)] 
dt 


(9) 


dt  NJm  Tf 

where  .V„  is  the  number  of  phonon  modes.  V  is  the  volume 
of  the  quantum  well.  is  the  occupaQon  number  at 
thermodynamic  equilibrium,  and  is  the  phonon  lifetime. 
We  assume  all  phonons  involved  have  the  same  frequency 
and  decay  at  the  same  rate,  which  is  very  nearly  true.  To 
establish  the  number  of  modes  we  observe  that  momentum 
and  energy  conservation  restricts  the  phonon  wave  vector, 
q,  to  lie  between  and  q,^  where  for  emission: 

for  absorption:  ( 10) 

where  E,  k  are  the  energy  and  in-plane  wave  vector  of  the 
electron.  In  our  model  of  degeneracy  we  can  for  simplicity 
talte  0<q<  2k f,  where  k^  is  the  wave  vector  on  the  Fermi 
surface.  The  number  of  modes  is  then  given  by 


AT.. 


y=iny. 


(U) 


dn(cii)  Af  [/i(ai)— e,(a*)] 

— [«.(a.)-e(<u)]  - - - - 


(12) 


The  phonon  number  in  the  steady  state  is  then  given  by 
(.Vfta/2n)  (Tj, /r) 


n(u>)=- 


(13) 


l-!-(iVAu/2n)(r^/T) 

and  n{cii)  =n,(ai)  only  when  /i,(<a)  =n,(Q)).  In  extremely 
degenerate  systems  such  that  (iVAu/lnKr^/rl^l,  hot- 
phonon  effects  are  negligible.  But  typicaUy,  r^/T>I,  and 
so  hot-phonon  effects  can  occur  even  when,  as  we  have 
been  assuming,  (fi(a/£f)<l. 

We  have  implicitly  assumed  that  the  normalization 
volumes  of  the  carriers  and  phonons  are  idendcal,  but  in 
general  this  will  not  be  so.  For  e.xample,  vibradons  may  be 
e.xcited  in  the  barriers  in  addidon  to  being  excited  in  the 
wells.  Following  an  analogy  with  the  situadon  for  photons 
in  a  semiconductor  laser,  we  can  introduce  a  phonon  con¬ 
finement  factor,  which  is  generally  less  than  unity,  and 
interpret  the  scattering  rate,  r“',  as 


(H) 


where  Tq  is  the  scattering  time  for  r^  =  1- 

The  linear  response  of  the  coupled  system  can  be  de- 
senbed  by  introducing  a  power  source  P  on  the  right-hand 
side  of  Eq.  (6)  with  a  small  moduladon  at  an  angular 
frequency  H.  We  obtain 

n^rH  1  +a-<nrp  {SPr/Nfiu) 

n; - n V /ft [  ( l/rr)  -r  ( l/r,)  ( 1  +a )  ]  ’  ' 

where 

a=(iVftu/2n)(r,/r). 


and 
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The  dc  response  is 

•I] 


(16) 


*V,Ar* 


(17) 


(18) 


where  is  the  well  width  and  K  is  the  well  volume.  [Note 
that  Eq.  (11)  overestimates  the  number  of  modes  if  degen¬ 
eracy  is  weak  and  hence  leads  to  an  underestimation  of 
hot-phonon  effects.]  Thus, 
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kg5T=TT(\-ha)(8Pr/N^). 

Normalizing  to  the  dc  response  leads  finally  to  the  re¬ 
sponse  funedon  £(ft)  given  by 

ftjd-f/ftrf) 

/?(ft)=7^T - ^ 


(19) 


'ft;-ft^-Hft((l/rr)-(-(i/T,)l’ 

where  tJ  =  r/(  1  a).  In  Eq.  ( 17).  /i,(«)  is  given  by 
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n,(cij)=n,(u;) 


?r(l-a) 
iV(Ai»)'  ■ 


(20) 


where  P  is  the  steady-state  power  density  input. 

The  time  constants  characterizing  the  response  are  tj-, 

Ir,  and  r^.  Typically,  at  room  temperature,  each  of  the 
three  time  constants  is  of  order  10  X  *  i.e.,  of  the  order  of 
1  ps.  In  other  words,  degeneracy  and  hot-phonon  effects 
reduce  the  energy  reiaxation  rate  by  an  order  of  magni¬ 
tude.  This  reduction  is  due  partly  to  the  reduced  number  of 
electrons  which  can  participate  and  partly  due  to  a  hot- 
-  phonon  effect.  It  should  be  noted,  however,  that  the  hot- 
phonon  effect  weakens  with  increasing  degeneracy  as  more 
and  more  phonon  modes  becomes  available  and  the  occu¬ 
pation  of  each  becomes  progressively  less  (a— 0).  In 
choosing  the  degeneracy  limit  we  underestimate  the  hot- 
phonon  effects  somewhat. 


ill.  LASER  DYNAMICS 

In  a  semiconductor  quantum  well  laser  earner  heating 
occurs  through  the  act  of  injecting  electrons  and  holes 
across  band-edge  discontinuities.  If  the  latter  are  and 
for  the  discontinuities  in  the  conduction  band  and 
valence  band  then  an  amount  of  energy  A£j=  A.Ej  — A£,is 
injected  for  every  electron-hole  pair.  In  the  presence  of  a 
substantial  photon  ffux  there  will  also  be  heating  associated 
with  free-camer  absorption,  but  for  simplicity  we  will  as¬ 
sume  that  this  is  small  compared  with  that  unavoidably 
associated  with  the  injection  process.  The  questions  to  be 
answered  are:  how  is  that  injected  energy  dissipated?  and 
how  does  heating  affect  performance  of  the  laser? 

To  answer  these  questions  we  will  assume  for  simplic¬ 
ity  that  complete  carrier  thermalization  occurs  virtually 
instantaneously  so  that  a  single  temperature,  instanu- 
neously  characterizes  the  quasi-Fermi  distributions  of  elec¬ 
trons  and  holes.  For  300  meV  electrons  in  bulk  GaAs  in  a 
concentration  of  order  lO”  cm“"  the  mean  free  path  asso¬ 
ciated  with  plasmon  emission  and  single  particle  processes 
is  about  SOO  A  according  to  the  calculations  of  Rorison 
and  Herbert,"  corresponding  to  a  time  constant  of  about 
40  &.  Monte  Carlo  estimates  of  thermalization  dme  in  a 
quantum  well  of  GaAs  gives  200  fs."  These  times  are  cer¬ 
tainly  short  enough  to  warrant  our  assumption  as  regards 
electrons.  This  will  also  be  true  for  holes  provided  the  hole 
mass  is  similar  to  that  of  the  electron  which  it  is  in  strained 
layer  lasers  of  GaAs/lnGaAs.  In  bulk  material,  the  dispar¬ 
ity  between  electron  and  hole  masses  means  that  equaliza¬ 
tion  of  electron  and  hole  temperatures  is  slow,  full  ther¬ 
malization  taking  times  of  order  10  ps."''^  But  in  the  case 
of  strained  layer  lasers  thermalization  times  will  be  of  or¬ 
der  100  fs.  Thus,  our  assumption  will  be  reasonably  valid 
for  (working)  strained  layer  lasers,  and  may  be  a  reason¬ 
able  approximation  for  quantum  well  lasers  in  general,  but 
it  is  certainly  invalid  for  bulk  lasers. 

Our  assumption  of  instantaneous  thermalization  al¬ 
lows  us  to  make  a  funher  simplification,  over  and  above 
that  concerning  a  common  temperature,  regarding  veil 
capture.  In  detail  this  process  is  difficult  to  describe  accu¬ 
rately,  but  if  it  is  assumed  that  carrier-carrier  interactions 


are  strong  enough  to  establish  a  common  electron-hole 
temperature  there  is  no  need  to  consider  well  capture 
explicitly — it  occurs  as  rapidly  as  the  thermalization  pro¬ 
cess. 

Energy  relaxation  of  the  strongly  coupled  electron-hole 
degenerate  gases  is  via  the:  emission  of  optical  phonons  (or, 
more  accurately,  hybridons,  the  quanu  of  the  hybrid  LO, 
TO,  and  IP  modes  which  occur  in  layered  systems,  but  we 
will  continue  to  use  the  more  familiar  term),  through  the 
Frohlich  interaction  for  both  types  of  earner  and  through 
the  deformation  potential  interacuon  for  holes.  In  general, 
different  phonon  modes  are  involved  because  of  the  differ¬ 
ence  in  mass  of  the  electrons  and  holes  and  because  the 
deformation  potential  interaction  involves  TO  as  well  as 
LO  modes.  We  should  therefore  consider  four  separate 
phonon  populations  viz.  LO  (electrons),  LO  (holes),  TO 
(holes),  TO  (holes),  and  more  if  IP  modes  are  involved. 
To  avoid  this  complication  we  assume  that  we  can  work 
with  a  single  set  of  phonon  modes  whose  occupation  is  an 
average  of  the  true  LO  and  TO  occupation  numben  and 
whose  associated  emission  and  absorption  rates  are  char¬ 
acterized  by  a  single  time  constant  r.  Such  an  assumption 
is  relatively  crude,  but  it  does  not  seriously  falsify  the  basic 
physics.  We  come  back  to  this  point  later. 

The  equation  describing  energy  relaxation  is  now 

dU  fao 

(21) 

at  r 


where  U  is  the  total  energy  density  of  the  electron-hole  gas 
in  a  single  well,  cu  is  an  average  phonon  frequency,  N/t 
=  {N/r)f-^{N/r)^  and  Jg  is  the  energy  current  density. 
The  latter  may  be  related  to  the  current  density,  J,  as 
follows.  The  panicle  currents  J,  and  7*  are  related  to  J  via 
and  thus  7,-r  A£,.7*),  where  A£j. 

and  are  the  relevant  band-edge  discontinuities  in  the 
wells  minus  the  subband  energies.  At  steady  state  electrical 
neutrality  entails  that  J,—J„  and  hence 
where  t}  is  the  injection  efficiency.  If  d  is  the  cavity  width, 
the  divergence  of  the  energy  current  is  —  i7A£j7/2de.  and 
this  corresponds  to  minus  the  rate  of  increase  in  energy 
density  in  the  whole  cavity.  If  there  are  N„  wells,  the  rate 
of  increase  in  energy  density  in  one  well  is  therefore 

For  both  degenerate  electrons  and  holes,  we  have  from 
Eq.  (5). 


+ 


(22) 


In  the  case  of  electrical  neutrality  and  no  background  dop¬ 
ing  of  significance.  n=p,  where 


(23) 


where  IV,  is  the  reduced  density  of  states  per  unit  energy 
interval,  and  N/,  are  the  density  of  states  per  unit  cn- 
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ergy  interval  in  the  conduction  and  valence  subbands. 
Equadon  (21)  becomes 


n  dn  ir 


dt 


Hci) 

=  — —  [/i,(<u)  —n(ci)) \N-fia-rFJ, 


(24) 


where  we  have  introduced  an  effecdve  electric  field  F 
=  and  n  is  the  density  of  electrons  (or 

holes)  in  a  single  well. 

The  phonon  occupadon  is  to  be  described  by  the  fol¬ 
lowing  equadon 

dn{o})  -hbi  ,  (rt(<u)-/ra(«)l 

T— =r, —  n,{<u)-rt((u)  — - 

dt  ''  r  2n 

(25) 

in  which  is  a  factor  less  than  unity  introduced  to  allow 
for  the  diversity  of  phonon  modes  involved,  and  it  can 
serve  also  to  include  effects  due  to  differences  of  confine¬ 
ment.  It  will  henceforth  be  referred  to  as  the  phonon  mode 
factor. 

Addition  of  the  rate  equations  describing  the  photon 
fiua  density,  5,  and  the  electron  density  complete  the  set. 
These  are 


(26) 


dn  n  TfJ 


(27) 


where  u  is  the  group  velocity  of  the  photon,  T  is  the  photon 
confinement  factor,  g  is  the  gain,  r,  is  the  photon  lifetime, 
and  r,  is  the  spontaneous  recombinadon  lifetime.  We  have 
ignored  the  contribudon  of  the  spontaneous  emission  to 
the  photon  density  rate.  The  gain  is  given  by 


i=S^lAEe)-AEh)].  (28) 


TfJ  n 


(32) 


The  linear  response  of  the  coupled  canier-phonon- 
photon  system  can  be  obtained  in  the  usual  way,  and  not¬ 
ing  that  dg/dT,=  —C dg/dn.  The  response  fiincdon  turns 
out  to  be 


i2:(  I -i-idfl-fin^) 

^ -ricn-idO.^  -reCl*  ' 

where  the  relaxadon  angular  frequency  is  given  by 

2  n5(i+/) 


(33) 


(34) 


with 


,  rjcgFJr. 

n^=ug^/T„  (35a) 

2C(n(u)*n,(a))  l/i,(ia)  -I- 1] 


V.+rr(r,-l-^)  ( I  +r) 


Hereg,=dg/dn.  C=n/T^  r^ld-Zp-lCgS^^  and 
3  \  fiW  j  NNr  «,(&))  [/J,(li))-(- I]  ’ 


’=r-t- 


(35b) 


(36) 


(37) 


2T,C  ■ 

Tj-  is  just  a  generalizadon  of  ±e  temperature  relaxadon 
time  of  Eq.  ( 17)  and  r*  is  essentially  t(  I  -ha).  The  other 
parameters  in  Eq.  (33)  are  given  by 

-h  (  Tt/t)  [  r*  ( I  -i-  r-p )  -p(  rr^  /t,)  ] 

“  l+/-p(rrT-/rr,)  ’ 


b= 


fTt-pil+r-p) 


(39) 


2e^m*(p^/2m) 

?m«i“  r,  /-t  ’  (29) 

€r/nfrcti,r\,cL,^ 

where  m*  is  the  reduced  electron-hole  mass.  p„  is  the  mo¬ 
mentum  matrix  element  including  spin  selecdon.  eg  is  the 
permittivity  of  free  space,  m  is  the  free  electron  mass,  <a,  is 
the  photon  angular  frequency,  t\,  is  the  refracdve  index, 
and  c  is  the  velocity  of  light  in  vacuo.  For  simplicity,  line 
broadening  has  been  ignored.  For  thermalized  Fermi- 
Dirac  distribudons  and  n=>p. 


/(£,)  -/(£,)  =  I  -exp( -n/iV,  kgT,] 
-exp[-n/N^gT,] 


(30) 


and  thus  a  given  gain  determines  the  rado  n/T,^C.  At 
steady  state,  the  gain  is  given  in  Eq.  (26)  by 


I 


g-- 


vrr. 

and  from  Eq.  (27), 


(3l) 


where  p=rN,  ^jZT,C,  which,  like  r  and  /,  is  dimension¬ 
less  and 


I 


^  +  n^(r,+r,(l^/)-h^(l-hr)Jj,  (40) 
r*  +  r(5^n^r,r,(l+r)l||. 


fT 

r,+— 


(41) 

(42) 


Several  new  factors  enter  the  above  analysis.  First,  to 
calculate  the  response  funcdon  [Eq.  (33)],  one  notes  that 
the  cavity  losses  determine  the  gain  according  to  Eq.  (31). 
which  in  turn  determines  C—n/T,  for  all  currents  above 
threbold.  A  choice  of  T,  then  determines  n  and  all  the 
other  steady  state  quantities  including  the  current  density 
and  the  differential  gain.  Since  g,  decreases  with  increasing 
this  phenomena  is  a  source  of  nonlinear  gain.  The  oc¬ 
cupadon  numbers.  n[u)  and  n,{a>),  are  determined  under 
steady-state  condidons  from  Eqs.  (13)  and  (21),  respec- 


2583 


J.  Aopl.  Phys..  Voi.  72.  No.  7.  1  Oclooar  1992 


L  F  Lester  and  8.  K.  RicSey 


2583 


tivcly.  The  new  facton,  f.  r.  and  p  are  all  dimensionless. 
The  first,  /.  may  be  termed  the  injection  hot-phonon  fac¬ 
tor;  it  is  a  measure  of  the  energy  injected  into  the  phonon 
system  and  hence  of  the  increase  of  phonon  number  in  a 
.mode.  The  second  factor,  r.  is  the  recombination-heating 
coefficient;  it  is  a  measure  of  the  rise  in  temperature  asso¬ 
ciated  with  a  fall  in  earner  density.  The  third  factor,  p,  is 
the  direct-heating  coefficient,  which  is  a  measure  of  the  rise 
m  temperature  associated  with  the  injected  energy.  A  no¬ 
table  pattern  in  the  expressions  for  a  and  b  [Eqs.  (38)  and 
(39)]  is  the  term  (l-t-r— p).  This  reflects  the  variation  of 
photon  density  with  earner  density  and  temperature,  as 
these  vary  with  current  density.  Thus,  increasing  current 
density  results  in  increasing  carrier  density,  which  mani¬ 
fests  itself  as  the  unity  m  the  expression.  Increasing  tem¬ 
perature  reduces  photon  density  hence  the  negative  sign  of 
p.  and  the  positive  sign  of reflects  the  fact  that  increasing 
carrier  density  induces  cooling.  These  are  thus  competing 
mechanisms  at  work  here. 

The  above  response  is  significantly  more  complex  than 
the  standard  one  obtained  by  neglecting  carrier  heating 
altogether.  To  obtain  the  standard  result  it  is  necessary  to 
eliminate  both  sources  of  heating — recombination  and  in¬ 
jection.  The  latter  is  zero  if  F=Q,  where  /=0  and  p=Qi.  It 
is  interesting  that  this  does  not  produce  much  simplifica¬ 
tion.  This  points  to  the  main  source — recombination  heat¬ 
ing.  To  eliminate  this  means  making  both  r  and  api- 
proach  zero— letting  only  approach  zero  merely 
diminishes  hot-phonon  eflects.  and  letting  only  r  approach 
zero  merely  emphasizes  the  role  of  hot  phonons.  If  this  is 
done,  a-=b—d=e=Q,  and  the  standard  result  is  obtained 
viz. 


^(n)  = 


ns 


ftj— n*-ricn 


(43) 


with  c=r~'-rn5rp 

The  free  oscillations  of  the  system,  corresponding  to 
the  response  to  a  sudden  increase  in  current  density  when 
a  voltage  step  is  applied,  is  described  by  entailing  that  the 
denominator  in  Eq.  (33)  vanish.  Replacing  il  by  il  +  iy 
gives  two  equations  which  determine  the  relaxation  fre¬ 
quency  and  the  damping  rate,  viz.: 


—  1  3y<f)  — n*  — ycr-y^— dy^-i-ey*=0. 

(44a) 


4«7^-3dr-f2r(l-2en^)-c-dn*=0.  (44b) 

These  equations  are  most  easily  solved  numerically.  At  low 
drives  there  exists  a  fast  solution  which  is  approximately 
given  by 


I  ,  3i/ 

n-=-.^2r.  y=-.  (45) 


which  describes  the  temperature  relaxation  process.  There 
is  also  a  slow  solution,  which  is  very  roughly  given  by 


(2-t-3cd)n;-r 

2~cd 


c-dCll 
^“2(1-Hcd)  ■ 


(46) 


The  departure  from  ideal  behavior  is  often  described  in 
terms  of  a  nonlinear  gain  coeffiaent  such  that  g^gQ/(  i 
The  parameter  is  found  from  the  ratio  of  the 
damping  rate  y  to  the  square  of  the  resonance  frequency, 
f\.  This  rauo.  usually  named  K,  is  proporuonal  to  at 
low  drives.  From  Eq.  (46) 

,  (c  —  d£l\)[2  —  cd) 

K=2-r. - - - T-.  (47) 

(l^cd)[(2*3rf)n;-rl 

In  the  limit  of  small  current  densities 


A  = 


2irc 

i^-<r/2 


,,  ’■/D'O 

=  2ir [  Tj-i-T^-i — ^  ( 1  -i-r) 


(48) 


where  Tq  and  are  the  time  constants  when  the  carrier 
temperature  equals  that  of  the  lattice.  Small-signal  analysis 
gives  A=2rr^(r,-re„|/ug„„)  and  so.  comparing  with  Eq. 
(48),  we  obtain 


eni=«g,|r»-(-Tn,-^  (1-l-r)  j,  (49) 

where  g„  is  the  linear  differential  gain.  At  room  tempera¬ 
ture  the  total  time  constant  in  the  brackets  is  typically 
about  8  ps,  and  so  with  g,~10~'^  cm^  cm’, 

which  is  the  order  typically  observed.  However,  at  high 
dnves  the  nonlinear  effects  cannot  be  characteriz'd  by  a 
single  constant  parameter,  as  a  perusal  of  Eq.  (47)  indi¬ 
cates.  and  it  is  not  possible  to  distinguish  fast  solutions  and 
slow  solutions  as  clearly.  Generally  we  can  expect  the  fre¬ 
quency  observed  in  relaxation  oscillations  to  be  different 
from  that  of  the  peak  of  the  frequency  response. 


IV.  APPUCATION  TO  GaAs/Ga^.j^n^  STRAINED- 
LAYER  LASERS 

We  consider  the  system  in  which  x=0.3.  Nothing  def¬ 
inite  is  known  about  the  strength  of  electron  and  hole  in¬ 
teractions  with  optical  modes  in  this  system.  An  estimate 
can  be  made  along  the  following  lines.  Band-structure  cal¬ 
culations.  supported  by  the  results  of  magneto-optical  ex¬ 
periments.  suggest  that  the  electron  and  hole  masses  in  the 
lowest  subbands  are  about  equal  and  thus  the  strengths  of 
the  electron  and  hole  Frohlich  interactions  should  be  about 
equal  if  the  hole  overlap  integral  is  unity.  The  latter  should 
be  true  since  the  hole  subband  has  in-plane  symmetry  at 
energies  sufficiently  removed  from  the  next  higher  sub¬ 
band.  In  GalnAs  lattice-matched  to  InP  (x=0.53)  it  is 
known  that  the  Frohlich  interaction  is  much  stronger  for 
the  GaAs-like  mode  than  for  the  InAs-like  mode,’’  so  we 
ignore  the  effect  of  this  latter  mode.  For  x=0.3  the  GaAs- 
like  LO  mode  has  an  energy  of  about  34  meV  compared 
with  36  meV  in  GaAs,  and  thus  will  be  essentially  uncon¬ 
fined.  (GaAs  modes  with  energies  between  34  and  36  meV 
will  be  confined  in  the  bairieis  to  some  extent. )  Thus  the 
interaction  can  be  approximated  by  a  2D  carrier  interact¬ 
ing  with  bulk-like  LO  modes.  (Complications  associated 
with  interface  polariton  and  hybridization  will  be  ignored. ) 
From  the  results  of  Nash  er  aH”  we  deduce  a  coupling 
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parameter  a^o  of  0.055  (compared  with  0.068  for  GaAa) 
and  a  charactensdc  tMsic  dme  constant  Tq  of  17S  fs 
where 


i- =  2at.(>u.  (50) 

'o 

Comparison  with  bulk  rates  suggests  that  for  an  electron 
with  energy  2  Au  and  higher.  This  is  the  corresponding 
scattering  tune  and  this  is  also  true  for  quantum  wells.'* 
In  a  working  laser  with  carrier  densities  exceeding  10“ 
cm~\  it  is  necessary  to  take  screening  into  account.  The 
2D  plasma  frequency  is  approximately 


01 


? 


(ln)rqL.„ 

2m*e 


(51) 


where  2n  is  the  electron  plus  hole  densities.  if 

>  ct»Lo>  tmd  s  =  €5  if  <  coiQ,  e*.  are  the  high-frequency 
and  static  permittivities,  and  ?  is  the  wave  vector  of  the  2D 
plasma  wave.  Taking  ?  to  be  the  wave  vector  q„=(2m*a)/ 
I?)''*,  about  which  the  allowed  in-plane  wave  vectors  of 
LO  modes  distribute  themselves,  Q)^=<ttuo  for  £.^,=50  A 
when  n~2x  10'^  cm~\  Since  typical  densities  are  over 
10'*  cm~^  we  can  assume  and  as  a  crude  approx¬ 

imation  assume  that  the  screening  is  stadc,  which  implies 
reducing  the  Frdhlich  rate  by  a  factor  involving  and  the 
screening  vector,  q,=e~m*/-ff€/^  (normal  e,xpression  x2 
for  electrons  -i-  holes )  viz.; 

1  1  1 

where  the  factor  ^  takes  into  account  that  only  phonons 
traveling  in  the  plane  are  screened.  This  increases  the  time 
constant  to  368  fs. 

Being  short-range,  the  deformation  potential  interac¬ 
tion  is  not  as  afected  by  screening.  The  scanering  rate  in 
an  inhnitely  deep  well  is'’ 

1 

where  Dq  is  the  deformation  potendal  constant  and  p  is  the 
mass  density.  With  Dti=lX  10*  eV  cm~'  and  £„=50  A, 
To=2.(t\  ps.  Unfortunately  an  accurate  expression  for  Dfj 
is  not  available — the  value  chosen  here  is  that  for  the  L 
valleys  in  Ge — but  clearly  some  inaccuracy  here  will  be 
relatively  unimportant  in  view  of  the  relative  weakness  of 
this  interaction. 

In  Eq.  (21)  we  can  now  express  (he  total  scattering 
rate,  taking  into  account  our  assumption  that  as 

follows: 


1  _2 

To’ 


(54) 


whence  r=  170  fs.  In  these  processes  the  same  set  of  pho¬ 
non  modes  is  involved  in  both  the  electron  and  hole 
Frdhlich  scattering,  and  is  roughly  one  third  of  the  non¬ 
polar  scattering.  Thus  a  fracnon  [(2/r^)  -¥  (  1/3to)]t,  that 
is.  0.946,  of  scattering  processes  involve  the  same  set  of 
phonon  modes,  i.e.,  the  LO  modes,  and  so  we  can  regard 


Table  L  VUtenal  panmeten. 


ao7 

♦  p* 

Au, 

34  meV 

Sm.. 

6700  cm-' 

r 

170  fs 

Ir 

3.6 

30  A 

r 

0.086 

230  A 

1  ps 

■V, 

5 

r. 

2  ns 

r. 

0.1 3S 

the  ejects  of  hot  TO  phonons  as  being  negligible.  As  re¬ 
gards  coniinement  it  is  not  possible,  in  the  absence  of  any 
theoretical  study  of  optical  modes  in  this  system,  to  do 
anything  other  than  offer  a  crude  estimate.  Although  we 
e.xpec:  the  GaAs-like  LO  modes  to  be  relatively  uncon¬ 
fined.  this  does  not  mean  we  can  confidently  give  equal 
weight  to  well  and  barrier  and  accept  a  confinement  factor 
of  the  form  L,^[L^-rLg)\  where  L.^,  are  the  widths  of 
the  well  and  barrier,  because  the  frequency  mismatch  be¬ 
tween  well  and  barrier  will  result  in  there  being  different 
amounts  of  oscillator  strength  in  the  well  and  barrier.  Nev¬ 
ertheless.  for  consistency  of  the  model — the  scattering  has 
been  assumed  to  be  via  bulk-like  modes — we  will  ignore 
this  aspect  and  use  the  simple  form  mentioned  above.  This 
means  that  the  phonon  mode  factor  introduced  in  Eq.  (25) 
is 


V— ^ 

^  Vf  '  ^Td!  i-j.-rLg 


(55) 


The  phonon  lifetime  has  been  taken  to  be  the  same  as 
in  GaAs  at  room  temperature.'*  Other  parameters  are 
listed  in  Table  I.  We  consider  the  response  of  a  laser  con- 
tairung  five  quantum  wells.  1(X)  pm  long.  The  maximum 
gain  (cm*')  was  obtained  from  Eq.  (29)  with 
=  2.3  eV,  being  the  value  for  GaAs  multiplied  by  the  ratio 
of  the  photon  energy  to  the  GaAs  band  gap.  This  calcula¬ 
tion  assumed  an  injection  efficiency  of  100%  and  a  thresh¬ 
old  current  density  of  280  A/cm".  Figure  I  shows  the 
square  modulus  of  the  response  function  [Eq.  (33)]  for  the 
ideal  case  (no  carrier  heating)  and  in  Fig.  2  for  different 


FIG.  I.  Ideal  frequency  response.  The  current  density  ( 10*  A  cm  and 
light  intensity  (MW  cm'*)  respectively  for  each  curve  are.  in  order  of 
peak  frequency:  (1)  0.J75.  1.19;  (2)  0.72J.  2.39;  (3)  1.7a  3.74;  (4)  3.62, 
12.3. 
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FIG.  2.  Frequency  response  with  i£,=0.30  eV  The  earner  lemperaiure 
(K).  current  density  (Iff*  Acm"),  and  light  intensity  (MW  cm**), 
respectively,  are  in  order  of  peak  frequency:  ( 1 )  310,  0.t>48.  2.12;  (2)  325. 
1.70,  5.74;  (3)  350.  3.62.  12.3;  (4  )  400.  8.11.  2'.7. 

Steady-State  carrier  temperatures  when  the  lattice  temper¬ 
ature  is  300  K  and  A£j=300  meV.  Comparison  with  ob¬ 
served  responses  in  a  100  fxm  cavity  length.  5  well  laser  in 
Fig.  3”  shows  that  the  model  qualitatively  describes  the 
trends  reasonably  well.  (We  discuss  the  quantitative  as¬ 
pects  later.)  Figure  shows  the  square  of  the  peak  fre¬ 
quency  (w^)  as  a  function  of  photon  density  for  different 
levels  of  injected  energy.  From  this  plot,  we  can  see  that  as 
the  injection  energy  rises,  the  peak  frequency  is  increas¬ 
ingly  suppressed  at  high  photon  densities.  Modeling  these 
responses  according  to  the  simple  nonlinear  gain  model 
described  above,  we  obtain  nonlinear  gain  coefficients  of 
3.5  X  10-  and  8.3  X  10" ' '  cm^  for  A£,=0.3  and  0.75  eV, 
respectively.  Withou  earner  heating  these  responses 
would  be  a  straight  hne  as  demonstrated  by  the  ideal  per¬ 
formance.  In  fact,  at  low  photon  densities,  a  straight  line  is 
obtained  at  low  current  densities  for  all  injected  energies, 
which  goes  some  way  in  explaining  the  linear  behavior 
observed  by  Lester  er  ai We  should  further  point  out  that 
the  relaxation  frequency  of  £q.  (34)  is  not  a  good  guide  to 
what  frequencies  are  attainable.  Its  variadon  with  photon 
density  is  vinually  independent  of  the  injected  energy,  be¬ 
ing  principally  determined  by  the  characteristic  time  con¬ 
stants  of  the  system. 


FIG.  3.  Observed  frequency  response  of  s  lOOx  5  five  quantum  well 
strained  layer  laser  (see  Ref.  19). 


FIG.  4.  Squared  angular  frequencies  as  a  function  of  photon  density  ( 10” 
cm*')-  The  curves  are  (from  top  down)  (1)  ideal.  (2)  A£,=0.30  eV. 
and  (3)  A£,=0.75  eV. 

Figure  5  shows  the  respo..se  when  5l£j=0,  i.e.,  zero 
energy  input.  In  this  case  the  steady-state  temperature 
suys  at  300  K  for  all  input  current  densities,  but  recom¬ 
bination  heating  and  hot-phonon  effects  are  still  present. 
Note  that  the  sag  at  low  frequencies  has  been  eliminated. 
Figure  6  shows  the  response  for  il£y=0.75  eV,  corre¬ 
sponding  to  higher  barriers,  which  turns  out  to  be  similar 
to  that  for  A£j=0,  though  the  low-frequency  sag  has  re¬ 
turned.  but  in  a  much  reduced  form.  .A.t  first  sight  this  is  a 
surpnsing  result — one  expects  that  pumping  m  more  en¬ 
ergy  would  make  things  worse — but  this  is  not  the  case.  In 
fact  the  situation  with  A£j=0.3  eV  is  worse  than  either  of 
the  others.  The  cause  of  this  peculiar  behavior  is  to  be 
found  in  the  factor  (1-rr—p)  discussed  earlier.  Basically, 
the  response  is  damped  by  the  processes  of  photon  loss, 
phonon  loss,  and  recombination  heating,  all  of  which  con¬ 
tribute  an  effective  time  constant  of  about  10  ps.  and  this 
determines  the  magnitude  of  the  factor  c  in  the  denomina¬ 
tor  of  /{(tn).  This  damping  is  offset  to  some  degree  by  the 
factors  a  and  b  in  the  numerator.  W'here  c  is  only  weakly 
dependent  on  via  the  factor  /,  this  is  not  so  for  a  and 
b,  both  of  which  are  heavily  dependent  on  the  factor  ( 1  -f-r 
—p).  When  the  latter  is  negative,  the  frequency  depen¬ 
dence  of  the  numerator  under  some  circumstances  can  be 
small,  and  so  the  damping  effects  are  maximum.  This  cu- 


FTC.  5.  Frequency  response  with  a£,w0.  Current  densities  and  light 
intensities  as  for  the  ideal  case  (Fig.  I ). 
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Frequency  (GHz) 

FIG.  4.  Frequency  response  with  A£,«0.75  eV.  The  earner  tetnperattOT 
,JC;,  current  density  '10*  A  cm'*),  and  lisht  intenaty  (MW  cm'*), 
resoectively.  and  ate  in  order  of  peak  frequency;  ( 1 )  322.  0.373.  1.08:  ;2) 
350.  0.723.  116;  :3)  aOO.  1.63.  5.69;  (4)  450.  166.  9.27. 

nous  behavior  redeca  the  competitioa  between  recombina¬ 
tion  cooling  and  injection  heating  in  a  degenerate  gas. 

Because  of  the  simpiihcadon  and  approximaaons 
adopted  in  the  model  and  in  its  application  to  the  GaAs/ 
Gao.7lno jAs  system,  it  is  interesting  to  see  the  e£ect  of 
varying  some  of  the  magnitudes  which  enter  the  response. 
Figure  '  shows  the  result  of  doubling  the  scattering  time 
constant  and  the  phonon  lifetimes  for  the  case  when 
=0.3  eV.  Although  some  reduenon  in  frequency  occurs 


FIG  '  Frequency  raponsc  for  a£, >0.30  cV  with  longer  duM  cousunu 
and  a  lower  differenuai  gain,  (ai  Scuteruif  ume  rw340  h  and  r,w8  pa. 
‘  6 1  Vs  a  I  but  in  addiuon.  a  lialved  differcnual  gam.  The  dashed  cu.wes 
rcpieieiii  responses  ai  :hc  tane  curreni  density,  but  with  the  ume  con- 
xnnts  and  differential  gam  as  in  Fig.  1 
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the  mam  cdec:  is  to  enhance  the  low  frequency  droop  and 
generally  to  damp  the  response.  .Many  factors  can  reduce 
the  didierential  gam  and  Fig.  7  also  illustrates  the  effect  of 
reducing  by  a  factor  of  2  and  mamtainmg  the  longer 
time  constants.  The  frequency  at  peak  response  is  reduced 
as  c.spected  by  a  factor  of  roughly  v2.  and  the  droop  is 
greatly  enhanced,  giving  response  curves  much  nearer  to 
those  observed.” 

V.  DISCUSSION 

The  model  chosen  by  as  to  illustrate  the  effects  of  ear¬ 
ner  heating  contains  many  simplifying  features.  These  in¬ 
clude  equal  electron  and  hole  masses,  parabolic  bands,  in- 
liniteiy  deep  wells,  instantaneous  thermalization  within  the 
gas  of  earners,  the  applicability  of  degenerate  statistics, 
and  the  neglect  of  level  broadening  and  lattice  heating. 
Close  quantitative  agreement  with  published  work  was  not 
obtained  and  was  not  e.tpected.  but  we  hopw  we  have 
shown  that  earner  heatmg  produces  complex  effects  that 
simulate  the  main  elements  of  nonlinear  gam.  although  it  is 
pomted  out  that  these  effects  cannot  be  quantihed  in  terms 
of  a  smgle  phenomenological  factor,  e^.  Thus,  the  rela.x- 
ation  frequency  observed  in  a  pulse  e.xpenment  is  affected 
differently  by  nonlineanty  than  is  rhe  frequency  of  maxi¬ 
mum  response.  .Moreover,  the  droop  in  response  at  low 
frequencies  caused  by  injection  heating  :s  not  described  by 
the  usual  phenomenological  expression  for  nonlinearity. 
Tnis  IS  not  surprising  in  view  of  the  intricate  interplay  of 
the  dynamics  of  injection  heatmg.  recombination  heating, 
and  phonon  heating.  Interesungly,  the  frequency  response 
remains  useful  out  towards  1(X)  GHz.  which  suggests  that 
earner  heating  itself  need  not  be  a  limitation.  This  is  im- 
ponant  since  the  use  of  quantum  weils  entails  the  fact  of 
injection  heatmg,  although  as  we  have  illustrated  m  Fig.  5. 
even  without  injection  heating  che  process  of  recombina¬ 
tion  heating  with  hot  phonons  is  present  and  mtroduces 
■‘nonlinear  gam"  effects.  Nevertheless,  the  amount  of  in¬ 
jection  heatmg  is  a  function  of  the  matenal  structure  of  the 
laser  and  is  therefore  m.'mipulabie  in  pnnaple. 

In  our  illustration  ’.he  differential  gam  {g,)  was  be¬ 
tween  2  and  3x  10“”  cm  “7  Although  differential  gams  of 
this  magnitude  have  been  reponed.  the  more  usual  figures 
found  are  a  factor  of  2  or  3  smaller.  Lower  gains  can  easily 
come  about  through  the  comoined  effect  of  nonparabolic- 
ity,  lifetime  broadening,  Snite  well  depth,  well  width  fluc¬ 
tuations.  and  hole  burning.  Another  important  variable  is 
the  injection  efficiency  (tjj.  There  can  be  significant  loss  of 
comers  in  ndge  lasers  by  sideways  diffusion  and  in  mesa 
lasers  by  surface  recombination.  Probably  v  in  small  lasers 
does  not  exceed  0.30  and  may  be  as  low  as  0.5. 

Another  sunplifying  assumption  is  that  of  instanta¬ 
neous  thermalization.  As  long  as  the  electron  and  hole 
masses  are  about  equal,  this  assumption  is  not  likely  to  lead 
to  major  error  as  far  as  mtrasubband  processes  are  con¬ 
cerned.  However,  the  large  momentum  change  involved  in 
•ntersubband  transitions — including  the  process  of  well 
capture — may  introduce  further  time  constants  of  the  or¬ 
der  of  1-10  ps.  The  population  in  each  subband  would  then 
have  to  be  considered  individually,  leading  to  a  signifi- 
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cantly  enhanceti  compiexity.  This  aspec:  needs  further  in¬ 
vestigation.  but  It  IS  unlikely  that  the  elective  relaxation 
tune  constant — m  our  present  model,  of  magnitude  20  ps — 
would  be  increased  beyond  20  ps.  Figures  7(a)  and  7(b) 
illustrate  the  degradauon  m  the  frequency  response  when 
time  constants  are  long  and  the  differential  gam  is  weak  for 
the  case  m  which  A£j=0.c.  As  can  be  seen  in  Fig.  7(a). 
longer  time  constants  do  not  produce  a  shift  in  the  reso¬ 
nance  peak,  but  do  significantly  e.tacerbate  the  sag  in  the 
modulation  response.  Figure  shows  the  e.tpected  re¬ 
sult  that  the  resonance  peak  changes  with  the  square  root 
of  the  differential  gam. 

There  are  several  factors  like  those  mentioned  above 
that  make  realistic  modeling  difficult,  but  there  is  one  fac¬ 
tor  we  have  yet  to  mention  which  is  of  prime  importance  to 
the  achievement  of  high-frequency  performance,  and  that 
IS  the  photon  flu.x  intensity  which  the  laser  can  support 
without  being  damaged.  Because  we  have  wished  to  illus¬ 
trate  earner  heating  effects,  we  have  chosen  high  current 
densities  with  corresponding  photon  flux  intensities  m  the 
range  1-30  MW  cm~*,  but  it  is  rare  to  find  real  lasers  with 
damage  thresholds  much  a'oove  10  MW  cm~’.  It  is  likely 
that  this  damage  threshold,  rather  than  earner  (or  even 
lattice)  heating,  which  will  determine  tiie  bandwidth  of  the 
laser. 

In  summary,  we  have  presented  a  model  of  carrier 
heatmg.  that  incorporates  recombination  heatmg.  injection 
heatmg,  and  hot-phonon  effects,  and  have  applied  the 
model  to  the  case  of  a  quantum  well  laser  and  its  frequency 
response.  Carrier  heatmg  was  shown  to  lead  to  nonlinear 
gains  and  to  a  modified  frequency  response.  Application  to 
the  GaAs/TnCaAs  strained  layer  laser  showed  that  a  char¬ 
acteristic  relaxation  time  constant  of  order  10  ps  deter- 
mmed  the  behavior,  and  this  time  constant  was  derived 
from  the  phonon  lifetime  and  the  temperature  relaxation 
time  m  the  degenerate  population  in  addition  to  the  more 
familiar  photon  lifetime.  It  was  pointed  out  that  a  distinc¬ 
tion  between  the  frequency  of  relaxation  oscillations  and 
the  frequency  of  maximum  response  must  be  drawn,  and 
that  the  various  nonlinear  effects  could  not  be  quantita¬ 


tively  described  satistactonly  with  a  single  nonlinear  gain 
parameter.  It  was  shown  that  the  combination  of  weak 
gam  and  long  relaxation  umes  can  have  a  devastating  effect 
on  the  frequency  response. 
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